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SESSION  1 

Modeling  of  Ultrashort-Pulse  Interaction 
with  Matter 


Ultrashort  Pulse  Lasers: 

A  New  Tool  for  Biomedical  Applications 


Joseph  Neev 

950  Acapulco  St.  Laguna  Beach,  CA.  92651 

INTRODUCTION 

The  past  decade  has  brought  about  significant  advances  in  Ultrashort 
pulse  lasers  technology.  The  development  of  broad-band  solid  state 
gain  media  opened  up  new  possibilities  for  ultrashort  pulse 
generation.  In  particular,  the  development  of  all-solid-state 
ultrashort  pulse  devices  promise  to  make  such  devices  rugged  and 
reduce  their  cost. 

Capitalizing  on  the  evolving  technology  of  ultrashort  pulse  lasers 
could  result  in  many  advantages  for  biomedical  applications.  Tissue 
interaction  characteristics  that  are  superior  to  conventional 
surgical  technologies  and  other,  longer  pulse  laser  systems  can  be 
offered.  The  major  advantages  of  the  ultrashort  pulse  laser  (USPL) 
tissue  ablation  method  are:  1)  efficient  ablation  due  to  the  small 
input  of  laser  energy  per  ablated  volume  of  tissue  and  the  resulting 
decrease  of  energy  density  needed  to  ablate  material;  2)  minimal 
collateral  mechanical  damage  due  to  the  efficient  ablation  and  the 
short  duration  of  the  stress  impulse;  3)  minimal  collateral  thermal 
damage  due  to  the  extremely  short  deposition  time  and;  4)  the 
ablation  threshold  and  rate  are  less  dependent  on  tissue  type  and 
condition;  5)  high  precision  in  ablation  depth  is  achievable 
because  only  a  small  amount  of  tissue  is  ablated  per  pulse;  6)  low 
acoustical  (operating)  noise  level  (as  compared  to  the  acoustical 
noise  produced  by  other  laser  systems ) ;  7 )  minimized  pain  due  to 
localization  of  energy  deposition  and  damage;  8)  ability  to  texture 
surface  by  controlled  beam  profile  and  rastering;  9)  precise  spatial 
control:  the  intensity-dependent,  multiphoton  process  self-ensures 
that  tissue  below  or  laterally  removed  from  the  beam  focus  will  not 
experience  ablative  interaction,  and  finally,  10)  since  ultrashort 
pulses  interact  strongly  with  all  matter  regardless  of  specific 
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linear  absorption  characteristics,  efficient  processing  of  many 
tissue  types  is  possible.  Finally,  the  extremely  short  pulse 
duration  allows  ultrashort  pulse  lasers  to  become  useful  in  a  range 
of  biomedical  applications  involving  diagnostics  and  biomedical 
imaging . 


Current  Status  of  Ultrashort  Pulse  Laser  Development 

Titanium-doped  sapphire  (TitSapphire)  is  the  most  successful  laser 
medium  used  for  ultrashort  pulse  application  because  it  possesses  a 
broad  gain  bandwidth  of  approximately  200nm  along  with  favorable 
mechnical  and  thermal  properties. 

The  Generation  of  Ultrashort  Pulses  is  Currently  Most  Commonly 
Achieved  Using  Chirped  Pulse  Amplification  (CPA)  of  Ultrashort 
pulses.  Devices  utilizing  CPA  allow  increase  in  peak  power  to  the 
Hundreds  of  Terawatt  (TW)  level  and  result  in  decrease  of  pulse 
duration  down  to  <  10  fs.  Diode-pumped  fiber  oscillator  allow 
further  reduction  in  size  cost. 

Types  of  CPAs 

Several  types  of  Chirped  pulse  oscillator  are  now  available 
commercially : 

Nd: glass  amplify  Pulses  to  the  1  ps  range. 

Ti: Sapphire  or  CriLiSAF  result  in  amplification  of  pulses  as  short 
as  lOOfs  range 

and  Ti: Sapphire  -  High  average  power  oscillators  allow 
amplifications  to  the  10  to  20  fs  range 

Engineering  and  Design  Considerations 

Engineering  considerations  should  emphasize  the  following  operating 
parameters : 

Stable,  Consistent  operation 
Compactness  and  Portability 
Durability  and  ruggedness 

•  Delivery 

•  Control 


Commercially  Available  System  include,  for  excimple  - 
A  mode-locked,  diode-pumped  erbium  fiber  laser  repetition  rate  of 
37MHz  and  average  powers  of  50mW,  (InJ/pulse)  Turn-key,  self-start 
mode-locking;  High  power  ultrashort  pulses  at  1550nm  and  775nm; 

Fiber  gain-guided  cavity  for  ultrastable  operation  with  superb 
reliability;  Wall-plug  or  battery  powered  with  no  water;  and  a 
highly  compact  head  (12"  x  24") 

Presently,  an  amplified  Ti: Sapphire  Laser  System  is  characterized  by 

•  Energy  of  mJ/pulse 

•  KHz  repetition  rates 

•  Pulse  widths  ranging  from  sub-50fs  to  over  lOps 

•  High  contrast  ratio  and  near  bandwidth  limited 

•  Wavelengths  extending  over  the  entire  Ti: Sapphire  tuning  range 
Alternatively,  Ti: Sapphire  laser  systems  operating  at  energy  levels 
of  up  to  50mJ  at  10-20  Hz  are  also  available. 

APPLICATIONS 

Scientific  and  non-biomedical  Applications: 

•  Ultrafast  Pulses  allow  the  study  of  physical  phenomena  that 
occur  in  extremely  short  time  scales  (picosecond  and  femtoseconds) 
and  the  study  of  high-speed  dynamics 

•  Snapshots  in  this  time  domain  reveal  the  most  fundamental 
mechanisms  of  molecular,  atomic,  and  electron  interactions.  Sub  100 
fs  pulses  can  probe  the  rotational  and  vibrational  dynamics  of 
atomic  systems  such  as  molecules  and  condensed  mater. 

Monitoring  of  molecular  motion  in  liquids  or  gases,  and 
studying  of  how  electrons  collide  in  semiconductors  and 
superconductors,  electrons,  excitons,  and  phonons  in  various 
condensed-matter  media 

•  Quantum  Control  of  Nal  Photo-dissociation  Reaction  Product 
States  by  Ultrafast  Tailored  Light  Pulses 

Studies  of  high  intensity  field  phenomena 
XUV  Pulse  generation 

Ultrafast  incoherent  X-ray  generation 
High  order  harmonic  generation 


Additional  applications  were  provided  by  the  ultrashort  pulse 
generation  of  Terahertz  Electromagnetic  Pulse  Generation  for 
Mil lime ter- Wave  Spectroscopy 

Terahertz  radiation  allows: 

differentiation  between  various  materials,  chemical 
compositions,  or  environments 

•  determination  of  fat  Content,  water  content  in  vegetation,  and 
internal  content  of  containers - 

•  biomedical  imaging  of  tissue 

•  chemical-reaction  analysis,  environmental  and  pollution 
control,  materials  inspection,  fault  Detection 

•  profiling  of  doping  and  defects  in  semiconductors,  and 
packaging  inspection 

Additional  biomedical  applications  include: 

Optical  coherence  tomography 

Optical  coherence  tomography  (OCT)  is  an  optical  imaging  technique 
that  uses  low  coherence  interferometry  to  obtain  micron  scale, 
cross-sectional  images  of  biological  systems.  OCT  has  been  used  to 
provide  topographic  images  of  the  transparent  structures  in  the  eye. 
Clinical  studies  have  shown  that  OCT  provides  high  resolution, 
cross-sectional  images  of  the  retina  and  can  be  used  to  diagnose  a 
wide  range  of  retinal  macular  diseases.  OCT  imaging  in  other  human 
tissues  is  more  difficult  due  to  optical  scattering.  However, 
recent  in  vivo  studies  have  shown  that  OCT  can  image  architectural 
morphology  in  highly  optically  scattering  tissues.  Performing 
optical  biopsy  with  OCT  requires  powerful  sources  with  broad 
frequency  bandwidth.  Solid  state  lasers  based  on  Ti:Al203  and 
Cr:Mg2Si04  and  rare-earth  doped  fiber  lasers  are  capable  of 
providing  hundreds  of  milliwatts  of  single-trans verse  mode  light 
with  coherence  lengths  as  short  as  a  few  microns. 
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Cornea  Surgery 

Corneal  PhOto-disruption  with  ultrashort  pulse  lasers  utilizes 
laser-induced  optical  breakdown  (LIOB)  which  take  place  with  the 
pulse  energy  density  reaches  a  plasma  formation  threshold.  The 
reduce  thermal  and  mechnical  collateral  effects  associated  with 
ultrashort  pulse  interactions  offer  superior  performance  for  cornea 
reshaping  as  well  as  for  intrastromal  surgery. 

Related  topics  include; 

•  Laser-induced  breakdown  in  the  eye  at  pulse  durations  from  100  ns 
to  100  fs 

•  Retinal  damage  mechanisms  from  ultrashort  laser  exposure 

•  Refractive  surgical  applications  with  ultrashort  pulsed  lasers, 

Dental  Hard  Tissue  Removal 

While  several  systems  are  currently  being  considered  for  the  purpose 
of  processing  of  hard  dental  tissue,  undesired  thermal  collateral 
damage  and  the  lack  of  speed  and  efficiency  are  major  obstacles. 
Picosecond  and  femtosecond  pulse  duration  have  shown  themselves  able 
to  produce  very  precise  cavities  with  only  negligible  thermal  or 
shock  wave  collateral  damage.  Surface  quality  and  morphological 
characteristics  also  appear  to  be  superior  to  those  of  other  laser 
systems.  If  a  proper  deliver  system  can  be  identified  and  developed 
and  if  the  cost  of  ultrashort  pulse  laser  system  will  continue  to 
drop,  such  devices  may  become  valid  alternative  to  existing 
technology . 

Additional  topics  include 

•  Microsurgical  Effect  of  laser  irradiation  on  the  functional 
activity  of  enzymes  with  different  structural  complexity 

•  Influence  of  ultrashort -pulse  laser  on  gandioblockade  compounds, 

•  Feedback  system  for  precision  ultrashort  pulse  biotissue  ablation 

•  Light  transmission  imaging  in  tissue-like  phantom 

•  Real-time  2-photon  confocal  microscopy 


AREAS  OF  ACTIVE  RESEARCH  LEADING  TO  ENABLING  ULTRA  SHORT  PULSE 
LASER  TECHNOLOGY  INCLUDE: 

•  Ultrashort  Pulse  Lasers :  Instrumentation  and  Measurements 
Compact  ultrafast  sources  for  biomedical  applications 

•  Ultrashort  laser  pulse  propagation  through  hollow  core  fibers 

•  Review  of  ultrashort  pulse  measurement:  changing  the  basic 
ultrashort  pulse 

•  Feedback  system  for  precision  ultrashort  pulse  biotissue 
ablation 

•  Modeling  of  Ultrashort  Pulse  Interaction  with  Matter 
Ultrashort  Pulse  Propagation  in  Transparent  media 

•  Radially  varying  dispersion  in  high-numerical-aperture  focusing 

•  Pressure  and  temperature  evolution  induced  by  ultra-short  laser 
pulse  ablation. 
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Radially  varying  dispersion  in  h^h-numerical-aperture  focusing 
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®BioCentrum  Amsterdam 

Department  of  Molecular  Cytology,  University  of  Amsterdam. 
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^Electrical  and  computer  engineering  department 
University  of  California  San  Diego, 
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ABSTRACT 

Over  the  last  few  years  a  number  of  microscopical  techniques  have  been  developed  that  take  advantage 
of  ultrashort  optical  pulses.  All  these  techniques  rely  on  temporal  pulse  integrity  at  the  focal  point  of  a 
high-numerical  aperture  (NA)  focusing  system.  We  have  investigated  the  dispersion  induced 
broadening  for  pulses  on  the  optical  axis,  using  the  two-photon  absorption  autocorrelation  (TPAA) 
technique.  We  demonstrate  that  the  induced  broadening  can  be  pre-compensated  for  by  a  properly 
designed  dispersion  pre-compensation  unit  for  pulses  as  short  as  15  femtoseconds. 

Another  source  of  pulse  broadening  in  high-NA  focusing  systems  is  due  to  radial  variations  in  the 
dispersion  over  the  pupil  of  the  objective.  This  may  cause  differences  in  the  group  delay  between  on- 
axis  and  outer  ray  wave  packets,  as  well  as  differences  in  tire  broadening  of  the  wave  packets 
themselves.  In  this  paper  we  present  experimental  results  on  the  measurement  of  these  radial  variations 
in  the  dispersion  characteristics  over  the  aperture  of  high-NA  microscope  objectives,  using  a  slightly 
modified  TPAA  technique. 

Keywords:  dispersion,  group  delay,  multi-photon  microscopy 
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1.  INTRODUCTION 


The  advances  in  recent  years  in  ultrashort  laser  pulse  technology  -  notably  all-solid  state  tunable  lasers 
which  routinely  produce  sub-100  fs  optical  pulses  -  has  opened  the  way  to  the  application  of  non-linear 
optical  techniques  in  new  fields,  such  as  high  resolution  multi-photon  microscopy  [Denk,  et  al.,  1990; 
Webb,  1990;  Piston  &  Webb,  1991;  Hell,  et  al.,  1996;  Szmacinski,  et  al.,  1996].  Since  the 
fluorescence  intensity  in  two-  and  three-photon  absorption  depends  non-linearly  on  the  excitation 
intensity,  maximum  efficiency  is  attained  when  both  the  spatial  focusing  conditions  are  diffraction 
limited  and  the  optical  pulses  have  a  minimum  pulse  duration  at  focus.  This  is  especially  important  for 
microscopy  of  biological  objects  where  the  radiation  load  on  the  specimen  is  to  be  kept  at  a  minimum  to 
prevent  biological  damage,  while  at  the  same  time  extracting  the  maximum  of  information  in  terms  of 
high  signal-to-noise  microscopic  images. 

The  application  of  femtosecond  pulses  -  which  inherently  have  a  substantial  bandwidth  -  to 
microscopy  leads  to  a  new  regime  of  imaging  where  the  effects  of  (temporal)  dispersion  should  be 
taken  into  account.  In  general  this  dispersion  will  cause  broadening  and  even  distortion  of  the  optical 
pulses,  leading  in  turn  to  a  reduced  efficiency  in  multi-photon  absorption.  Especially  since  the 
microscope  objectives  consist  of  a  large  number  of,  highly  dispersive,  glass  elements,  the  effects  of 
dispersion  are  generally  non-negligible.  Previously  we  investigated  the  influence  of  on-axis  dispersion 
on  the  optical  pulse  duration  at  the  focal  point  of  high-numerical  aperture  (NA)  microscope  objectives 
and  we  developed  methods  to  measure  and  pre-compensate  the  induced  dispersion  [Brakenhoff,  et  al., 
1995;  Muller,  et  al.,  1995;  Muller,  et  al.,  1997].  In  this  paper  we  extend  the  analysis  of  dispersion 
induced  by  high-NA  microscope  objectives  to  include  dispersion  which  is  a  function  of  the  normalised 
radius  with  respect  to  the  entrance  pupil  of  the  objective. 


2.  DISPERSION  IN  HIGH-NA  OBJECTIVES 


Dispersion  of  femtosecond  optical  pulses  is  the  result  of  die  functional  dependence  of  the  refractive 
index  on  the  wavelength.  This  results  in  a  frequency  dependent  phase  shift  experienced  by  the  optical 
pulse.  To  separate  out  the  various  contributions  to  the  total  dispersion  we  expand  the  phase  ^(to)  in  a 

Taylor  series: 


(p(a))  =  (p{(Oo)  + 


(  d(p\ 
\do}) , 


(CO  -  COg)  + 


l(^) 

2l[da}V 


(co -a)g)^+... 


^  (Po  +  tp,  +  (P2  +  ... 


(1) 
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Here  cOg  denotes  the  centre  frequency  of  the  pulse.  q)g  is  merely  a  phase  offset  and  q)^  is  the  group 
delay  (GD)  experienced  by  the  pulse.  The  higher  order  terms  -  9)2  •  group  delay  dispersion  (GDD),  : 
third-order  dispersion  (TOD),  etc.  -  lead  to  a  broadening  of  the  pulse  and  to  a  distortion  of  the  pulse 
shape.  The  total  amount  of  dispersion  induced  on  a  pulse  propagating  through  a  material  depends  both 
on  the  material  properties  (through  n{(o)),  the  pathiength  and  the  frequency  bandwidth  of  the  pulse. 

In  the  case  of  dispersion  induced  by  high-NA  microscope  objectives  we  should  discriminate 
between  two  possible  sources  of  dispersion.  The  first  is  the  ordinary  dispersion  resulting  from  the 
propagation  of  the  pulse  through  the  various  glass  elements  in  the  objective.  We  will  refer  to  this  as  the 
on-axis  dispersion.  The  second  source  for  possible  temporal  broadening  and  distortion  of  the  pulse  at 
the  focal  point  is  dispersion  which  is  a  function  of  the  relative  position  of  the  line  of  propagation  with 
respect  to  the  objective’s  entrance  pupil.  In  other  words,  differences  in  dispersion  between  rays 
propagating  through  the  centre  of  the  lens  or  through  the  outer  parts.  We  will  refer  to  this  as  radially 
varying  dispersion.  Clearly  for  the  on-axis  dispersion  only  the  second  and  higher  order  dispersion 
terms  are  important  since  they  will  cause  a  general  broadening  and  distortion  of  the  pulse.  For  the 
radially  varying  dispersion  however,  we  also  need  to  consider  differences  in  the  GD  experienced 
between  on-axis  and  off-axis  rays  (i.e.  GD(r)),  as  well  as  radially  varying  higher  orders  of  dispersion 
GDD{r),  TOD{r),  etc.. 


3.  MEASURING  THE  DISPERSION  AT  THE  FOCAL  POINT  OF  A  HIGH-NA  LENS 

In  the  measurement  of  the  dispersion  induced  by  high-NA  microscope  objectives  we  discriminate 
between:  on-axis  dispersion  and  radially  varying  dispersion.  Two  techniques  have  been  developed  to 
measure  each  of  these  separately. 

3.1  Two-photon  absorption  autocorrelation  (TPAA) 

Since  the  TPAA  technique  has  been  discussed  in  detail  elsewhere  [Brakenhoff,  et  al.,  1995;  Muller,  et 
al.,  1995;  Muller,  et  al.,  1997],  we  summarise  the  main  characteristics  here.  The  general  scheme  of  the 
measurement  is  depicted  in  figure  1.  The  Iciser  beam  first  passes  a  dispersion  pre-compensation  unit, 
based  on  the  double  prism  pair  compensation  configuration,  permitting  the  addition  of  a  variable 
amount  of  dispersion  -  from  positive  to  negative  -  to  the  pulses.  The  pulses  are  then  split  in  two  parts, 
one  passing  a  variable  delay  line  before  being  recombined  with  the  other  part.  The  recombined  parts 
then  propagate  collinearly  either  to  a  conventional  second  harmonic  generation  (SHG)  set-up  for  a 
reference  pulse  duration  measurement,  or  are  focused  by  the  objective  under  investigation  in  a  dye 
solution  (10-3  M  Rhodamine  6G  in  water)  to  generate  two-photon  absorption  (TPA).  The  fluorescence 


is  detected  in  the  backscattering  direction.  The  autocorrelation  signal  consists  of  the  measurement  of 
either  the  SHG  signal  or  the  TPA  fluorescence  signal  as  a  function  of  the  time  delay,  induced  by  the 
variable  pathlength  in  the  variable  delay,  between  the  pulses. 


Figure  1.  General  schematic  of  the  TPAA  set-up.  Symbols  used:  DPC:  dispersion  pre¬ 
compensation  unit;  VD:  variable  delay;  D:  dichroic  mirror;  BBO:  type  I  BBO  doubling 
crystal;  PMT:  photo  multiplier  tube. 

The  procedure  for  the  TPAA  measurements  is  then  as  follows.  First,  the  system  is  optimised  to 
produce  the  minimal  pulse  width  in  the  SHG  autocorrelation.  For  the  laser  system  used  in  these 
experiments  this  resulted  in  near  transform  limited  (Au  At  ~.37,  assuming  a  Sech^  input  pulse)  pulses 
of  15.5  fs.  This  provides  the  “zero-setting”  of  the  disjjersion  pre-compensation  unit.  Next,  the  beam  is 
directed  through  the  microscope  objective  to  be  investigated.  Again  the  dispersion  pre-compensation 
unit  is  optimised  to  provide  the  minimal  pulse  duration  in  the  TPA  autocorrelation.  Note  that  this 
corresponds  to  a  maximum  in  the  TPA  fluorescence.  From  the  adjustment  required  in  the  dispersion 
pre-compensation  unit  the  induced  dispersion  by  the  objective  can  be  calculated. 

A  typical  result  of  the  TPAA  signal  obtained  in  this  manner  is  shown  in  figure  2.  The  objective  under 
study  in  this  case  was  a  Zeiss  CP-Apochromat  l(X)x/l,25  oil.  The  configuration  of  the  dispersion  pre¬ 
compensation  set-up  corresponded  to  an  induced  dispersion  of  GDD  =  -579  fs^,  TOD  =  -209  fs^  and 
FOD  =  -156  fs4  A  dispersive  ray-tracing  calculation  of  the  same  system  -  based  on  the  optical  data 
provided  for  this  objective  -  gives:  GDD  =  +579.4  fs^,  TOD  =  +129.9  fs^  and  FOD  =  +24.3  fs^, 
demonstrating  excellent  agreement  between  die  experimental  results  and  the  theoretical  modeling.  Note 
that  we  ensured  that  the  entrance  pupil  of  the  microscope  objective  was  significantly  under-filled  to 
eliminate  any  possible  effects  from  radially  varying  dispersion. 
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Figure  2.  TPAA  signal  for  a  dispersion  pre-compensation  setting  (579  fs^)  giving  minimal 
pulse  width  for  the  Zeiss  CP-Apochromat  lOOx/1,25  oil  microscope  objective.  The  best  fit 
was  obtained  assuming  a  Sech^  input  pulse. 

Similar  measurements  have  been  done  for  a  number  of  microscope  objectives,  the  results  of  which  are 
summarised  in  table  I.  These  results  clearly  show  that  it  is  possible,  using  a  properly  designed 
dispersion  pre-compensation  unit,  to  focus  15  fs  optical  pulses  with  high-NA  microscope  objectives.  It 
is  important  to  note  that  the  slight  increase  in  pulse  width  for  the  higher  GDD  pre-compensation  values 
results  from  residual  third-order  dispersion  from  the  pre-compensation  unit.  This  shows  that  it  is 
possible  -  using  a  specizdly  designed  dispersion  pre-compensation  set-up  -  to  obtain  the  original  pulse 
width  at  the  focal  point  of  any  of  these  microscope  objectives. 


Table  1 

Objective 

Pulse  width 

(fs) 

GDD 

(fs2) 

SHG 

15.5 

0 

Zeiss  C-Apochromat  63x/l,2W  Korr 

18.4 

-1140 

Zeiss  C-Apochromat  40x/l,2W  Korr 

18.4 

-1104 

Zeiss  Plan  Neofluar  lOOx/13  oil 

14.7 

-778 

Zeiss  Plan  Neofluar  63x/l,25  oil 

16.4 

-887 

Zeiss  Plan  Neofluar  40x/l,3  oil 

17.5 

-1104 

Zeiss  CP-Achromat  lOOx/1,3  oil 

15.9 

-579 

Table  /.  The  minimum  pulse  width  obtained  in  the  TPAA  measurements  for  a  series  of 
microscope  objectives  and  the  dispersion  pre-compensation  value  required  to  obtain  this 
pulse  width. 
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3.2  Two-photon  absorption  cross-correlation  (TPAC) 

Since  the  pathlengths  through  the  various  glass  elements  in  a  microscopic  objective  are  different  for  the 
chief  and  marginal  rays,  we  need  to  consider  dispersion  variations  eis  a  function  of  the  radial  position 
with  respect  to  the  entrance  pupil.  To  be  able  to  measure  this,  the  general  TPAA  set-up  is  modified  to 
cross-correlate  the  outer  rays  with  the  central  rays.  Figure  3  shows  a  schematic  of  the  modifications. 


Obj. 


Figure  3.  Schematic  of  the  modifications  to  the  TPAA  set-up  (figure  1)  to  enable  TPAC.  (a) 
In  the  “reference”  arm  of  the  autocorrelator  an  iris  is  placed  j)assing  only  the  centre  of  the 
beam.  In  the  “object”  arm  of  the  autocorrelator  an  annulus  of  variable  radius  (r)  is  placed, 
(b)  When  focused  in  a  non-linear  medium  a  TPA  cross-correlation  signal  is  produced  as  a 
function  of  the  relative  delay  between  object  and  reference  pulse. 


In  one  of  the  beams  in  the  autocorrelator  -  to  be  denoted  by  “reference”  beam  -  an  iris  is  positioned, 
passing  only  the  centre  of  the  beam.  In  the  other  beam  -  denoted  by  “object”  -  an  annulus  of  variable 
radius  is  positioned.  The  TPA  signal  will  thus  consist  of  a  cross-correlation  between  the  central  rays 
and  the  rays  passed  by  the  annulus.  The  relative  time  position  of  the  maximum  of  this  cross  correlation 
signal  determines  the  relative  delay  between  the  reference  and  object  rays.  The  delay  equals  zero  value 
is  taken  for  the  (auto)  cross-correlation  with  two  irises  in  both  arms.  The  principle  of  the  TPAC 
measurement  is  to  measure  the  relative  delay  value  for  the  maximum  of  the  cross-correlation  signal  as  a 
function  of  the  radius  of  the  aimulus  in  the  object  beam. 

A  typical  result  of  such  a  measurement  is  shown  in  figure  4  for  the  Zeiss  CP-Achromat  lOOx/13  oil 
microscope  objective.  The  central  rays  are  delayed  more  than  the  outer  rays.  In  case  of  the  15  fs  pulses 
used  in  these  experiments  this  effectively  means  that  there  is  no  interaction  between  pulses  propagating 
through  the  centre  of  the  lens  and  pulses  propagating  through  the  outer  part.  Clearly  this  temporal 
effect  is  coupled  to  the  spatial  domain  and  consequently  may  affect  the  point  spread  function  of  the 
system. 


Figure  4.  TPAC  signal  for  the  Zeiss  CP-Apochromat  lOOx/1,25  oil  microscope  objective, 
showing  an  increasing  group  delay  as  a  fimction  of  normalised  radius. 


Figure  5.  TPAA  signal  as  a  function  of  normalised  radius  at  a  dispersion  pre-compensation 
setting  giving  the  minimal  pulse  width  at  r/R=0  for  the  Zeiss  CP-Apochromat  l(X)x/l,25  oil 
microscope  objective.  The  best  fit  was  obtained  assuming  a  Sech^  input  pulse. 

The  methods  of  TPAA  and  TPAC  can  be  combined  to  measure  the  pulse  broadening  as  a  function  of 
radius.  In  this  case,  equal  annuli  are  placed  in  both  arms  of  the  autocorrelator  and  the  TPAA  signal  is 
measured  in  a  similar  manner  as  in  section  3. 1.  As  is  shown  in  figure  5  there  is  negligible  pulse 
broadening  for  this  objective  as  a  function  of  radius. 
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4.  CONCLUSION. 


Significant  amounts  of  dispersion  are  induced  by  microscope  objeetives  on  femtosecond  optical  pulses. 
There  are  two  contributions  to  this  dispersion;  (i)  an  on-axis  contribution;  and  (ii)  a  radially  varying 
contribution.  The  on-axis  dispersion  causes  a  substantial  broadening  of  the  optical  pulse,  which  can 
however  be  eompensated  for  by  a  properly  designed  dispersion  pre-compensation  unit.  It  has  been 
shown  that  it  is  possible  to  pre-compensate  the  dispersion  for  pulses  as  short  as  15  femtoseeonds.  The 
radially  varying  dispersion  is  primarily  due  to  a  radiedly  varying  group  delay.  This  causes  wave  paekets 
propagating  through  the  centre  of  the  lens  to  experience  a  larger  delay  than  wave  packets  propagating 
through  the  outer  part  of  the  objective.  This  effect,  which  is  independent  of  the  pulse  duration  itself, 
becomes  relatively  more  important  for  pulses  of  deereasing  pulse  duration.  The  effect  of  a  radially 
varying  group  delay  may  potentially  reduce  the  numerical  aperture  of  the  system  and  eonsequently 
degrade  the  resolution  in  microscopic  applications.  Experiments  are  under  way  to  determine  how  the 
resolution  is  affeeted.  We  have  shown  that  the  pulse  broadening  itself  -  due  to  ODD  -  is  approximately 
constant  over  the  radius  of  the  pupil  of  the  objective. 
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ABSTRACT 

Pulse  broadening  of  ultrashort  optical  pulses,  as  short  as  15  femtoseconds  (fs),  due  to  the  propagation  through  high- 
numerical-aperture  (NA)  microscope  objectives  can  be  pre-compensated  to  ensure  temporal  pulse  integrity  at  the  focal  point. 
The  predictions  from  dispersive  ray-tracing  calculations  show  excellent  agreement  with  the  experimental  results  from  two- 
photon  absorption  autocorrelation  (TPAA)  for  the  Zeiss  CP-Achromat  lOOx/1,25  oil  microscope  objective.  From  this, 
general  predictions  can  be  inferred  for  dispersion  in  most  types  of  microscope  objectives.  Key  element  to  the  work  is  a 
carefully  designed  dispersion  pre-compensation  configuration,  which  minimizes  pulse  broadening  due  to  residual  third  order 
dispersion.  The  capability  to  focus  these  ultrashort  pulses  with  control  of  the  pulse  definition  at  the  focal  point  is  important 
for  two-photon  absorption  (TP A)  and  time-resolved  microscopy. 

Keywords:  High  numerical  aperture  objectives,  ultrashort  pulses,  dispersive  ray-tracing,  two-photon  absorption 
autocorrelation 

1.  Measurement  and  modeling  the  propagation  of  ultrashort  pulses  in  a  high  NA  system 

Laser  systems  which  routinely  produce  sub-100  fs  optical  pulses  have  recently  become  widely  available.  These  systems 
have  inspired  the  application  of  ultrashort  optical  pulses  in  various  areas  of  microscopy,  including  two-photon  absorption 
(TPA)  microscopy  (e.g.;  [1,  2]),  multi-photon  absorption  imaging  (e.g.  [3,  4]),  fluorescence  lifetime  imaging  ([5,  6])  and 
optical  coherence  tomography  (e.g.  [7,  8]).  For  effective  application,  all  these  techniques  rely  on  proper  focusing  of 
femtosecond  optical  pulses  with  a  high  numerical  aperture  (NA)  microscope  objective.  Because  these  objectives  consist  of 
multiple  lens  elements,  they  induce  substantial  dispersion  to  -  inherently  broad  band  -  femtosecond  pulses.  Proper 
understanding  of  this  induced  dispersion  facilitates  the  application  of  dispersion  pre-compensation  techniques  to  ensure 
optimal  temporal  -  and  spatial  -  focusing  of  the  femtosecond  pulses. 

In  earlier  work  [9-12]  the  general  characteristics  of  focusing  >  50  fs  pulses  has  been  investigated.  It  has  been  found  that 
for  these  pulses  the  induced  broadening  is  primarily  caused  by  group  delay  dispersion.  The  effect  of  phase  time  delay  was 
found  non-negligible  only  for  lenses  with  substantial  chromatic  aberration  over  the  band  width  of  the  pulse  [13].  In  this 
paper  we  investigate  in  detail  the  induced  dispersion  when  focusing  15  fs  optical  pulses  with  a  Zeiss  CP-Achromat 
lOOx/1,25  oil  microscope  objective. 

The  induced  dispersion  by  the  NA  objective  is  measured  using  the  TPAA  technique,  which  has  been  described  in  detail 
elsewhere  [9].  Figure  1  shows  a  schematic  of  the  experimental  set-up.  The  output  from  a  home-built  Ti:Sappbire  oscillator, 
operating  at  800  nm,  80  MHz,  passes  a  dispersion  pre-compensation  unit  based  on  a  double  prism  set-up.  The  dihedral 
reflector  Ri  is  mounted  on  a  translation  stage  to  permit  varying  from  negative  to  positive  induced  dispersion.  The  fused 
silica  prisms,  with  an  apex  angle  of  69°,  are  used  at  55°20’  angle  of  incidence.  The  beam  then  passes  a  dispersion  balanced 
Mach-Zehnder  interferometric  autocorrelation  unit.  Insertion  of  mirror  Mj  directs  the  beam  to  an  off-axis  parabolic  mirror 
which  focuses  it  in  a  100  pm,  type  I KDP  crystal,  producing  a  second  harmonic  generation  (SHG)  signal,  which  is  defected 
through  a  BG39  filter  with  a  photomultiplier  tube  (PMT).  Alternatively,  the  mirror  Mj  can  be  removed  and  the  beam  is 

focused  by  a  high  NA  objective  (Zeiss  CP-Achromat  lOOx/1,25  oil)  into  a  10"^  M  solution  of  Rhodamine  6G  in  water.  The 
fluorescence  generated  in  the  induced  TPA  process  is  detected  in  the  backscattering  direction,  reflected  off  a  dichroic  mirror, 
through  a  BG39  filter  with  a  photomultiplier  tube  (PMT).  The  TPAA  amplitude  showed  a  quadratic  dependence  on  the  input 
pulse  power,  demonstrating  that  the  absorption  is  indeed  a  two-photon  process. 
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Figure  1. 

Schematic  of  TPAA  set-up.  (a)  DPC:  dispersion  pre-compensation  unit; 
BD:  beam  dump;  D:  dichroic  mirror;  O:  objective,  S:  sample;  C;  100  Jim 
type  I KDP  crystal;  PMT:  photomultiplier  tube.  Insertion  or  removal  of 
mirror  Mj  permits  switching  between  SHG  and  TPA  autocorrelation 
respectively,  (b)  Detail  of  DPC.  R2  is  a  dihedral  step-up  mirror  which  retro 
reflects  the  beam  at  an  elevated  plane  with  respect  to  the  incoming  beam. 
The  pick-off  mirror  Mp  reflects  the  beam  at  this  elevated  plane  out  of  the 


Figure  2b  shows  the  TPAA  signal  for  the 
dispersion  pre-compensation  unit  set  to 
minimum  pulse  width.  The  best  fit  to  data 
was  obtained  for  a  sech^  input  pulse,  yielding 
an  almost  transfoim  limited  (A\)At  ~.39) 
pulse  width  of  15.9  (+  0.1)  fs.  Note  that  all, 
noise  induced,  errors  in  the  pulse  width 
measurements  reported  here  are  based  on  the 
change  in  pulsewidth  measurement  from 
acquisition  to  acquisition.  Compared  to  the 
dispersion  setting  to  obtain  a  minimum  in  the 
SHG  autocorrelation  signal  of  15.5  (+  0.1)  fs 

(figure  2a),  an  additional  -579  fs^  of 
dispersion  was  required  to  compensate  for  the 
dispersion  induced  by  the  objective.  Figure  2c 
shows  the  dependence  of  the  pulse  width  and 
TPAA  amplitude  as  a  function  of  induced 
dispersion  pre-compensation.  As  expected,  the 
TPA  efficiency  increases  close  to  linearly  with 
decreasing  pulse  width.  The  fact  that  the 
minimum  TPAA  pulse  width  almost  equals 
the  minimum  pulse  width  obtained  with  SHG 
autocorrelation  indicates  the  absence  of  sample 
induced  broadening.  This  was  further 
substantiated  in  independent  experiments, 
where  the  KDP  crystal  was  substituted  for  a 
solution  of  Rhodamine  6G,  and  the  all 
reflective  optic  geometry  was  used. 


DPC.  The  choice  of  prism  material  is  critical 

for  optimal  compensation.  Traditionally, 
prism  pairs  have  been  made  from  dispersive  glasses  such  as  SFIO  in  order  to  produce  large  amounts  of  dispersion  with 
modest  prism  spacings.  The  difficulty  with  the  more  dispersive  glasses,  however,  is  that  the  rate  of  group  delay  dispersion 
(GDD)  to  third  order  dispersion  (TOD)  produced  by  the  prism  pair  is  a  poor  match  to  the  material  path  length  in  the 
objective  it  is  attempting  to  compensate.  Thus,  the  residual  TOD  due  to  the  prism  compressor  quickly  becomes  dominating, 
and  it  is  this  term  that  ultimately  limits  the  compressed  pulse  width.  Indeed,  attempting  to  compress  the  pulse  with  SFIO 
prisms  resulted  in  minimum  compressed  pulsewidths  no  better  than  30  fs.  For  this  reason,  fused  silica  prisms  are  used.  The 
prism  compressor  separation  increases  due  to  the  lower  dispersion  of  the  fused  silica,  but  the  residual  phase  errors  are 
smaller  enabling  the  pulse  to  be  recompressed  to  essentially  its  transform  limit. 


Using  the  optical  data  for  the  CP-Achromat  lOOx/1,25  oil  microscope  objective,  the  dispersion  characteristics  were 
theoretically  modeled.  The  geometric  optical  pathlength  propagated  by  a  certain  ray  is  calculated  using  straight  forward  ray¬ 
tracing.  The  dispersion  of  each  element  is  then  calculated  using  a  Taylor  expansion  of  the  frequency  dependent  phase 
((piO)))  [14]: 


<p(t»)=(p(fflo)  +  (-^l  (CO-O) 

By  expressing  the  derivatives  {d" (p j dG)’'^  in  terms  of  {d^nj dX'^,  multiplied  by  the  geometric  pathlength  through  each 

lens  element  in  the  objective,  the  total  dispersion  of  the  microscope  objective  can  be  calculated  to  any  desired  order.  The 
first  two  terms  in  equation  (1)  represent  a  phase  shift  and  the  group  delay  and  do  not  alter  the  shape  of  the  pulse.  The  higher 
order  terms  denoted  by  group  delay  dispersion  (GDD),  third  order  dispersion  (TOD)  and  fourth  order  dispersion  (FOD) 
respectively  cause  pulse  broadening  and  distortion.  A  dispersive  ray-tracing  calculation  for  the  on-axis  ray  of  the  objective 

gives  a  relative  contribution  of  these  terms  of:  GDD  =  579.4  fs^;  TOD  =  129.9  fs^  and  FOD  =  24.3  fs^.  Note  that  the 
factorial  preceding  each  coefficient  is  already  included  in  these  dispersion  numbers,  and  those  reported  throughout  this  paper. 


)  +  -f^ 
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Figure  2. 

(a)  SHG  interferometric  autocorrelation  signal,  (b)  TPAA  signal  for  a  dispersion  pre-compensation  setting  (-579  fs^)  giving 
minimal  pulse  width  for  the  CP-Achromat  lOOx/1,25  oil  microscope  objective.  The  best  fit  was  obtained  assuming  a  sech^ 
input  pulse,  (c)  TPAA  pulse  width  and  amplitude  as  function  of  induced  dispersion  pre-compensation. 


Based  on  the  same  principles  of  dispersive  ray-tracing,  the  induced  dispersion  at  each  setting  of  the  dispersion  pre- 
compensation  unit  (figure  2c)  can  be  calculated.  The  setting  for  the  minimum  TPAA  pulse  width  (figure  2b)  for  the  prism 

compressor  is  calculated  to  correspond  to  GDD  =  -579  fs^,  TOD  =  -209  fs^  and  FOD  =  -156  fs^.  This  includes  the 
dispersion  of  the  coatings  in  the  prism  compressor.  When  the  net  dispersion  of  the  system  is  used  (prism  compressor 

2 

dispersion  plus  objective  dispersion)  it  is  calculated  that  a  15.5  fs  sech  pulse  should  only  broaden  by  -.8  fs,  which  is  in 
excellent  agreement  with  the  measured  pulse  broadening. 

At  these  pulse  durations,  optical  coatings  can  also  produce  significant  effects.  All  turning  optics  in  this  experiment 
were  silvered  mirrors  (Denton  Vacuum,  Inc. ,  coating  FSS-99),  with  the  exception  of  the  dihedral  mirror  R2  used  to  change 
the  beam  height  in  the  prism  compressor.  This  is  a  CVI,  TLMl  coating  for  s-polarization,  45  degree  angle-of  incidence. 
The  effect  of  this  coating  was  investigated  by  taking  multiple  reflections  off  the  mirror,  and  measuring  the  induced  pulse 
broadening  as  a  function  of  number  of  mirror  bounces.  A  19.7  (+  0.1)  fs  pulse  was  found  to  broaden  to  22.3  (±  0.1)  fs  after 
40  reflections  off  this  coating.  This  is  the  expected  pulse  broadening  given  the  calculated  dispersion  values  for  this  coating 
found  in  the  CVI  catalog  [15].  Thus,  for  our  application,  this  is  an  excellent  coating  when  using  sub-20  fs  pulses. 

In  these  experiments,  the  beam  diameter  (FWHM  =  5  mm)  just  fills  the  pupil  of  the  objective  (4.1  mm).  The  drop  in 
intensity  of  the  beam  near  the  pupil  edges,  substantially  reduces  the  contribution  of  any  radially  varying  dispersion  terms, 
such  as  group  delay,  to  the  TPAA  signal.  The  combination  of  a  radially  varying  group  delay,  and  the  fact  that  the 
autocorrelation  is  measured  under  high  NA  conditions  leads  to  a  slight  decrease  in  the  normally  expected  8:1  contrast  of  the 
interferometric  autocorrelation  trace. 


In  conclusion,  we  have  shown  that  with  proper  dispersion  compensation  it  is  possible  to  produce  15-fs  pulses  at  the 
focus  of  a  high  numerical  aperture  objective.  The  entire  system,  prism  dispersion  compensator  plus  objective,  has  been 
modeled  using  a  straight  forward  dispersive  ray-tracing  analysis.  Model  and  experiment  are  in  excellent  agreement.  We  have 
shown  that  by  the  appropriate  choice  of  prism  material  the  pulse  width  can  be  optimized,  below  the  20-30  fs  pulse  width 
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limit  shown  for  other  materials  such  as  SFIO  glass  [12].  Further,  with  proper  accounting  of  all  higher  order  dispersion, 
including  the  quartic  phase,  it  is  reasonable  to  expect  that  systems  can  be  designed  for  even  shorter  pulse  durations.  A 
radially  varying  group  delay  will  also  play  a  more  significant  role  as  the  pulse  duration  is  decreased.  This  must  be  accounted 
for,  by  some  type  of  radially  varying  correction,  or  the  effective  numerical  aperture  over  which  the  pulse  duration  is 
essentially  constant  will  be  reduced. 
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Ultrashort  Pulse  Experiment 
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ABSTRACT 

Until  recently,  experiments  in  ultrashort  pulse  science  have  involved  measuring  the  spectrum  and  autocorrelation  of  the 
input  pulse(s)  and  only  measuring  the  integrated  energy  or  perhaps  time-resolved  energy  of  the  output  signal.  These 
experiments  ignored  the  information  contained  in  the  input  and  output  pulse  phases  and  intensity  profiles.  New  pulse 
measurement  techniques  such  as  frequency-resolved  optical  gating  (FROG),  when  combined  with  older  techniques  such  as 
spectral  interferometry,  now  allow  the  complete  characterization  of  the  pulses.  These  techniques  allow  measurements  of  the 
intensity,  phase,  and  polarization  state  of  ultrashort  pulses  as  functions  of  time  (or  frequency)  and  space.  These  techniques 
work  for  wavelengths  from  the  UV  to  the  IR  and  for  extremely  weak  pulses  and  very  high  power  pulses.  They  also  allow 
entirely  new  classes  of  experiments  for  measuring  ultrafast  phenomena.  Now  the  phases  and  temporal  profiles  of  the  input 
pulses  may  be  measured  and  controlled,  and  the  intensity  and  phase  of  the  output  pulses  can  also  be  measured.  These  new 
measurement  techniques  have  thus  greatly  increased  the  obtainable  information  in  ultrafast  experiments.  This  paper  reviews 
current  pulse  measurement  methods  including  frequency-resolved  optical  gating  and  spectral  interferometry  and  describes  how 
they  are  changing  the  way  that  ultrashort  pulse  experiments  are  performed. 

Keywords:  frequency-resolved  optical  gating,  spectral  interferometry,  autocorrelation,  ultrashort  pulses,  ultrafast  science 

1.  INTRODUCTION 

From  their  first  use,'’^  autocorrelators  were  the  main  device  used  for  measxuing  the  temporal  characteristics  of  short  laser 
pulses.  When  combined  with  a  measured  spectrum,  the  autocorrelation  gave  an  estimate  of  the  pulse  length,  and  of  how  close 
the  width  of  the  pulse  was  to  the  spectral  limit  (the  limit  set  by  the  Fourier  transform  of  the  spectrum  assuming  a  constant 
phase).  As  these  were  the  best  measurements  available,  the  generic  experiment  in  ultrashort  pulse  science  involved  measuring 
the  spectrum  and  autocorrelation  of  the  input  pulse(s)  and  the  integrated  energy  or  perhaps  time-resolved  energy  of  the  output 
signal.  Figure  1  shows  a  schematic  of  this  type  of  experiment.  The  main  difficulty  with  this  generic  experiment  is  that  it 
ignores  most  of  the  information  contained  in  the  input  and  output  pulses.  If  we  ignore  the  spatial  behavior  of  the  pulse  and 
assume  that  it  is  linearly  polarized,  a  pulse  is  defined  by  its  electric  field  as  a  function  of  time,  E(t),  where 

E{t)  =  Rej.^^  expiico^t- (1) 

Here  I(t)  and  ())(t)  are  the  time-dependent  intensity  and  phase  of  the  pulse,  and  cOo  is  a  carrier  frequency.  The  intensity 
contains  information  about  the  relative  amounts  of  energy  contained  in  Ae  pulse  at  given  times,  while  the  time-dependent 
phase  contains  the  frequency  versus  time  information.  The  instantaneous  frequency  of  the  pulse,  ®(t),  is 
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Figure  1  Schematic  of  the  old  generic  experiment  in  ultrafast  science.  The  input  pulse  has 
intensity  Ii(t)  and  phase  (t)i(t),  and  the  output  pulse  has  intensity  Io(t)  and  phase  (l)i)(t).  The 
only  characteristics  of  the  input  pulse  that  are  measured  are  the  autocorrelation  and 
spectrum,  and  only  the  energy  (or  time  resolved  energy)  is  measured  for  the  output  pulse. 
The  input  pulse  and  the  output  pulse  are  not  completely  characterized. 


It' 


(2) 


The  pulse  may  also  be  described  in  the  frequency  domain  as  E(co),  the  Fourier  transform  of  Eft).  Ignoring  the  negative- 
frequency  term  of  the  Fourier  transform  gives 


e{co)  =  Rej.^l(fi)-(0o)  exp(z^(£0  -  Oo))}- 


(3) 


Here  IfocOg)  is  the  spectrum,  and  tpfta-cOc)  is  the  spectral  phase.  While  the  spectrum  is  directly  measurable,  an 
autocorrelation  gives  little  information  about  the  spectral  phase,  so  measuring  the  autocorrelation  and  spectrum  allows  for  a 
large  uncertainty  in  the  behavior  of  the  pulses  that  are  being  measured. 

New  pulse  measurement  techniques,  however,  such  as  frequency-resolved  optical  gating  (FROG)  combined  with  older 
techniques  such  as  spectral  interferometry,  now  allow  the  complete  characterization  of  the  intensity  and  phase  of  ultrashort 
pulses.  These  techniques  work  for  wavelengths  from  the  UV  to  the  IR  and  for  extremely  weak  pulses  and  very  high  power 
pulses.  They  also  allow  entirely  new  classes  of  experiments  for  measuring  ultrafast  phenomena.  Now  the  phases  and  temporal 
profiles  of  the  input  may  be  measured  and  controlled,  and  the  intensity  and  phase  of  the  output  pulses  can  also  be  measured. 
These  new  measurement  techniques  have  thus  greatly  increased  the  information  that  can  be  obtained  in  ultrafast  experiments. 
This  paper  reviews  current  pulse  measurement  methods  including  frequency-resolved  optical  gating  and  spectral 
interferometry  and  describe  how  they  are  changing  the  way  that  ultrashort  pulse  experiments  are  performed. 

2.  AUTOCORRELATIONS  AND  SPECTRA 


2.1.  Autocorrelations 

When  measuring  the  duration  of  an  event,  the  method  normally  used  is  to  compare  it  to  a  faster  event.  We  measure  years 
by  the  passing  of  seasons  or  days.  We  measure  days  or  hours  by  the  passage  of  the  sun,  or  the  swing  of  a  pendulum.  As 
technology  has  advanced,  we  have  compared  fractions  of  seconds  to  vibrations  of  quartz  crystals  or  other  very  rapid  events.  In 
measuring  an  ultrashort  laser  pulse,  however,  the  problem  is  that  the  pulse  is  by  far  the  fastest  available  signal  in  the 
laboratory.  How,  then,  do  we  measure  a  pulse  when  we  don’t  have  anything  faster  than  it  for  a  reference?The  autocorrelation 
was  one  of  the  first  real  attempts  to  answer  this  question.  The  key  to  the  technique  is  to  allow  the  pulse  to  measure  itself. 
Specifically,  autocorrelation  involves  splitting  the  pulse  into  two  replica  pulses,  introducing  a  time  delay  between  the 
replicas  and  spatially  overlapping  them  in  an  instantaneously  responding  nonlinear  medium.  The  autocorrelation  is  then  the 
signal  produced  by  the  nonlinear  interaction  of  the  two  replicas  in  the  medium  as  a  function  of  the  time  delay. 
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Figure  2  SHG  autocorrelations  for  a  spectrally 
limited  10  fs  Gaussian  pulse  (constant  phase), 
and  the  same  pulse  after  adding  17.5  fsVradian  of 
GDD,  250  fsVradian^of  cubic  phase,  and  1400 
fs^adian^  of  quartic  phase.  The  full  width  at  half 
maxima  of  the  autocorrelations  are  14.1  fs,  15.7 
fs,  15.5  fs,  and  15.3  fs,  respectively.  Since  the 
non-spectrally  limited  pulses  have  approximately 
11  fs  full  width  at  half  maxima,  fte  error  in 
estimating  the  pulsewidth  based  on  the 
spectrally  limited  pulseshape  is  small.  Note, 
however,  that  the  different  phase  distortions 
produce  essentially  identical  shapes  of  the  SHG 
autocorrelation. 
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Figure  3  The  SHG  autocorrelations  of  an  initially 
10-fs  spectrally  limited  Gaussian  after  adding  850 
fsVradian^  and  7000  fsVradian^  quartic  phase.  Both 
of  these  pulses  have  14.3  fs  intensity  FWHM,  and 
the  autocorrelations  are  virtually  identical  except 
for  small  differences  in  the  size  of  the  wings  of  the 
autocorrelations.  These  autocorrelations  are  nearly 
identical  for  most  experimental  purposes,  and 
provide  essentially  no  information  about  the  order 
of  the  phase. 


The  most  common  nonlinear  interaction  used  for  autocorrelation  is  second-harmonic  generation  (SHG).  The  signal 
produced  by  an  SHG  autocorrelator  occurs  at  twice  the  center  frequency  of  the  input  pulse  and  has  a  form  given  by 


(4) 


Here  T  is  the  time  delay  between  the  pulses.  Looking  at  Equation  4,  we  see  that  the  signal  is  essentially  obtained  by  using 
one  of  the  replica  pulses — ^here  I(t-'c) —  to  gate  out  part  of  tihe  other  pulse.  Thus  the  autocorrelation  uses  the  pulse  to  measure 
itself  Unfortunately,  this  measurement  yields  a  smeared  out  version  of  I(t)  that  often  masks  complicated  pulse  structure. 
Satellite  pulses,  for  instance,  often  look  exactly  the  same  as  a  long  smooth  tail  on  the  pulse.  Worse,  while  an  estimate  of  the 
pulsewidth  may  be  made  by  assuming  a  specific  pulseshape,  the  estimate  obtained  is  highly  dependent  on  the  assumed 
shape.  This  makes  predictions  based  on  SHG  autocorrelation  unreliable.  There  are  a  number  of  different  autocorrelation 
methods  based  on  other  nonlinearities  or  specific  interferometric.^'*  None  of  these  autocorrelations  eliminate  all  of  these 
problems. 

For  the  moment,  let  us  consider  pulses  with  durations  close  to  the  spectral  limit.  By  measuring  the  spectmm  of  the 
pulse  and  assuming  a  constant  phase,  the  spectrally  limited  pulsewidth  can  be  determined.  Then  direct  calculation  gives  the 
autocorrelation  of  the  spectrally  limited  pulse.  If  the  measured  autocorrelation  of  the  actual  pulse  matches  that  of  the  spectrally 
limited  pulse,  then  the  pulse  must  be  near  the  transform  limit.  Consider,  for  example,  making  an  SHG  autocorrelation  of  a 
spectrally  limited  10  fs  Gaussian.  As  shown  in  Figure  2,  adding  17.5  fsVradian  of  GDD,  250  fsVradian^  of  cubic  phase  or 
1400  fsVradian^  of  quartic  phase  produces  an  only  marginally  different  autocorrelation.  The  autocorrelation  full  width  at  half 
maxima  for  these  four  cases  are  14.1  fs,  15.7  fs,  15.5  fs,  and  15.3  fs,  respectively.  Since  the  non-spectrally  limited  pulses 
have  approximately  1 1  fs  FWHM,  the  error  in  estimating  the  pulsewidth  based  on  the  spectrally  limited  pulse  shape  is 
small.  Unfortunately  differentiating  between  thesedifferentphases  based  on  these  autocorrelations  is  impossible.  This  is  not 
true  just  for  pulses  that  are  close  to  the  spectral  limit.  Figure  3  shows  the  autocorrelations  of  an  initially  10-fs  spectrally 
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limited  Gaussian  after  adding  850  fsVradian^  of  cubic  phase  or  7000  fs^'/radian^  of  quartic  phase.  Both  of  these  pulses  have 
14.3  fs  intensity  FWHM.  Again  the  autocorrelations  are  very  similar  except  for  small  changes  in  the  wings  of  fte 
autocorrelations.  For  experimental  purposes,  these  autocorrelations  are  nearly  identical,  and  provide  almost  no  information 
about  the  phase  distortion  that  is  present.  Even  the  parity  of  the  distortion  cannot  be  determined.  Other  forms  of 
autocorrelation  are  all  relatively  insensitive  to  the  exact  phase  as  well;  although,  some  third-order  autocorrelators  are  sensitive 
to  the  parity  of  the  dominant  phase  term. 

2.2.  The  Spectrum 

As  mentioned  above,  the  spectrum  is  directly  measurable,  and  entire  books  have  been  written  on  how  to  build  optical 
spectrometers.®  Thus  the  measurement  is  generally  a  matter  of  buying  the  appropriate  spectrometer.  The  spectmm  is  related 
to  the  field  by 


/»  = 


J  j  exp(-f£0t)(fi 


(5) 


The  definition  given  in  Equation  3  shows  that,  when  combined  with  the  appropriate  spectral  phase,  the  spectmm  gives  the 
Fourier  transform  of  the  electric  field  of  the  pulse.  Without  the  phase,  however,  nothing  is  known  about  the  temporal 
behavior  of  the  pulse.  Even  though  there  have  been  attempts  to  deduce  the  phase  from  the  spectmm  and  autocorrelation, 
taking  the  two  measurements  at  the  same  time  still  does  not  remove  all  the  uncertainties  in  the  behavior  of  the  pulse.'®’'' 

3 .  FREQUENCY-RESOLVED  OPTICAL  GATING 


3.1.  Introduction 


The  real  goal  for  developing  ultrashort  pulse  measurements  has  been  full  characterization  of  the  intensity  and  phase  of  the 
pulse,  and  recently  a  burst  of  new  techniques  operating  in  the  time-frequency  domain  have  achieved  full  characterization.'^.  Of 
these  techniques,  the  most  widely  applied  is  frequency-resolvedoptical  gating  (FROG)'®.  FROG  involves  measuring  a 
spectrogram  of  the  pulse.  Mathematically,  the  spectrogram  of  a  pulse  is  given  by 


exp^-icotldn 


(6) 


Here  g(t-Tj  is  a  variable-delay  gate  function.  The  spectrogram  then  is  the  set  of  all  parts  of  £ftj  that  are  gated  out  as  t  is 
varied.  To  within  an  absolute  phase  factor,  the  spectrogram  of  the  pulse  determines  £(tj  completely. 

What,  however,  is  the  gate  function  in  the  laboratory?  The  gate  pulse  must  be  shorter  than  the  pulse  to  be  measured  or 
spectral  information,  and  hence  phase  information,  will  be  lost.  In  the  field  of  acoustics  where  the  spectrogram  has  been  used 
for  decades.  The  sound  pulses  being  measured  are  relatively  long,  and  it  is  a  trivial  matter  to  produce  electronic  pulses  can  be 
used  as  a  gate.  For  ultrafast  pulses,  however,  the  pulse  to  be  measured  is  usually  by  far  the  fastest  event  in  the  laboratory. 
Thus,  as  in  an  autocorrelator,  the  pulse  must  be  used  to  gate  itself. 

3.2.  Combining  the  Autocorrelator  and  Spectrometer 

Since  the  pulse  must  be  used  to  gate  itself,  a  nonlinear-optical  process  must  be  used  as  the  gating  interaction.  Consider  a 
SHG  autocorrelator.  The  signal  field  produced  by  the  autocorrelator  is 

=  (7) 

where  as  before  t  is  time  and  ris  the  time  delay  between  the  two  pulses  as  they  arrive  at  the  second-harmonic  material.  If 
that  signal  field  is  then  spectrally  resolved  in  a  spectrometer,  the  result  is 
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Figure  4  Schematic  of  an  SHG  FROG  device.  The  SHG  FROG  is  simply  a  second-harmonic  generation 
autocorrelator  followed  by  a  spectrometer  that  spectrally  resolves  the  autocorrelation  signal  as  a  function  of 
the  time  delay,  T,  between  the  input  pulse 
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(8) 


This  measurement,  also  called  the  SHG  FROG  trace,*"*’**  is  a  spectrogram  with  the  gate  pulse  equal  to  the  signal  pulse. 
Experimentally,  the  old  methods  of  characterization  involved  taking  separate  autocorrelations  and  spectra.  FROG  is  not 
particularly  more  complicated  because  it  simply  involves  moving  the  spectrometer  after  the  autocorrelator  and  taking  spectra 
of  autocorrelations.  Figure  4  shows  a  schematic  of  an  SHG  FROG.  Note  too  that  other  FROG  beam  geometries  are  possible 
based  on  polarization  gating  (PG)’*,  self  diffraction  (SD),  third-harmonic  generation  (THG)*’,  and  many  other  near 
instantaneous  nonlinear  interactions.*^  The  variety  of  possible  beam  geometries  makes  FROG  useful  fi’om  the  ultraviolet  to 
the  infrared  and  for  multi-shot  or  single-shot  experiments.*^  While  the  forms  of  the  signal  field  Esig(t,T)  may  be  differentffom 
that  given  above  for  SHG,  the  FROG  measurements  are  all  given  by 
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3.3.  Phase  Retrieval 

Mathematically,  using  the  pulse  to  gate  itself  complicates  the  problem.  The  FROG  trace  must  be  inverted  to  find  the 
intensity  and  phase  of  the  pulse.  Spectrogram  inversion  algorithms  used  to  obtain  the  field  (or  intensity  and  phase)  require 
knowledge  of  the  gate  function.  Since  the  gate  is  related  to  the  pulse  and  is  unknown,  such  algorithms  cannot  be  used  for 
FROG.  Other  methods  are  needed.  The  solution  to  the  inversion  of  FROG  data  is  to  recast  the  problem  as  a  two-dimensional 
phase  retrieval  problem.*"*"**-**  Instead  of  thinking  of  the  problem  as  a  gating  problem,  consider  the  signal  field  Esig(t,x)  to  be 
the  Fourier  transform  with  respect  to  T  of  a  new  quantity  Esig(t,Q).  We  note  that  once  Esig(t,Q)  is  found,  the  pulse  field  is 
easily  obtained  by  setting  Q  equal  to  zero.  That  is  E(t)  =  Esig(t,£2=0)  to  within  a  complex  multiplicative  constant  that  is 
typically  unimportant.  Thus  finding  Esig(t,Q)  is  equivalent  to  finding  E(t).  Rewriting  the  expression  for  the  FROG  trace  in 
terms  of  E,ig(t,Q)  gives 


(^co,t^=  J  J  £„.j(/,£2)exp(-z£at-j£2r)drdO 


(10) 


Thus  the  measured  quantity  lfrog(co,t)  is  the  squared  magnitude  of  the  two-dimensional  Fourier  transform  of  Esig(t,  £2.).  The 
FROG  trace  thus  yields  the  magnitude,  but  not  the  phase,  of  the  two-dimensional  Fourier  transform  of  Esig(t,£2).  Finding  the 
phase  ofthe  Fourier  transform  of  £'j4(if,f2^  is  the  two-dimensional  phase  retrieval  problem,  which  is  a  solved  problem  from 
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Figure  5  Generic  experiment  based  on  frequency-resolved  optical  gating  (FROG).  Instead  of  measuring  only 
the  spectrum  and  autocorrelation  of  the  input  pulse  and  the  energy  (or  time-resolved  energy)  of  the  output 
pulse  as  in  Figure  1,  the  experiment  measures  the  full  intensity  and  phase  of  both  the  input  and  output 
pulses.  Thus  essentially  all  information  contained  in  the  temporal  or  spectral  behavior  of  the  pulse  is 
obtained. 

image  science.'®  The  measured  trace  and  knowledge  of  the  mathematical  form  of  Esig(t,Q.)  are  sufficientto  retrieve  the  full 
intensity  and  phase  of  a  pulse  by  using  an  iterative  Former  transform  algorithm'^-^®. 

3.4.  Generic  Experiments  using  FROG 

Rather  than  considering  the  FROG  phase  retrieval  algorithm,  which  has  been  discussed  in  previous  publications, '^’2° 
consider  the  impact  of  knowing  how  to  measure  the  full  intensity  and  phase  of  a  pulse  on  potential  experiments  in  ultrafast 
science.  Instead  of  relying  on  partial  information  such  as  that  provided  by  an  experiment  of  the  form  shown  in  Figure  1,  now 
the  full  intensity  and  phase  may  be  measured  in  a  new  generic  experiment  shown  in  Figure  5. 

Using  such  an  experiment,  it  is  now  possible  to  measure  the  full  intensity  and  phase  (or  equivalently  the  electric  field)  of 
both  the  input  and  output  pulses.  Thus,  for  instance  if  the  experiment  were  to  simply  replace  the  block  labeled  “Experiment” 
in  Figure  5  with  a  length  of  some  material,  the  complex  transmission  coefficientof  the  material  could  be  obtained  over  the 
wavelength  of  the  pulse  by  simply  dividing  the  measured  electric  field  of  the  input  pulse  by  the  measured  electric  field  of  the 
output  pulse.  Thus  essentially  all  information  contained  in  the  temporal  or  spectral  behavior  of  the  pulse  is  obtained. 
Moreover,  by  using  pulse  shaping  techniques^',  it  is  also  now  possible  to  control  the  pulse  shape  of  the  input  pulse.  These 
pulse  shaping  techniques  have  already  been  applied  to  a  variety  of  fields  including  nuclear  magnetic  resonance  in  dipolar 
broadened  spin  systems^^  and  quantum  control  of  wavepackets^^,  and  ^eatly  increase  the  number  of  possible  experiments. 
Measurement  techniques  such  as  FROG  make  using  pulse  shaping  techniques  much  more  feasible. 

4.  SPECTRAL  INTERFEROMETRY 


4.1.  Need  to  Measure  Weak  Pulses 

Despite  the  great  breakthrough  in  pulse  measurement  that  FROG  represents,  one  characteristic  of  FROG  can  make  it  less 
than  ideal  for  some  experiments.  Specifically,  FROG  requires  a  nonlinear  interaction,  which  limits  its  use  to  pulses  strong 
enough  to  produce  a  measurable  nonlinear  signal.  For  SHG  FROG,  this  limits  the  measurement  to  optical  pulses  with 
energies  of  about  a  picojoule.  As  shown  by  Walmsley  and  Wong^®,  all  techniques  for  measuring  an  isolated  pulse  must  use 
a  nonlinear  interaction,  not  just  FROG.  Thus  all  techniques  for  measuring  isolated  pulses  have  similar  energy  requirements. 
Many  ultrafast  nonlinear-optical  material-characterization  experiments,  however,  yield  signal  pulse  energies  of  femtojoules  or 
less.  In  general,  only  the  signal-pulse  energy  has  been  measured  in  such  experiments.  Many  researchers,  however,  have 
realized  that  significantly  more  material  information  is  available  if  the  signal  pulse  is  better  characterized.  Chemla  and 
coworkers  and  Patkar  and  coworkers^®  have  use  several  pulse-characterizationmethods  to  characterizethe  signal  pulses  in 
four-wave-mixing  experiments  performed  on  semi-conductors.  Consequently,  they  have  significantly  improved  on  our 
understanding  of  the  physics  of  multiple-quantum  wells,  for  example.  Measuring  an  isolated  ultraweak  pulse  by  itself  is  a 
daunting  task,  but  by  combining  spectral  interferometry  and  FROG  it  is  now  possible  to  completely  characterize  the  intensity 
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and  phase  of  pulses  with  energies  in  the  zeptojoule  range  (1  zeptojoule  =  10'^'  joules).  This  combination  makes  possible 
experiments  that  rely  on  understanding  the  intensity  and  phase  behavior  of  ultraweak  pulses. 

4.2.  Spectral  Interferometry 

Spectral  interferometry  (SI),  a  technique  for  measuring  the  phase  differencebetween  two  pulses — say  a  referenceand  an 
unknown  pulse — ^was  first  introduced  by  Froehly  and  coworkers,^’^’^*  and  it  has  been  used  for  several  applications  since 
then.^^"^*  It  involves  simply  combining  at  a  beamsplitter  the  two  pulses  whose  relative  phase  is  desired  and  directing  the 
two  pulses  collinearly  into  a  spectrometer.  The  spectrum  of  the  two  pulses  is  simply  given  by: 

Is, (CO)  =  I,,f((o)  +  /„„* (CO)  +  2.J/„/(m)  ^|I,„t((») cos[<^„„t (m) -  t^((o) -cot]  (11) 

where  Ire/co)  and  Iu„k(co)  are  the  spectra  of  the  reference  and  unknown  (weak)  pulses,  respectively,  and  T  is  the  delay  between 
the  two  pulses.  The  delay  (only  one  delay  is  necessary)  is  chosen  to  yield  fringes  in  the  sum  spectmm.  The  two  individual 
pulse  spectra  are  measured  by  blocking  one  beam  and  measuring  the  spectrum  of  the  other  and  vice  versa.  The  only 
remaining  unknown  in  Equation  11,  the  phase  difference,  <punk((o)-(pref(ffl)-t£tt,  is  then  determined  by  the  SI  spectrum,  Isi(ffl). 
It  is  easily  extracted  non-iteratively  using  one  of  several  well-known  fringe-inversion  techniques. In  fact,  the  relative 
delay  also  emerges  from  the  analysis  (it  is  the  linear  term)  and  hence  does  not  need  to  be  independently  measured.  An 
advantage  of  SI  is  that  it  is  a  type  of  heterodyne  technique,  which  can  act  to  amplify  the  weak  pulse.  For  example,  the 
reference  pulse  may  be  chosen  to  be,  say,  M  times  more  intense  than  the  unknown  pulse,  yielding  fringes  that  are  4M  times 
as  intense  as  the  spectrum  of  the  unknown  pulse.  The  only  requirement  of  SI  is  that  the  spectrum  of  the  unknown  pulse  lies 
within  that  of  the  reference  pulse. 

Note  that  SI  is  a  linear  technique  (the  reference  pulse  eliminates  the  need  for  a  nonlinear  interaction),  so  it  is  highly 
sensitive.  Also,  if  the  phase  of  the  reference  pulse  is  known,  SI  may  be  used  to  obtain  both  the  spectrum  and  spectral  phase  of 
the  pulse.  Using  Equation  3,  this  then  gives  the  electric  field  of  the  pulse  and  completely  characterizes  the  pulse.  During  the 
majority  of  the  time  that  SI  has  existed,  however,  measuring  the  reference  pulse  was  impossible.  Thus  Ae  technique  was 
limited  to  measuring  only  phase  differences  or  pulses  wiA  phase  distortions  Aat  were  large  compared  to  the  phase  distortions 
of  the  reference,  and  making  Ae  assumption  that  the  reference  pulse  was  spectrally  limited. 

4.3.  Combining  FROG  and  Spectral  Interferometry 

With  the  invention  of  techniques  such  a  FROG,  it  is  now  a  simple  matter  to  measure  ultraweak  pulses  because  in  most 
cases  Ae  ultraweeik  pulse  is  produced  from  a  strong  pulse  directly  from  a  laser.  Thus  in  almost  all  cases  a  reference  pulse  is 
available.  We  first  use  FROG  to  characterize  Ae  reference  pulse  directly  from  Ae  laser  (although  any  method  Aat  yields  Ae 
reference-pulse  phase  vs.  frequency,  <Pref(<B),  will  suffice).  By  characterizing  the  pulse  using  FROG,  we  avoid  Ae  necessity  of 
assuming  a  spectrally  limited  pulse.  We  Aen  use  SI  to  measure  Ae  frequency-domain  phase  difference  between  the  two 


Figure  6  Generic  experiment  based  on  TADPOLE  (Ae  combination  of  FROG  and  spectral  interferometry). 
The  experiment  measures  Ae  All  intensity  and  phase  of  boA  Ae  input  and  output  pulses,  even  for 
extremely  weak  pulses.  TADPOLE  measurements  have  been  performed  on  pulses  wiA  energies  as  low  as 
42  zeptojoules  (42x10'^’  joules).  Moreover,  by  using  Ae  non-dispersed  dimension  in  the  spectrometer,  it  is 
possible  to  measure  Ae  intensity  and  phase  as  a  Anction  of  position  along  the  entrance  slit  of  Ae 
spectrometer  as  well  as  time. 
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pulses.  The  FROG  and  SI  measurements  together  yield  the  full  intensity  and  phase  of  the  unknown  ultraweak  ultrashort 
pulse  in  the  frequency  domain.  Thus,  the  combination  of  FROG  and  SI,  which  is  sometimes  called  temporal  analysis  by 
dispersing  a  pair  of  light  e-fields  (TADPOLE),  provides  a  nearly  general  technique  for  measuring  even  the  weakest  ultrashort 
laser  pulses.  So  far,  this  technique  has  been  used  to  measure  pulses  with  energies  as  low  as  42  zeptojoules,  or  42  x  10 
joules.^^  Moreover,  since  only  one  physical  dimension  of  the  detector  is  used  for  the  spectral  measurement,  spatial 
information  may  be  obtained  by  using  the  other  dimension  in  an  imaging  spectrograph.  Thus  TADPOLE  allows 
measurements  of  the  intensity  and  phase  as  a  function  of  position  as  well  as  time. 

4.4.  Polarization  Measurement 

In  addition  to  measuring  intensity  and  phase  as  a  function  of  time  (or  frequency)  and  space,  TADPOLE  may  also  be  used 
to  measure  the  time-dependent  polarization  state  of  a  pulse  by  using  a  dual-beam  geometry  to  measure  the  time  dependent 
phase  of  two  orthogonal  polarizations.^^-^'*  Consider  a  pulse  propagating  in  the  z  direction.  If  we  define  two  components  x 
andy  of  the  electric  field  that  are  orthogonal  to  the  z  and  each  other,  we  may  perform  a  spectral  interferometry  measurement  on 
each  component  using  a  linearly  polarized  referencepulse  at  45°  with  respect  to  the  x-axis,  so  that  it  has  equal  x  and  y 
components.  The  experimental  pulse  x  andy  components  may  have  arbitrarily  varying  intensities  and  phases  as  a  function  of 
time.  One  performs  a  spectral  interferometric  measurement  on  the  x  andy  components  of  the  pulse  by  introducing  a  fixed 
delay  between  the  referencepulse  and  the  signal,  combining  the  referenceand  signal  collinearly,  separating  the  combined 
reference  and  signal  into  x  andy  components,  and  measuring  the  x  andy  components  separately  in  a  spectrometer. 

The  measured  signals  have  the  same  form  as  for  single  chaimel  spectral  interferometry 

4  (m)  =  4(<d)  +  4t(ft))  +  2  .^4/6))  V4/c(<«)  cos[4,(<»)  -  4/ffl)  -  tor  J  (12) 

Here  i  indicates  the  x  ory  component.  Since  the  same  pulse  is  used  as  a  reference  for  the  two  polarizations,  the  relative  phase 
of  the  two  components  is  obtained  as  well  as  the  full  intensity  and  phase  (electric  field)  of  the  two  components  individually. 
Thus  having  measured  two  orthogonal  electric  fields  of  the  pulse  and  the  phase  relationship  between  those  fields,  the  full  time 
dependent  field  may  be  obtained  including  the  polarization  state.  This  method  has  already  been  used  to  study  coherent 
processes  in  quantum  confined  structures  by  measuring  the  intensity,  phase  and  polarization  state  of  coherent  four-wave¬ 
mixing  signals  from  GaAs  quantum  wells.^”* 

5.  CONCLUSION:  THE  NEW  GENERIC  EXPERIMENT 

Until  recently,  experiments  in  ultrashort  pulse  science  have  involved  measuring  the  spectrum  and  autocorrelation  of  the 
input  pulse(s)  and  only  measuring  the  integrated  energy  or  perhaps  time-resolved  energy  of  the  output  signal.  These 
experiments  ignored  the  information  contained  in  the  input  and  output  pulse  phases  and  intensity  profiles.  New  pulse 
measurement  techniques,  however,  such  as  frequency-resolved  optical  gating  (FROG)  combined  with  older  techniques  such  as 
spectral  interferometry,  now  allow  the  complete  characterization  of  the  intensity,  phase  and  polarization  state  of  ultrashort 
pulses  as  functions  of  time  and  position.  These  techniques  work  for  wavelengths  from  the  UV  to  the  IR  and  for  extremely 
weak  pulses  and  very  high  power  pulses.  They  also  allow  entirely  new  classes  of  experiments  for  measuring  ultrafast 
phenomena.  Now  the  phases  and  temporal  profiles  of  the  input  may  be  measured  and  controlled,  and  the  intensity  and  phase 
of  the  output  pulses  can  also  be  measured.  These  new  measurement  techniques  have  thus  greatly  increased  the  information 
that  can  be  obtained  in  ultrafast  experiments.  This  talk  will  review  current  pulse  measurement  methods  including  frequency- 
resolved  optical  gating  and  spectral  interferometry  and  describe  how  they  are  changing  the  way  that  ultrashort  pulse 
experiments  are  performed. 
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ABSTRACT 

Nonlinear  absorption  through  laser-induced  breakdown  (LIB)  offers  the  possibility  of  localized  energy  deposition  in 
linearly  transparent  media  and  thus  of  non-invasive  surgery  inside  the  eye.  The  general  sequence  of  events  -  plasma 
formation,  stress  wave  emission,  cavitation  -  is  always  the  same,  but  the  detailed  characteristics  of  these  processes  depend 
strongly  on  the  laser  pulse  duration.  The  various  aspects  of  LIB  are  reviewed  for  pulse  durations  between  80  ns  and  100  fs, 
and  it  is  discussed,  how  their  dependence  on  pulse  duration  can  be  used  to  control  the  efficacy  of  surgical  procedures  and 
the  amount  of  collateral  effects. 

Key  Words:  Optical  breakdown,  intraocular  microsurgery,  ultrashort  laser  pulses,  breakdown  thresholds,  -  plasma 
formation,  moving  breakdown,  plasma  transmission,  shock  waves,  cavitation,  self  focusing. 

1.  INTRODUCTION 

Nonlinear  absorption  through  laser-induced  breakdown  (LIB)*  ^  can  occur  at  material  surfaces  as  well  as  inside  media 
which  are  transparent  at  low  light  intensities.  Nevertheless,  the  characteristics  of  LIB  are  very  different  in  both  cases.  At 
material  surfaces  LIB  competes  with  material  ablation  based  on  linear  absorption  and  has,  besides  some  characteristic 
features,  many  aspects  in  common  with  the  latter^’"*.  LIB  inside  of  linearly  transparent  media  offers  a  possibility  of 
localized  energy  deposition  which  can  be  achieved  by  no  other  optical  means.  This  unique  feature  enables  non-invasive 
surgery  inside  the  eye’’^,  and  it  has  been  suggested  to  apply  it  for  the  design  of  3-dimensional  storage  elements®.  The 
results  described  in  this  review  were  obtained  with  distilled  water  used  as  a  model  substance  for  the  transparent  media  of 
the  eye.  This  simplification  guaranteed  reproducibility  of  the  exjrerimental  results  and  is  justified  by  the  fact  that  the 
thresholds  for  optical  breakdown  in  distilled  water  are  very  similar  to  the  breakdown  thresholds  in  the  ocular  media  ’ .  The 
review  focuses  attention  on  the  dependence  of  the  breakdow'n  events  on  laser  pulse  duration  and  discusses,  how  that 
dependence  can  be  used  to  control  the  efficacy  of  surgical  procedures  and  the  amount  of  collateral  effects. 


2.  OPTICAL  BREAKDOWN 

Optical  breakdown  consists  of  the  generation  of  large  amounts  of  free  electrons  (S  10*®  cm"^  ^■’°)  in  the  medium  by 
multiphoton  ionization  (MPI)  and  cascade  ionization  via  inverse  bremsstrahlung  absorption  (Cl).  The  rate  of  MPI  has  a 
very  strong  intensity  dependence  (oc  7^  where  K  is  the  number  of  photons  required  for  ionization),  whereas  the  intensity 
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dependence  of  the  rate  of  Cl  is  much  weaker  (« 1  when  electron  losses  from  the  focal  volume  are  neglected)''.  When  the 
laser  pulse  duration  is  reduced,  a  higher  intensity  of  the  laser  radiation  is  required  for  the  breakdown  to  be  completed 
within  the  pulse  duration.  Since  the  MPI  rate  increases  more  strongly  with  I  than  the  Cl  rate,  MPI  gains  ever  more 
importance  with  shorter  pulse  durations.  For  pulses  below  »  1  ps.  Cl  is  furthermore  limited  by  time  constraints,  since  one 
doubling  sequence  in  the  ionization  cascade  lasts  at  least  6  fs  (for  water  and  'k  =  1064  mn)'^,  whereas  MPI  can  occur 
"instantaneously"  The  changing  interplay  between  MPI  and  Cl  with  decreasing  pulse  duration  is  one  main  reason  for  the 
changes  in  the  breakdown  characteristics  as  a  whole. 

Another  major  factor  is  the  changing  energy  distribution  between  electrons  and  heavy  particles  which  influences 
the  energy  density  in  the  plasma  and  thus  the  radiant  energy  threshold  for  breakdown.  Optical  energy  is  deposited  into  the 
medium  through  generation  of  free  electrons.  With  ns-pulses,  a  temperature  equilibrium  between  electrons  and  heavy 
particles  is  achieved  during  the  pulse  through  recombination  processes,  and  therefore  the  average  energy  density  is  high  . 
With  ultrashort  pulses,  however,  very  little  energy  has  been  transferred  to  the  heavy  particles  at  the  end  of  the  laser  pulse. 
An  equilibrimn  temperature  develops  only  after  the  laser  pulse.  The  equilibrium  temperature  will  thus  be  much  lower  than 
in  the  case  of  the  ns-pulses,  particularly  because  the  specific  heat  of  the  electrons  is  much  smaller  than  that  of  the  ions  and 
other  heavy  particles'^.  Due  to  the  smaller  average  energy  density  in  plasmas  produced  with  ultrashort  laser  pulses,  the 
radiant  energy  threshold  for  breakdown  is  smaller  than  with  ns-pulses,  and  the  mechanical  effects  after  breakdown  (shock 
wave  emission  and  cavitation)  are  less  extensive. 

A  third  reason  for  changes  of  the  breakdown  characteristics  with  decreasing  pulse  duration  is  the  increasing  role  of 
self  focusing  effects  which  goes  along  with  the  higher  laser  powers  required  to  complete  breakdown  in  a  shorter  period  of 
time.  Self  focusing  changes  the  breakdown  threshold,  the  plasma  transmission,  and,  probably  most  importantly,  the  shape 
of  the  breakdown  region  as  well  as  the  energy  distribution  within  that  region  (Fig.  I). 

The  statements  made  above  about  the  three  main  factors  influencing  the  pulse  duration  dependence  of  breakdown 
are  in  the  following  discussed  in  more  detail  and  supported  by  experimental  evidence. 

3.  BREAKDOWN  THRESHOLDS 

The  evolution  of  the  free  electron  density  under  the  influence  of  a  laser  pulse  can  be  described  by  the  rate  equation^’ 

Yt  =  +  '/case  P-SP-V^P^  (1) 

The  first  two  terms  describe  the  creation  of  free  electrons  through  multiphoton  and  cascade  ionization,  whereas  the  losses 
by  electron  diffusion  out  of  the  focal  volume  and  by  recombination  are  represented  by  the  last  two  terms.  In  order  to  predict 
LBB  thresholds,  several  authors  have  solved  this  equation  neglecting  different  terms  and  using  different  expressions  for  the 
individual  terms®  ' Recently  Noack  et  al.'®  presented  an  analysis  based  on  all  four  terms  including  the  recombination 
term,  which  was  based  on  experimental  data'^.  For  the  other  terms,  the  expressions  of  Kennedy’s  first  order  model'^  were 
used.  When  solving  Eq.  (1),  it  has  to  be  considered  that  Cl  can  only  take  place  when  free  electrons  are  already  present  in 
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Figure  1:  Optical  breakdown  at  various  pulse  durations,  (a)  Pulse  duration  =  76  ns,  wavelength  X  =  750  nm,  focusing 
angle  19°,  £  =  50  mJ,  EIE,h  =  9,  picture  taken  120  ns  after  breakdown  (the  line  is  a  slit  used  for  simultaneous  streak 
photography,  and  has  otherwise  no  physical  meaning);  (b)  =  6  ns,  A  =  1064  nm,  9=  22°,  E  =  8.2  mJ,  E/Eth  =  60,At  = 

10  ns;  (c)  fjr  =  30  ps,  X  =  1064  nm,  9  =  14°,  E  =  740  pj,  EIE,h  =  150,  ^  r  =  8  ns;  (d)  =  100  fs,  2  =  580  nm,  6>  =  16°, 

£  =  35  pj,  EIE,h  =  200,  zl  r  =  3  ns.  The  light  is  always  incident  from  the  right,  the  bars  represent  a  length  of  100  pm. 


the  interaction  volume.  In  pure  media  like  distilled  water,  MPI  is  required  to  provide  the  initial  electrons.  Therefore,  Cl 
was  only  included  in  the  calculations  after  the  probability  of  finding  an  electron  in  the  interaction  volume  V  had  risen 
above  50%  through  MPI.  The  time  constraint  for  Cl  given  by  the  time  At  for  an  electron  to  gain  enough  energy  through 
inverse  bremsstrahlung  absorption  to  ionize  another  electron  was  considered  by  evaluatimg  the  contribution  of  Cl  at  a 
retarded  time  (r- At).  Following  Keimedy’^’^’,  we  chose  A/  =  30  fs.  The  breakdown  threshold  was  defined  as  the 

minimum  irradiance  where  a  critical  electron  density  of  p„  =  10^®  cm“^  was  exceeded  during  the  laser  pulse.  This 
threshold  was  obtained  by  solving  Eq.  (1)  iteratively  for  different  intensities. 

Table  1  summarizes  various  experimental  threshold  values  Iso  (50%  breakdown  probability)  obtained  at  different 
laser  wavelengths"’’®’’^  together  with  the  calculated  threshold  values  /rate  for  the  same  laser  parameters'®.  For  all  pulse 
durations  except  60  ps  the  calculated  and  experimental  thresholds  deviate  by  less  than  a  factor  of  1.4. 
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pulse  duration 


wavelength 

[nm] 


measured  spot 
diameter  [pm] 


ho 

[10”  Wcm'^] 


76ns 

750 

6  ns 

532 

532 

3ps 

580 

580 

580 

E50 

[|iJl 


Eso 

[Jcm'^] 

/rate 

[10^^  Wcm'^] 

1750 

0.27 

174 

0.35 

16.8 

0.62 

2.6 

7.10 

1.4 

41.0 

1.1 

84.0 

20 

5.3 
5.6 
5.0 
5.0 

4.4 


5500 

39 

4.1 

0.51 

0.29 

0.17 


0.23 

0.29 

2.8 

8.5 

47.6 

111.0 


Table  1:  Measured  breakdown  thresholds  Eso,  Iso  and  F5o,and  calculated  breakdown  thresholds  hate  for  various  laser  pulse 
durations. 


The  numerical  solution  of  Eq.  (1)  does  not  only  allow  the  calculation  of  breakdown  thresholds,  but  also  provides 
insight  into  the  differences  between  the  breakdown  mechanisms  at  different  pulse  durations.  To  discuss  these  differences, 
Figure  2  presents  the  evolution  of  the  free  electron  density  for  laser  pulses  of  different  duration,  and  Figure  3  shows  the 
predicted  threshold  as  a  function  of  laser  pulse  duration  for  two  different  wavelengths. 


Figure  2:  Evolution  of  the  free  electron  density  at 
threshold  irradiance  for  laser  pulses  of  different 
durations  .  All  curves  were  calculated  for  h  =  580  nm 
and  6  pm  focus  diameter.  The  irradiance  maximum  of 
the  Gaussian  laser  pulse  is  reached  at  t  =  0. 


"Cl 

Figure  3:  Irradiance  threshold  /rate  as  a  function  of  laser 
pulse  duration  calculated  for  2  =  1064  nm  and  1  =  580 
nm  (6  pm  focus  diameter). 


In  pure  media  as  distilled  water,  the  initial  electrons  for  the  breakdown  process  must  be  generated  by  MPl,  even  at 
pulse  durations  where  breakdown  is  otheivvise  dominated  by  Cl.  For  nanosecond  pulses,  the  electron  density  initially  rises 
only  slowly  through  multiphoton  absorption,  because  jj^^p  is  small  at  the  threshold  irradiance  for  breakdown.  After  the 

creation  of  the  first  electron  in  the  interaction  volume  at  /  «  0.8  r^,  the  beginning  Cl  raises  the  electron  density  by  almost 


1 1  orders  of  magnitude  within  a  small  fraction  of  the  laser  pulse  (»  1  ns)  (Fig.  2;  76  ns,  6  ns).  The  high  rate  of  Cl  is  due  to 
the  relatively  high  irradiance  required  for  MPI  of  the  initial  electron.  The  growth  of  the  free  electron  density  is  slowed 
down  at  high  electron  concentrations  by  the  rapidly  increasing  influence  of  recombination  (oc/?^).  Because  of  the  strong 
influence  of  recombination,  is  reached  near  the  peak  of  the  laser  pulse.  As  long  as  this  is  the  case,  the  threshold 
intensity  required  to  overcome  the  recombination  losses  is  nearly  independent  of  the  pulse  duration  (Fig.  3). 

For  shorter  laser  pulses,  the  ionization  cascade  becomes  slower  with  respect  to  the  pulse  duration.  Around  1  ns,  the 
pulse  duration  becomes  comparable  to  the  rise  time  of  the  electron  density  through  Cl,  and  for  shorter  pulses  a  higher 
^casc  is  required  to  complete  the  cascade  during  the  pulse.  The  rising  threshold  intensity  implies  that  Cl  can 
compensate  the  recombination  losses  even  during  the  second  half  of  the  laser  pulse  when  the  intensity  has  dropped  below 
its  maximum  value.  For  60  ps  pulses,  is  therefore  reached  after  the  peak  of  the  laser  pulse  (Fig.  2;  60  ps). 

Whereas  for  nanosecond  pulses  MPI  was  only  required  to  produce  seed  electrons  for  the  cascade,  it  gains  ever  more 
importance  with  decreasing  pulse  duration  (i.e.  increasing  threshold  intensity)  because  of  its  strong  intensity  dependence 
(oc/^).  For  pulse  durations  around  and  below  1  ps,  MPI  provides  a  large  number  of  seed  electrons  for  the  ionization 
cascade  and  contributes  significantly  to  the  generation  of  fi'ee  electrons,  particularly  during  the  first  half  of  the  laser  pulse 
(Fig.  2:  3  ps,  100  fs).  For  <  Ips,  the  critical  electron  density  could  not  be  reached  by  Cl  alone,  regardless  of  irradiance, 
because  one  doubling  sequence  in  the  ionization  cascade  takes  >  6  fs.  Nevertheless,  Cl  is  still  the  dominant  mechanism  at 
all  pulse  durations  investigated  in  our  study,  only  at  pulse  durations  below  »  40  fs,  the  majority  of  the  free  electrons  is 
produced  by  MPl’^’^°. 

The  above  discussion  shows,  how  the  interplay  of  MPI,  Cl  and  recombination  changes  with  laser  pulse  duration. 
The  influence  of  recombination  has  been  underestimated  in  previous  investigations”’’^.  Its  neglegt  does  not  lead  to  very 
different  threshold  values’®,  but  it  obscures  the  role  played  by  recombination  in  the  energy  transfer  between  electrons  and 
heavy  particles.  The  recombination  losses  occurring  during  the  ns  pulses  lead  to  a  continuous  heating  of  the  heavy 
particles  and  thus  to  an  increase  of  the  average  energy  density  in  the  interaction  volume.  Therefore,  a  much  larger  radiant 
exposure  F50  is  required  for  breakdown  at  long  pulse  durations,  particularly  at  76  ns,  than  with  ps  and  fs  pulses  (Table  1). 
For  pulse  durations  <  3  ps,  F50  varies  only  slightly,  because  recombination  losses  during  the  laser  pulse  are  always  minor 
and  thus  almost  independent  of  pulse  duration.  Other  factors  which  could  also  contribute  to  the  dependence  of  F50  on  pulse 
duration  but  were  not  considered  in  our  calculations,  are  alterations  of  the  plasma  size  at  threshold  and  variations  of  p^^ . 
They  will  be  discussed  in  the  following  sections. 

4.  PLASMA  FORMATION  AT  SUPERTHRESHOLD  ENERGIES 

When  pulses  with  superthreshold  peak  irradiance  are  applied,  the  plasma  grows  during  the  laser  pulse.  For  ps-  and  ns- 
pulses,  the  plasma  is  first  formed  at  the  beam  waist,  and  then  grows  "upstream"  toward  the  incoming  laser  beam”’^° 
("moving  breakdown").  Hardly  any  plasma  is  formed  beyond  the  laser  focus,  because  most  of  the  laser  light  is  already 
absorbed  upstream”’^’’^^  ("plasma  shielding").  For  picosecond  pulses,  the  position  of  the  plasma  front  at  each  time  is 


38 


approximately  defined  by  iso-intensity  lines  with  1  =  Itk,  aird  the  plasma  length  at  the  end  of  the  pulse  is  proportional  to 
where  p  =  EIE,h  =  IHth  denotes  the  normalized  pulse  energy  or  peak  irradiance,  respectively".  For  nanosecond 
pulses,  the  irradiance  threshold  is  lowered  to  a  value  7,^  after  the  plasma  formation  has  started,  because  UV-plasma 
radiation  generates  free  electrons  in  the  vicinity  of  the  plasma  which  make  the  MPI-generation  of  seed  electrons 
superfluous.  The  movement  of  the  breakdown  wave  is  approximately  determined  by  iso-intensity  lines  with  /  =  7,^,  and  ns- 
plasmas  are  therefore  considerably  longer  than  ps-plasmas  at  the  same  value  of  ".  The  relative  deviation  between  the 
length  of  ns-  and  ps-plasmas  at  equal  p  is  largest  for  large  focusing  angles  and  near  the  breakdown  threshold  (up  to  six¬ 
fold  at  <9  =  22°).  With  ps-  pulses,  the  UV  plasma  radiation  does  not  modify  the  breakdown  threshold,  because  there  are 
always  enough  seed  electrons  produced  by  MPI. 

With  fs-pulses,  the  spatial  extension  of  the  laser  pulse  is  shorter  than  the  plasma  length  at  superthreshold  energies. 
Therefore,  the  breakdown  starts  in  front  of  the  laser  focus,  and  the  breakdown  wave  moves  with  the  laser  pulse,  leaving 
free  electrons  in  its  wake^^.  Despite  of  the  different  direction  of  the  breakdown  wave,  the  length  of  the  breakdown  region 
scales  in  a  similar  way  with  p  as  in  the  case  of  ps-pulses  . 

The  plasma  growth  in  upstream  direction  occurring  during  ns-  and  ps-  laser  pulse  reduces  the  light  absorption  at 
each  point  of  the  plasma  and  thus  limits  the  average  energy  density  within  the  plasma  to  values  which  do  not  differ  much 
from  the  energy  density  at  threshold.  With  fs-pulses,  the  energy  density  in  the  plasma  is  limited  by  the  fact  that  for 
increasing  yffthe  onset  of  MPI  starts  further  and  further  upstream.  This  reduces  the  peak  power  and  irradiance  nearer  to  the 
laser  focus,  like  in  the  case  of  the  longer  pulses. 


5.  ENERGY  DEPOSITION 

Only  the  energy  absorbed  in  the  sample  is  effective  for  plasma-mediated  surgery  (and  for  material  processing  in  general); 
light  transmission  through  the  plasma  as  well  as  scattering  and  reflection  by  the  plasma  reduce  the  efficacy  of  the  process. 
The  plasma  absorption  cannot  be  measured  directly,  but  can  be  deduced  by  measuring  the  transmission,  scattering,  and 
reflection:  A  =  \-  T-S-R.  Nahen  and  Vogef  ’  found  in  a  recent  study  performed  with  pulse  durations  of  6  ns  and  30  ps 
that  the  back  reflection  from  the  plasmas  into  the  focusing  lens  amounts  to  less  than  1.7%  for  both  pulse  durations,  and 
that  the  forward  scattering  into  an  angle  of  ±30°  around  the  optical  axis  amounts  to  less  than  0.5%  for  ns-pulses  and  to 
less  than  7.6%  for  ps-pulses.  Near  threshold,  considerably  more  light  was  transmitted  through  the  plasma  than  reflected  or 
scattered.  Well  above  threshold  (fix.  100)  the  relative  importance  of  transmission  decreased,  but  transmission  was  still 
twice  as  large  as  scattering  and  reflection  together.  For  all  laser  parameters,  considerably  less  light  was  reflected  or 
scattered  by  the  plasma  than  absorbed.  The  plasma  absorption  is  thus  approximately  given  by  A  w  (1-7)  and  its  assessment 
for  a  larger  parameter  range  is  possible  through  transmission  measurements  alone.  This  differs  from  plasma  formation  at 
surfaces,  where  reflection  plays  a  large  role^.  The  difference  is  due  to  the  fact  that  inside  a  homogeneous,  linearly 
transparent  medium  the  breakdown  front  moves  during  the  laser  pulse  (see  section  4).  The  moving  breakdown  limits  the 
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electron  density  at  each  location  within  the  plasma.  The  plasma  frequency  remains  therefore  smaller  than  the  frequency  of 
the  laser  light,  and  the  laser-plasma-coupling  is  not  impaired. 

Figure  4  shows  measured  transmission  values'*  as  a  function  of  laser  pulse  duration  for  two  different  values  of  the 
dimensionless  laser  pulse  energy  p.  The  transmission  is  small  in  the  nanosecond  range,  but  considerably  larger  for  ps- 
pulses.  When  the  pulse  duration  is  further  reduced  into  the  femtosecond  range,  the  transmission  decreases  again.  At  first 
sight,  the  experimental  data  are  not  easily  understood  when  the  maximum  electron  density  in  the  plasma  is  assumed  to  be 
independent  of  the  laser  pulse  duration.  This  assumption  {/?„  =  const.)  was  made  for  the  threshold  calculations  in 
section  3  and  led  to  a  very  good  agreement  between  predicted  and  measured  threshold  values. 


To  explain  the  T(r^)  dependence,  one  has  to  consider  that  the  measured  transmission  values  represent  a  time 
average  over  the  whole  laser  pulse  duration.  They  must,  hence,  be  compared  with  the  complete  time  evolution  of  the 
electron  density  during  the  laser  pulse,  and  not  just  with  the  maximum  density.  The  time  averaged  absorption  coefficient 
aLm  I'^^r  threshold  can  be  estimated  from  the  absorption  cross  section  for  inverse  bremsstrahlung  cr^  and  the  temporal 
evolution  of  the  free  electron  density  obtained  by  integration  of  Eq.  (1)  to  be'® 
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Figure  5  shows  the  calculated  values  for  as  a  function  of  pulse  duration'®.  The  absorption  coefficient  is  largest 
for  nanosecond  pulses  (800  cm  ')  because  of  the  high  electron  density  throughout  the  laser  pulse  (see  Fig.  2).  Since  the 
evolution  of  the  free  electron  density  is  very  similar  for  pulse  durations  between  1  and  100  ns,  the  absorption  coefficient  is 
almost  constant  in  this  regime.  With  decreasing  laser  pulse  duration,  high  electron  densities  are  only  encountered  during 
the  second  half  of  the  laser  pulse  (Fig.  2:  60  ps,  3  ps).  Therefore  the  time  averaged  absorption  drops  and  reaches  a 
minimum  value  of  150  cm  '  for  a  pulse  duration  around  2-3  ps.  For  even  shorter  pulses,  the  increasing  role  of  multiphoton 


Pulse  duration 

Figure  4:  Plasma  transmission  as  a  function  of  laser 
pulse  duration  for  /?  =  6  and  p  =  60.  The  measured  data 
points'*  are  connected  with  lines  to  facilitate 
orientation. 


Figure  5:  Calculated  time  averaged  absorption 
coefficient  auB  as  a  function  of  laser  pulse  duration  Tl 
(2.  =  580  nm,  =  6  pm). 
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absorption  leads  to  an  increasing  electron  density  in  the  first  half  of  the  pulse  (Fig.  2;  100  fs),  so  that  again  a  larger 
fraction  of  the  laser  pulse  energy  is  absorbed. 

The  qualitative  behavior  of  the  calculated  absorption  coefficient  near  threshold  (Fig.  5)  agrees  well  with  the  trends 
observed  in  the  measured  transmission  data  (Fig.  4):  The  transmitted  energy  shows  only  a  weak  dependence  on  pulse 
duration  in  the  nanosecond  range.  With  decreasing  pulse  duration  experiments  and  calculations  show  an  increase  in 
transmission  (a  decrease  in  absorption).  The  measured  transmission  coincides  with  the  predicted  minimum  of  cklib  at  2-3 
ps,  and  the  decreasing  transmission  for  femtosecond  pulses  can  be  explained  by  the  enhanced  generation  of  free  electrons 
during  the  first  half  of  the  laser  pulse. 

The  absolute  amount  of  the  calculated  absorption  coefficients  was  found  to  be  in  good  agreement  with  experimental 
data  obtained  from  the  measured  transmission  data  T  and  plasma  length  /  assuming  the  validity  of  Beers  law 
(a  =  -InT/  /  ).  It  is  interesting  to  note  that  the  average  absorption  coefficient  of  laser-produced  plasmas  in  water  is  for  all 
pulse  durations  considerably  smaller  than  the  linear  absorption  coefficient  of  water  at  >1  =  3  pm  (w  17700  cm" )  ,  or  the 
linear  absorption  coefficient  of  corneal  tissue  at  193  nm  (a  40000  cni"')^®.  This  means  that  energy  deposition  based  on 
linear  absorption  can,  under  certain  conditions,  be  more  localized  than  plasma-mediated  energy  deposition. 

The  pulse  duration  dependence  of  the  absorption  coefficient  explains  the  trends  in  the  T{  Ti  )  dependence,  but  alone 
it  does  not  suffice  to  explain  the  large  difference  between  the  transmission  of  ns-  and  ps-plasmas  (at  ^  =  6  ns-plasmas 
transmit  only  »1/10  of  the  light  which  passes  through  ps-plasmas,  see  Fig.  4).  The  second  factor  contributing  to  the  low 
transmission  of  ns-plasmas  is  that  they  are  longer  than  ps-  and  fs-  plasmas  at  equal  p  and  focusing  angle  (section  4).  The 
difference  in  length  is  particularly  pronounced  near  threshold.  This  is  the  reason  why  the  transmission  immediately  drops 
to  50%  once  a  ns-  plasma  is  formed,  but  it  only  drops  to  90-95%  for  ps-  and  fs-plasmas  at  threshold  energies^'. 

6.  MECHANICAL  EFFECTS 

6.1  Shock  wave  emission  and  cavitation 

Table  2  gives  an  overview  of  how  the  mechanical  effects  created  by  LIB  depend  on  the  laser  pulse  duration.  It  summarizes 
the  values  of  the  shock  wave  pressure  at  the  plasma  rim  and  at  12  mm  distance  from  the  plasma,  and  presents  the 
conversion  efficiency  of  absorbed  light  energy  into  cavitation  bubble  energy.  The  shock  wave  emission  in  the  vicinity  of 
the  plasma  was  investigated  by  streak  photography  with  high  spatial  and  temporal  resolution  (3  gm  and  <200  ps, 
respectively)'*’^^.  The  pressure  in  the  far-field  was  determined  by  hydrophone  measurements'*,  and  the  cavitation  bubble 
energy  was  obtained  through  a  determination  of  the  bubble  size  by  acoustic  means  ’  . 

The  data  of  Table  2  show  that  the  mechanical  effects  accompanying  breakdown  decrease  with  decreasing  pulse 
diuation.  This  is  true  not  only  in  absolute  terms  (due  to  the  lower  radiant  energy  threshold  for  breakdown),  but  also  in 
relative  terms:  a  smaller  percentage  of  the  incident  light  energy  is  converted  into  mechanical  energy.  This  becomes  most 
obvious  by  looking  at  the  reduction  of  the  conversion  efficiency  of  light  energy  into  cavitation  bubble  energy.  It  is. 
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pulse  duration 

P=EIE,h 

pressure  at 
plasma  rim 
[GPa] 

pressure  at 

12  mm  distance 
[MPa] 

conversion  efficiency  of 
absorbed  light  energy  into 
cavitation  bubble  energy  [%] 

76  ns 

6 

10 

4.0 

22.0 

6  ns 

60 

30 

2.5 

22.5 

60  ps 

60 

10 

0.65 

13.5 

3  ps 

60 

1.7 

0.23 

11.0 

100  fs 

60 

1.0 

0.06 

3.0 

Table  2;  Dependence  of  shock  wave  pressure  and  cavitation  on  laser  pulse  duration'*.  The  data  for  76  ns  pulses  are  given 
at  y0=  6  instead  of  >3=  60  because  of  the  limited  laser  pulse  energy  available. 

furthermore,  reflected  in  the  decrease  of  the  pressure  at  the  plasma  rim  for  shorter  pulse  durations.  The  smaller  peak 
pressure  indicates  a  smaller  energy  density  deposited  in  the  breakdown  volume.  The  pulse  duration  dependence  of  the  far- 
field  pressure  supports  the  above  observations.  It  should  be  noted  however,  that  the  far-field  pressure  depends  both  on  the 
pressure  at  the  plasma  rim  and  on  the  plasma  size.  The  similarity  law  states  that  the  shock  pressure  at  a  distance  r  from  a 
source  with  radius  depends  only  on  the  ratio  r//-o,  if  the  pressure  inside  the  source  is  constant^®.  Since,  at  equal  p,  the 
plasma  radius  ro  decreases  with  decreasing  laser  pulse  duration,  the  relative  distance  /V/'o  increases  for  r  =  constant.  Thus  a 
decreasing  shock  pressure  at  a  fixed  distance  r  would  be  expected  for  shorter  laser  pulses,  even  if  the  energy  density  in  the 
plasma  was  independent  of  the  pulse  duration. 

6.2  Energy  density  in  the  breakdown  region 

We  estimated  the  average  energy  density  in  the  breakdown  region  by  comparing  the  absorbed  energy  with  the  breakdown 
volume  deduced  from  photographs.  With  6-ns-pulses  at  1064  nm,  energy  densities  of  30-40  kJcm  *  were  reached  which  are 
far  beyond  the  evaporation  enthalpy  of  water  and  more  than  6  times  higher  than  the  energy  density  in  TNT^°.  The  high 
energy  density  in  ns-plasmas  explains  the  high  conversion  efficiency  into  cavitation  bubble  energy  as  well  as  the  high 
shock  wave  pressures  observed  (Table  2).  For  76-ns  pulses,  the  maximum  pressure  in  the  plasma  is  limited  by  the  fact  that 
the  plasma  and  bubble  expansion  start  already  during  the  laser  pulse  (see  Fig  la;  the  bubble  is  already  quite  large  near  the 
beam  waist  whereas  its  formation  has  just  started  further  upstream),  therefore  the  peak  pressure  is  lower  than  after  6-ns 
pulses.  With  fs-pulses,  the  energy  deposition  is  isochoric,  but  the  average  energy  density  in  the  breakdown  region  is  only 
about  1  kJ  cm  *  ’*  ,  i.e.  lower  than  the  evaporation  enthalpy  at  constant  pressure.  The  stress  wave  generation  features 
therefore  a  stronger  contribution  of  thermoelastic  effects,  and  the  peak  pressure  of  the  stress  transients  is  much  lower  than 
with  ns-pulses.  Correspondingly,  the  conversion  efficiency  into  bubble  energy  also  drops  to  only  3%. 

The  high  energy  density  in  the  ns-  plasmas  is  associated  with  an  easily  visible  plasma  radiation,  whereas  the  lower 
energy  density  in  the  fs-plasmas  goes  along  with  a  lack  of  a  visible  plasma  radiation:  the  occurrence  of  breakdown 
manifests  itself  merely  by  bubble  formation  at  the  laser  focus.  It  should  be  noted,  however,  that  the  plasma  radiation 
reflects  mainly  the  electron  temperatures  and  that  there  is  little  correlation  to  the  average  energy  density  in  the  sample 
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when  the  temperatures  of  the  electrons  and  heavy  particles  are  different,  i.  e.  with  ps-  and  fs-pulses.  The  different  visibility 
of  the  plasma  radiation  is  probably  largely  due  to  the  change  of  its  duration  and  size  with  different  laser  pulse  durations, 
and  not  so  much  a  consequence  of  a  different  intensity  or  of  changes  in  its  spectral  characteristics. 

The  reduction  of  energy  density  in  the  plasma  with  decreasing  pulse  duration  leads  to  a  change  of  energy  partition: 
whereas  a  smaller  percentage  of  the  incident  light  energy  is  converted  into  mechanical  energy,  a  larger  percentage  is  used 
for  evaporation  of  the  liquid  within  the  breakdown  volume^*.  We  found  that  the  ratio  (Emeoh/Ev)  at  60  is  about  12.5:1 
for  6-ns  pulses  and  0.25:1  for  100-fs  pulses'®'^*  Besides  the  reduction  of  the  radiant  energy  threshold  with  decreasing 
pulse  duration,  this  change  of  energy  partition  contributes  largely  to  a  less  disruptive  character  of  fs-breakdown. 

6.3.  Possible  explanations  for  the  pulse  duration  dependence  of  energy  density 

While  most  physical  mechanisms  involved  in  optical  breakdown  are  linked  to  the  irradiance  of  the  laser  light,  the  radiant 
energy  required  for  breakdown  and  the  resulting  energy  density  in  the  interaction  volume  are  of  more  practical  interest. 
The  change  of  the  radiant  energy  threshold  with  decreasing  pulse  duration  and  the  change  of  the  average  energy  density  in 
the  plasma  are  closely  linked,  but  not  exactly  the  same:  they  are  equivalent  only,  if  the  size  of  the  interaction  volume  and 
the  plasma  absorption  are  equal.  In  the  following,  we  shall  discuss  possible  factors  influencing  the  pulse  duration 
dependence  of  the  energy  density. 

Maximum  electron  density.  The  energy  density  in  the  plasma  volume  will  decrease  with  pulse  duration,  if  the 
maximum  concentration  of  free  electrons  reached  during  the  laser  pulse  is  smaller  for  shorter  pulse  durations.  This 
possibility  has  been  suggested  by  Kennedy  et  al.^*,  who  assume  that  p„  is  10^°  cm^  for  ns-  and  60-ps  plasmas,  but  lo'* 
cm‘^  for  3-ps  and  fs-plasmas.  It  has,  however,  not  yet  been  further  investigated  to  date;  measurement  data  of  the  free 
electron  concentration  after  optical  breakdown  in  water  are  only  available  for  30-ns  plasmas,  where  Barnes  and 
Rieckhotf  ®  found  p  »  10^®  cm'^  at  pulse  energies  well  above  threshold.  Our  threshold  calculations  presented  in  section  3 
yielded  good  agreement  with  experimental  values  under  the  assumption  of  a  constant  value  =  10  cm  for  all  pulse 
durations.  It  is  unlikely  that  p„  largely  exceeds  10^'  cm'^  at  any  pulse  duration  because  this  is,  for  X  =  1064  nm,  the 
critical  plasma  density  above  which  a  major  part  of  the  incident  laser  light  will  be  reflected.  At  ultrashort  pulse  durations, 
where  the  whole  energy  at  the  end  of  the  laser  pulse  is  stored  in  the  form  of  free  electrons  and  hardly  any  energy  has  been 
transferred  to  the  heavy  particles,  it  is  also  unlikely  that  p„  is  much  less  than  10^'  cm'^.  The  energy  of  6.5  eV  required  to 

produce  a  quasi-free  electron  in  water^^  corresponds  to  E  =  1.04  lO  '*  J.  An  electron  density  of  p=  2.5x10^' cm  ^ 
therefore  corresponds  to  an  energy  density  of  2.5  kJ  cm  ^  (if  the  free  electrons  do  not  have  very  much  kinetic  energy), 
equal  to  the  evaporation  enthalpy  of  water  at  room  temperature.  This  means  that  p„  must  be  2.5x10^'  cm  ^  to  produce  a 
pure  vapor  bubble.  If  thermoelastic  effects  take  part  in  the  generation  of  the  bubbles  produced  by  ultrashort  laser  pulses, 
p^  may  be  smaller,  but  certainly  not  by  several  orders  of  magnitude^^.  For  ns-pulses,  the  evaporation  enthalpy  can  be 

exceeded  with  maximum  electron  concentrations  below  10^’  cm”^  because  of  the  energy  transfer  from  electrons  to  heavy 
particles  during  the  pulse  through  recombination. 


43 


Energy  distribution  between  electrons  and  heavy  particles.  The  energj'  density  in  the  plasma  changes,  even  if  the 
maximum  concentration  of  free  electrons  is  the  same  for  all  pulse  durations:  with  long  pulses,  a  temperature  equilibrium 
between  electrons  and  heavy  particles  is  achieved  during  the  pulse  through  recombination  processes,  and  therefore  the 
energy  density  is  high.  With  ultrashort  pulses,  however,  very  little  energy  has  at  the  end  of  the  laser  pulse  been  transferred 
to  the  heavy  particles.  An  equilibrium  temperature  develops  only  after  the  laser  pulse,  and  the  equilibrium  temperature  will 
thus  be  much  lower  than  in  the  case  of  the  ns-pulses'^. 

Average  temperature  of  electrons  and  heavy  particles.  With  ns-  and  ps-pulses  at  superthreshold  energies,  the 
plasma  grows  towards  the  incoming  laser  beam  during  the  first  part  of  the  laser  pulse  and  is  heated  during  the  second  half 
of  the  pulse,  particularly  at  the  upstream  side  of  the  plasma’ This  heating  results  in  an  increase  of  the  average 
temperature  of  the  elctrons  and  heavy  particles.  It  is  evidently  most  pronounced  with  long  pulse  durations  where  the 
amount  of  surplus  energy  deposited  in  the  plasma  is  largest.  At  76  ns,  however,  the  increase  of  the  energy  density  is 
counteracted  by  the  expansion  of  the  plasma  during  the  laser  pulse. 

Plasma  size.  With  ns-pulses,  the  breakdown  threshold  decreases  during  the  laser  pulse  due  to  the  UV-radiation 
emitted  by  recombining  free  electrons  and  bremsstrahlung  (see  section  4).  This  leads  to  an  increase  of  the  plasma  size  and 
thus  to  a  reduction  of  the  average  energy  density  within  the  plasma. 

Sharpness  of  plasma  boundaries.  After  fs-breakdown,  heating  of  the  liquid  is  observed  in  front  of  the  actual 
breakdown  zone  where  a  bubble  is  produced,  particularly  for  large  j^values  (Figure  6).  Similar  observations  have  not  been 
made  for  ns-  or  ps-pulses.  This  means  that  in  fs-breakdown  a  larger  part  of  the  absorbed  energy  cannot  be  converted  into 
bubble  energy  than  with  the  longer  pulse  durations.  The  reason  for  this  phenomenon  is  probably  that  a  large  number  of 
free  electrons  is  produced  by  MPI  even  in  regions  where  is  not  completely  reached  (in  Figure  2  the  corresponding 
feature  is  that  a  large  number  of  free  electrons  is  already  present  at  the  beginning  of  tlie  100  fs  laser  pulse). 

We  can  conclude  that  there  are  several  factors  influencing  the  pulse  duration  dependence  of  energy  density  within 
the  breakdown  volume  which  are  partly  opposing  each  other.  The  most  important  factor  setting  the  general  trend  seems  to 
be  the  energy  distribution  between  electrons  and  heavy  particles  at  the  end  of  the  laser  pulse. 


7.  SELF  FOCUSING 


In  the  previous  sections  we  have,  for  the  sake  of  simplicity,  neglected  that  all  aspects  of  the  breakdown  process  with 
ultrashort  laser  pulses  are  actually  modified  by  self  focusing.  For  self  focusing  leading  to  beam  collapse  to  occur  with  a 
Gaussian  beam,  a  critical  power^® 


P„  =3.77- 


cA- 


(3) 
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has  to  be  surpassed,  where  iij  is  the  value  of  the  nonlinear  refractive  index.  Since  the  irradiance  and  power  thresholds  for 
optical  breakdown  increase  when  the  pulse  duration  is  reduced,  the  role  of  self-focusing  becomes  ever  more  important. 


Figure  6:  Optical  breakdown 
region  after  a  100  fs  pulse 
with  35  pJ  pulse  energy 
(y9=200),  photographed  at 
different  times  after 
breakdown.  The  photographs 
taken  1  and  2  ps  after  break¬ 
down  are  slightly  defocused 
to  visualize  the  zones  of 
elevated  temperature  outside 
the  cavitation  bubble.  The 
laser  light  was  incident  from 
the  right,  the  bar  represents  a 
length  of  100  pm. 


Self-focusing  will 

-  diminish  the  energy  threshold  Eso  for  breakdown,  because  it  reduces  the  cross  section  of  the  laser  beam, 

-  produce  an  elongated  breakdown  region,  because  it  causes  the  location  of  the  focus  to  move  during  the  laser  pulse  leaving 
a  trace  which  looks  like  a  filament  '  , 

-  change  the  energy  distribution  in  the  breakdown  region,  because  the  region  splits  into  the  part  of  the  "filament"  and  a 

"shoulder"  In  the  shoulder,  the  energy  density  in  the  plasma  is  reduced,  or  the  irradiance  may  even  stay  below  the 

breakdown  threshold,  and 

-  increase  the  light  transmission  through  the  breakdown  region^',  because  the  absorption  coefficient  is  smaller  in  the 
shoulder  region,  which  contains  a  fairly  large  part  of  the  laser  light  energy^®'^*. 

Theoretical  predictions  of  the  self-focusing  thresholds  for  ultrashort  laser  pulses  face  some  difficulties;  The  n2- 
values  are  probably  smaller  than  those  measured  for  longer  pulses,  because  the  laser  pulses  are  now  shorter  than  the  time 
constant  of  the  molecular  Kerr  effect  which  is  the  predominant  mechanism  responsible  for  the  refractive  index  change 
with  ns-and  ps  pulses^®  '”.  With  fs-pulses,  the  electronic  Kerr  effect  should  gain  relative  importance.  This  change  and  the 
n2-values  for  fs-pulses  have,  however,  not  yet  been  investigated  in  detail.  Furthermore,  self  focusing  is  modified  by  group 
velocity  dispersion  which  can  be  neglected  at  longer  pulse  durations’^’^'. 

Experimental  evidence  indicates  that  with  6-ns-pulses,  self  focusing  plays  a  role  at  threshold  energies  only  for  very 
small  focusing  angles  <  2°  and  at  infrared  wavelengths  where,  due  to  the  lower  rate  of  MPl,  the  required  breakdown  power 
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is  higher  than  at  visible  wavelengths’*.  With  30-ps  pulses  of  1064  nm  wavelength,  self  focusing  is  also  observed  at  larger 
focusing  angles  as  used  in  ocular  surgeiy,  but  only  at  energies  far  above  the  breakdown  threshold”.  With  fs-pulses,  self 
focusing  plays  a  role  for  all  focusing  angles  even  at  tlireshold  and  in  the  visible  wavelength  range’^.  We  could  identify  the 
occurrence  of  self-focusing  effects  with  fs  pulses  through 

-  the  change  of  the  shape  of  the  breakdown  region  observed  in  Figs.  6  and  Id  in  comparison  to  Figs,  la  -  c,  and  by 

-  the  observation  of  continuum  generation  at  all  energies  E  >  Eso-  Continuum  generation  is  like  self-focusing  a  self-phase 
modulation  effect^^  and  in  water  it  is  known  to  occur  always  at  powers  above  the  threshold  for  self-focusing'*^. 

8.  CLINICAL  CONSIDERATIONS 

The  radiant  exposure  threshold  Fsa  for  fs-breakdown  at  tissue  surfaces  and  the  ablation  rate  at  threshold^  '*  are  similar  to 
the  characteristic  parameters  of  ablation  using  excimer  laser  pulses  at  193  nm  wavelength"*^,  and  not  very  different  from 
those  of  Q-switched  IR-laser  ablation  with  wavelengths  near  3  pm'*'*.  This  means  that  plasma-mediated  surface  ablation  - 
even  with  ultrashort  laser  pulses  -  has  no  fundamental  advantage  as  compared  to  ablation  based  on  linear  absorption. 
Apparently,  it  has  also  little  advantage  as  compared  to  plasma-mediated  ablation  with  ps-pulses;  Stern  et  al.'*^  who 
investigated  plasma-mediated  corneal  ablation  with  pulse  durations  between  8  ns  and  65  fs  found  that  the  amount  of  tissue 
damage  produced  by  30-ps  pulses  at  50  pJ  was  so  slight  that  no  additional  improvements  in  excision  ultrastructure  resulted 
when  shorter  pulses  at  even  lower  pulse  energies  were  used  (at  a  tissue  surface,  the  plasma  is  not  confined,  and  its 
expansion  produces  therefore  little  disruptive  effects).  The  situation  is  completely  different,  when  the  application  site  is 
located  within  a  tissue  as,  for  example,  the  eye.  A  noninvasive,  localized  energy  deposition  inside  the  eye  can  only  be 
achieved  through  nonlinear  absorption  of  light  which  at  low  intensities  penetrates  the  tissue.  The  reduced  ablation 
threshold  of  ultrashort  laser  pulses  can  also  be  of  advantage,  because  the  plasma  confinement  leads  to  much  stronger 
mechanical  effects  than  at  surfaces. 

Evaporation  of  the  tissue  within  the  plasma  volume  is  the  primary  tissue  effect  at  all  pulse  durations'*^,  but  the 
degree  of  collateral  mechanical  effects  varies  with  pulse  duration.  Nanosecond  breakdown  features  a  very  high  efficacy  of 
energy  deposition  (Fig.  4),  and  a  very  high  conversion  rate  of  light  energy  into  mechanical  energy  (up  to  90%  It  is, 
hence,  well  suited  where  disruptive  effects  are  desired.  Ultrashort  laser  pulses  have  a  lower  energy  threshold  for 
breakdown!  than  ns-pulses  (F50  drops  by  a  factor  of  »  160  when  the  pulse  duration  is  reduced  from  6  ns-to  100  fs,  see 
Table  1),  a  lower  efficacy  of  energy  deposition  (Fig.  4),  and  a  lower  conversion  rate  into  mechanical  effects  (Table  2,  and 
Ref  28).  As  a  consequence,  et^aporation  is  more  important,  and  the  tissue  effects  are  much  less  disruptive.  Ultrashort 
pulses  allow  therefore  the  generation  of  much  finer  effects  than  can  be  achieved  with  ns-pulses. 

Ultrashort  laser  pulses  are  advantageous  where  precise  cuts  with  little  collateral  damage  are  desired,  as  in  laser- 
keratomileusis,  or  in  the  preparation  of  a  flap  for  LASIK'*'*'**.  For  many  tasks  in  intraocular  microsurgery,  however,  the 
laser  effects  need  not  be  so  fine  as  possible  with  femtosecond  laser  pulses.  In  tlrese  cases,  the  use  of  ns-pulses  may  even  be 
advantageous  as,  for  example,  in  posterior  capsulotomy.  Capsulotomy  with  ps-  or  fs-  pulses  requires  the  application  of 
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pulse  series  with  a  large  number  of  pulses  to  produce  several  long  cuts  of  different  directions  in  the  lens  capsule.  This 
technique  requires  a  fairly  large  amount  of  laser  energy  and  is  relatively  time  consuming.  The  cuts  imply,  furthermore,  a 
larger  risk  of  producing  damage  at  the  curved  surface  of  the  intraocular  lens  than  the  application  ns  pulses,  because  in  the 
latter  case  a  few  individually  aimed  pulses  suffice  to  achieve  the  surgical  aim.  The  ns-pulses  are,  hence  advantageous  for 
capsulotomies  just  because  of  their  disruptive  character"*’. 

The  intraocular  application  of  fs-pulses  is  complicated  by  self-focusing,  particularly  when  the  application  site  is 
located  deeply  inside  the  eye.  In  this  case,  the  localizability  the  laser  effects  will  suffer  from  an  elongation  of  the 
breakdown  region  (Liu  et  al.^  report,  for  example,  that  threshold  effects  in  silica  with  a  diameter  of  only  1  pm  had  a  length 
of  10-15  pm,  and  Figs.  Id  and  6  show  a  strongly  elongated  breakdown  region  at  energies  well  above  threshold).  Laser 
applications  in  the  cornea  as  in  laser  keratomileusis  will,  in  contrast,  hardly  be  influenced  by  self-focusing,  because  the 
propagation  distance  through  the  nonlinear  medium  is  very  short. 
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ABSTRACT 

Extensive  research  of  ultrashort  ocular  damage  mechanisms  has  shown  that  less  energy  is  required  for  retinal 
damage  for  pulses  shorter  than  one  nanosecond.  Laser  minimum  visible  lesion  (MVL)  thresholds  for  retinal 
damage  from  ultrashort  (i.e.  <  1  ns)  laser  pulses  occur  at  lower  energies  than  in  the  nanosecond  (ns)  to 
microsecond  (ps)  laser  pulse  regime.  We  review  the  progress  made  in  determining  the  trends  in  retinal  damage 
from  laser  pulses  of  one  nanosecond  to  one  himdred  femtoseconds  in  the  visible  and  near-infrared  wavelength 
regimes.  We  discuss  the  most  likely  damage  mechanism(s)  operative  in  this  pulse  width  regime  and  discuss 
implications  on  laser  safety  standards. 

Keywords:  eye,  safety,  laser,  nonlinear  optics,  retinal  damage 


1.  BmiODUCTION 

Several  enabling  laser  technologies  are  allowing  the  proliferation  of  ultrashort  laser  systems  in  the  visible  (VIS) 
and  near  infrared  (NIR)  spectral  regions  (i.e.  the  retinal  hazard  regime).  These  wavelengths  easily  propagate 
through  the  eye  and  can  result  in  retinal  damage  and  functional  vision  loss  if  proper  laser  protective  measures  are 
not  employed.  Unfortunately,  there  exists  no  current  maximum  permissible  exposure  limit  based  on  biological 
data  for  pulses  shorter  than  1  ns  in  diuration.  In  recent  works  by  Cain,  et  al.*’^’^  MVL  threshold  data  for 
ultrashort  single  laser  pulses  show  a  decrease  in  the  amount  of  pulse  energy  needed  to  cause  observable  retinal 
damage  below  1  ns.  Also,  laser  induced  breakdown  (LIB),  self-focusing  (SF),  supercontinuum  generation  (SC) 
and  stimulated  Brillouin  scattering  (SBS)  have  been  studied''  to  determine  their  influence  on  retinal  damage 
thresholds. 

In  this  paper  we  review  the  measured  thresholds  for  minimum  visible  lesion  (MVL)  and  for  nonlinear  optical 
phenomena  considered  in  previous  studies  and  the  recent  studies  of  melanin  mediated  damage.  We  will  consider 
visible  and  near  infrared  MVL  studies  and  their  impact  on  determining  new  damage  mechanisms.  A  review  of 
all  the  MVL  literature  from  the  past  couple  of  decades  indicates  that  the  single-pulse  MVL  data  exhibits  three 
broad  trends  (see  Figure  1).  For  exposures  longer  than  20  ps  there  is  a  region  where  near  constant  irradiance 
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(W/cm^)  is  required  for  retinal  damage.  For  pulses  from  1  ns  to  20  ps,  diffusion  of  heat  is  negligible  during  the 
exposure  and  the  fluence  (J/cm^)  required  for  retinal  damage  is  nearly  constant.  The  third  regime  occurs  for 
pulse  durations  shorter  than  one  nanosecond  where  the  data  shows  that  it  takes  less  energy  to  create  retinal 
damage  than  for  longer  pulses.  We  will  discuss  the  region  of  retinal  damage  for  pulse  durations  below  one 
nanosecond. 
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Figure  1.  MVL  data  collected  during  the  previous  thirty  years.  The  data  shows  three  trends  as  discussed  in  the  text. 


2.  ULTRASHORT  LASER  PULSE  RETINAL  DAMAGE 

We  have  published  several  studies  on  measurements  of  the  MVL  threshold  for  both  the  visible^  and  near- 
infrared^  spectral  regimes.  Figure  1  shows  our  measured  thresholds,  which  are  included  with  data  from  many 
other  groups  over  the  past  thirty  years.  Our  experiments  provide  the  majority  of  data  for  pulsewidths  below  one 
nanosecond,  and  concur  with  most  of  the  sub-nanosecond  data  points  taken  by  other  researchers. 

When  considering  ultrashort  retinal  damage  several  phenomena  should  be  considered  that  might  influence  the 
threshold  or  type  of  damage  seen  for  these  short  pulses  (see  Figure  2).  Because  ultrashort  laser  pulses  can 
possess  significant  peak  powers  and  irradiance,  nonlinear  propagation  effects  need  to  be  considered'^.  For  all 
pulse  durations  above  one  nanosecond,  the  propagation  of  laser  light  to  the  retina  can  be  considered  linear  and 
the  spot  size  can  be  calculated  from  the  wavelength,  aberrations  and  laser  divergence.  For  ultrashort  laser  pulses, 
the  retinal  fluence  and  wavelength  can  be  different  than  that  resulting  from  linear  propagation.  We  have  reported 
the  thresholds  measured  for  several  nonlinear  optical  phenomena^,  which  are  shown  in  Figure  3. 
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Figure  2.  This  figure  highlights  the  need  to  consider  both  propagation  and  laser  tissue  interaction  phenomena  for 
ultrashort  laser  retinal  damage.  For  pulses  with  significant  peak  irradiance,  laser  propagation  effects  should  be 
considered  to  determine  the  retinal  spot  characteristics. 
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Figure  3.  Plotted  is  the  peak  power  for  MVL  data  (triangles)  with  an  outline  of  the  peak  power  measured  for  various 
nonlinear  optical  effects  in  the  geometry  of  the  eye  (circle  and  diamond). 
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We  have  developed  an  artificial  eye*°  to  measure  nonlinear  phenomena.  In  Figure  3  the  measured  thresholds  for 
SBS'^,  and  are  shown  as  measured  in  our  artificial  eye.  Because  the  threshold  is  well  above  the  MVL 
threshold  SBS  is  not  expected  to  play  a  role  in  changing  the  MVL  threshold.  Self-focusing  is  known  to  change 
the  spot  size  in  certain  geometries.  Irradiance  enhancement  was  shown  to  occin  above  0.5  MW  of  peak  power 
in  the  geometry  of  the  artificial  eye*^  As  a  result  of  the  measured  MVL  thresholds,  self-focusing  is  not  expected 
to  significantly  affect  the  MVL  threshold  above  several  (tens  of)  picoseconds.  Below  this,  self-focusing  could 
affect  the  retinal  spot  size.  Well  below  several  himdred  femtoseconds  in  pulse  duration  additional  phenomena 
likft  wavelength  dispersion  would  also  influence  the  eventual  beam  characteristics.  Supercontinuum  generation  is 
seen  to  occur  nominally  near  the  MVL  level  for  visible  wavelengths  at  the  shortest  pulsewidths.  At  threshold, 
SC  shows  a  bandwidth  of  ten  nanometers  or  so,  and  probably  would  not  affect  the  MVL  level.  At  higher  energy 
levels,  SC  can  produce  wavelengths  from  UV  to  IR  light  as  the  pulse  broadens  spectrally.  This  deposition  of  UV 
on  the  retina  could  cause  interesting  retinal  effects,  but  the  effect  of  these  short  wavelengths  on  the  retina  have 
not  been  studied. 

The  thresholds  for  laser  induced  breakdown  in  the  artificial  eye  have  been  measured^’’®  and  compared  for  several 
different  focusing  geometries.  We  have  also  measured  the  threshold  for  LIB  creation  in-vivo^"^’*  through  the  use 
of  video  imaging.  Through  these  investigations,  we  have  found  that  LIB  is  a  MVL-threshold-level  damage 
mechanism  for  100  fs  visible  pulses,  but  occurs  at  superthreshold  levels  for  longer  pulse  durations. 

As  displayed  in  Figure  2,  laser  interaction  with  the  retina  plays  an  important  role  in  determining  the  threshold  for 
damage.  For  ultrashort  pulses,  stress  confinement  in  retinal  melanin**  has  been  hypothesized  to  accoimt  for  the 
dropping  MVL  threshold.  This  view  is  different  than  the  trend  seen  in  skin  melanin,  where  the  threshold  for 
immediate  dermal  whitening  was  shown  to  be  nearly  independent  of  pulse  width**  for  pulse  durations  between 
65  fs  and  40  ns.  In  collaboration  with  Charles  Lin  and  his  team  at  Wellman  Laboratories,  we  are  presently  doing 
experiments  to  determine  the  effect  of  smaller  spot  sizes  on  retinal  melanin  disruption  and  the  ultimate  damage 
response  of  retinal  tissues.  We  are  reporting  m  this  year’s  Laser-Tissue  Interaction  DC  the  thresholds  for 
damaging  an  artificial  retina  in  an  artificial  eye  using  both  visible  and  near-inffared  femtosecond  to  nanosecond 
pulses*^.  It  is  expected  that  the  interplay  between  the  propagation  and  tissue  interaction  issues  will  be  critical  in 
explaining  the  drop  in  MVL  threshold  for  ultrashort  pulses. 


3.  SUMMARY 

Self-focusing,  LIB  and  melanin-mediated  effects  are  proposed  as  possible  factors  affecting  sub-nanosecond 
MVL  thresholds.  Most  of  the  previous  research  has  not  considered  the  chromatic  aberration  in  the  eye  where 
NIR  light  is  focused  behind  the  absorbing  layer  of  the  eye  (i.e.  the  retinal  pigmented  epithelial  layer).  To  extend 
these  arguments  into  the  NIR,  one  needs  to  consider  that  light  at  NIR  wavelengths  focus  behind  the  retina*^.  As 
reported  by  Cain  et  al.**,  the  NIR  MVL  data  follow  a  more  dramatic  trend  than  that  seen  in  the  visible  MVL’s. 
The  relevant  mechanism(s)  responsible  for  lowering  the  MVL  threshold  will  also  need  to  explain  the  enhanced 
drop  for  the  NIR  wavelength  regime. 

There  are  numerous  implications  for  laser  safety  from  this  work.  The  drop  in  MVL’s  below  one  nanosecond 
suggests  that  the  current  constant-fluence  laser  safety  standard*®  (from  one  nanosecond  to  more  than  ten 
microseconds)  is  not  appropriate  for  shorter  pulsewidths.  There  exists  guidance  in  the  current  ANSI  standard*® 
that  recommends  using  a  constant  power  below  one  nanosecond.  Our  data  suggests  that  this  is  prohibitively 
conservative.  We  have  suggested  a  more  lenient  standard  for  ultrashort  laser  pulses  and  expect  that  in  the  next 
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ANSI  rewrite  the  maximum  permissible  exposure  limits  for  ultrashort  laser  pulses  will  be  changed  to  reflect  our 
data  and  analysis. 


The  use  of  ultrashort  laser  pulses  offers  the  possibility  of  unique  solutions  to  industrial  problems  such  as  remote 
sensing,  high  precision  LEDAR  and  laser  welding.  The  increased  interest  in  using  ultrashort  laser  pulses 
highlights  the  current  deficiencies  in  the  national  laser  safety  standard.  Our  work  is  intended  to  define  the 
interplay  between  the  nonlinear-optical  and  melanin-mediated  effects  that  determine  the  eventual  MVL 
threshold. 
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ABSTRACT 

We  characterized  the  effects  of  pulse  duration,  pulse  energy,  and  spot  separation  on  intrastromal  comeal 
photodisraption  to  determine  parameters  that  achieve  optimal  surface  quality  and  tissue  plane  separation. 
Experiments  utilized  two  laser  systems,  a  60  picosecond  Nd:YLF  laser  and  a  450  femtosecond  NdrGlass 
laser,  both  operating  at  1 .06  pm  wavelength.  Photodisruption  was  performed  by  tightly  focusing  the  laser 
beam  150  microns  below  the  tissue  surface  and  scanning  it  in  a  spiral  pattern  to  create  a  plane.  A  cut  to  the 
surface  was  made  with  the  laser  and  the  two  surfaces  separated  to  form  a  flap.  Tissue  plane  separation  was 
graded  according  to  the  additional  mechanical  dissection  required.  Internal  surfaces  were  analyzed  with 
standard  histologic  methods  and  scanning  electron  microscopy.  We  found  that  the  Nd:YLF  laser  required 
approximately  three  times  the  pulse  energy  to  achieve  intrastromal  cuts.  Picosecond  parameters  also 
required  more  mechanical  dissection  and  produced  lower  surface  quality  than  optimal  femtosecond 
parameters.  We  conclude  that  femtosecond  laser  pulses  offer  significant  advantages  that  make  them  ideal 
candidate  tools  for  high  precision  intrastromal  comeal  surgery.  The  flexibility  in  laser  pulse  delivery  opens 
up  a  number  of  potential  surgical  applications  not  possible  with  current  mechanical  or  laser  devices. 

Keywords:  cornea,  photodisruption,  ultrafast  laser,  femtosecond,  surgery 

2.  INTRODUCTION 

In  contrast  to  direct  photoablation  with  ultraviolet  light,  such  as  in  excimer  laser  photorefractive 
keratectomy  (PRK), '  comeal  photodisraption  with  transparent  wavelengths  is  based  on  laser-induced  optical 
breakdown  (LIOB).  LIOB  occurs  when  the  fluence  (energy/area)  at  the  laser  focus  reaches  a  threshold  that 
transforms  matter  of  normal  state  to  one  of  a  high  density  of  free  electrons  and  ions,  called  a  plasma. 
Absorption  of  laser  pulse  energy  by  the  opaque  plasma  leads  to  rapid  temperature  and  pressure  increases  in 
the  focal  spot.  The  expansion  of  the  hot  plasma  causes  a  shock  wave  and  development  of  a  cavitation 
bubble.  As  a  result  of  this  process,  termed  photodisraption,  tissue  in  the  focal  volume  is  destroyed.  ^ 

Pulsed  lasers  attain  the  high  field  intensities  required  for  LIOB  and  photodisraption  by  concentrating  light 
energy  into  a  short  time  duration  (since  field  intensity  is  inversely  proportional  to  the  pulse  duration). 
Historically,  ophthalmic  surgical  procedures  which  utilize  photodisraption  have  relied  on  the  Nd:YAG  laser, 
which  produces  pulses  of  approximately  10  nanosecond  and  1 .06  pm  wavelength.  At  this  pulse  duration, 
the  energy  required  for  beginning  LIOB  is  relatively  high,  which  then  is  transmitted  to  the  surrounding 
tissue  in  the  form  of  acoustic  (shock)  effects.  The  relatively  large  secondary  shock  waves  associated  with 
the  Nd;YAG  laser  limit  its  ophthalmic  applications  to  procedures  such  as  capsulotomy  and  iridotomy.  ^ 

Previous  studies  have  shown  that  lasers  operating  at  shorter  pulse  durations  require  less  energy  and  produce 
smaller  shock  waves  and  cavitation  bubbles  than  the  nanosecond  Nd:YAG  laser®"'®  The  Nd;YLF  laser, 
which  produces  pulses  in  the  30-60  picosecond  (picosecond)  range,  also  operates  in  the  near-infrared 
wavelength  (1.06  microns).  Due  to  its  more  localized  collateral  effects,  the  Nd:YLF  pulses  can  be  focused 
within  the  cornea  to  create  intrastromal  photodisraption.  By  scanning  the  laser  focus  at  a  high  repetition 
rate,  adjacent  photodisraptions  can  perform  near-contiguous  cutting. 
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The  Nd:  YLF  laser  has  been  evaluated  clinically  for  a  variety  of  refractive  and  therapeutic  corneal 
applications.  Despite  attempts  at  optimizing  laser  and  delivery  system  parameters,  including  the  pulse 
energy,  spot  size,  and  spot  separation,  picosecond  laser  intrastromal  procedures  have  had  limited  clinical 
success  due  to  an  inability  to  produce  fully  contiguous  intrastromal  photodisruption.^'  ^ 

Femtosecond  laser  pulses,  which  have  a  pulse  duration  approximately  100  times  shorter  than  Nd:YLF 
picosecond  pulses,  may  offer  improvements  that  could  make  intrastromal  photodisruption  clinically 
practical.  Compared  to  picosecond  pulses,  LIOB  energy  thresholds,  shock  waves  and  cavitation  bubbles  for 
femtosecond  pulses  are  significantly  reduced.  At  these  pulse  durations,  the  volume  of  tissue  affected  by  the 
secondary  shock  wave  and  cavitation  bubble  are  on  the  order  of  the  spot  size,  possibly  allowing  the 
contiguous  photodisruption  required  in  clinical  procedures*^ 

To  evaluate  femtosecond  laser  corneal  lamellar  dissections  and  surface  quality,  we  performed  patterns  of 
intrastromal  corneal  photodisruption  in  animal  cadaver  eyes,  using  a  range  of  laser  energies  and  spot 
separations.  Results  from  a  prototype  solid-state  femtosecond  laser  were  compared  with  similar  procedures 
performed  with  a  Nd:  YLF  picosecond  laser. 


3.  METHODS 


3.1  Cadaver  eyes 

Intrastromal  photodisruption  was  performed  on  porcine  and  primate  cadaver  eyes.  Porcine  eyes  were  obtained 
from  a  local  slaughterhouse  and  used  within  four  hours  of  death.  They  were  kept  on  ice  to  maintain  corneal 
clarity  and  thickness,  which  varied  between  750  and  850  p,m(Opthascan,  Teknar  Inc).  Primate  eyes  were 
enucleated  after  euthanasia,  performed  as  part  of  an  unrelated  experimental  protocol.  They  were  also  kept  on 
ice  and  used  within  four  hours  of  death.  Comeal  thickness  varied  between  400  and  500  pm. 

3.2  Lasers 

Femtosecond  pulses  were  produced  by  a  solid-state,  Nd:Glass  laser,  operating  at  1.06  mm  wavelength  and 
450  femtosecond  pulse  duration.  The  system  consists  of  a  laser  oscillator,  a  pulse  stretcher,  a  regenerative 
amplifier,  and  a  pulse  compressor.  The  oscillator  is  mode-locked  by  a  solid-state  saturable  absorber  and 
produces  a  train  of  250  femtosecond  pulses  at  a  repetition  rate  of  120  MHz.  Pulses  from  the  output  of  the 
oscillator  are  selected  at  a  few  KHz  repetition  rate  and  stretched  to  a  60  picosecond  pulse  duration.  The 
pulse  is  then  amplified  in  a  diode-pumped  regenerative  amplifier  and  then  recompressed  to  450  femtosecond 
pulse  duration.  Laser  energy  could  be  varied  from  1-25  pj.  The  beam  produced  from  the  regenerative 
amplifier  had  fundamental  Guassian  (TEM_00)  spatial  profile.  For  comparisons  with  picosecond  pulses,  a 
60  picosecond  Nd:YLF  laser  (Model  5000,  Escalon  Medical  Corp,  Skillman,  NJ)  was  passed  through  the 
same  delivery  system.  Pulse  energy  for  this  system  could  be  varied  from  5-50  pj. 

3.3  Corneal  delivery  system 

The  beam  is  passed  through  a  computer-controlled  comeal  delivery  system,  which  can  focus  to  a  spot  size 
less  than  10  pm  within  the  comeal  stroma.  The  focal  spot  can  be  scanned  within  a  disc  of  6  mm  and  a 
height  of  1  mm.  The  eye-laser  interface  is  a  piano  contact  lens  applanation  system,  with  globe 
stabilization  provided  by  a  suction  ring.  The  depth  of  the  scanned  pattern  is  measured  relative  to  the  lower 
surface  of  the  contact  lens  to  maintain  accuracy. 

3.4  Surgical  procedure 

Globes  were  placed  in  a  custom-built  suction  holder  which  maintained  an  intraocular  pressure  of 
approximately  30  mm  Hg.  The  modified  contact  lens  and  laser  delivery  system  was  then  applied  to  the 
comeal  surface.  A  spiral  pattern  of  laser  pulses,  focused  to  a  spot  size  of  8  pm,  150-200  pm  below  the 
comeal  surface,  was  performed  at  diameters  ranging  from  3-6  mm  (Figure  1).  Spot  separation  varied  from 
10-40  pm.  After  creating  a  subsurface  plane,  a  cut  to  the  surface  was  made  by  advancing  a  circle  pattern 
from  the  plane  through  the  epithelium  to  create  a  comeal  flap  (Figure  2). 
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Figure  1:  A  schematic  of  an  intrastromal  femtosecond  laser  cut.  (a)  The  front  view  of  the  intrastromal 
spiral  disk  pattern,  (b)  The  side  view  of  the  spiral  pattern  within  the  cornea.  Note  that  the  size  of  the 
bubbles  is  about  120  times  smaller  than  the  thickness  of  the  cornea. 


Figure  2:  A  schematic  of  the  partially  lifted  comeal  flap. 

3.5  Analysis 

Lamellar  dissections  were  graded  using  a  one  to  four  scale,  as  follows; 

Grade  1:  Lamellar  dissection  requires  extensive  free  blade  cutting 

Grade  2:  Free  blade  cutting  required  at  leading  edge  only  to  separate  surfaces 

Grade  3:  Significant  counter-traction  needed  to  break  residual  adhesions,  but  no  blade  required  to  cut 

adhesions  between  surfaces 

Grade  4;  Surfaces  separate  with  gentle  traction  on  anterior  surface  using  forceps  only. 

Following  the  procedure,  globes  were  fixed  in  2%  glutaraldehyde  with  the  photodisrupted  surfaces  exposed. 
Internal  surfaces  were  analyzed  by  scanning  electron  microscopy  using  standard  techniques. 


Orading 


4.  RESULTS 


4.1  Porcine  lamellar  dissections 

We  first  investigated  the  ability  of  the  femtosecond  laser  pulses  to  create  lamellar  dissections  that  did  not 
require  mechanical  assistance.  Based  on  our  previous  studies  evaluating  femtosecond  cavitation  bubble  size 
as  a  function  of  laser  pulse  energy,  values  from  2-8  pj  were  tested  at  spot  separations  ranging  from  10-30 
microns  (spot  size  was  constant  at  approximately  8  pm).  Dissection  quality  was  graded  on  one  to  four 
scale  as  described  in  the  methods. 

As  shown  in  Figure  3A,  the  ease  with  which  the  lamellar  dissection  could  be  accomplished  was  related  both 
to  the  pulse  energy  and  spot  separation.  For  low  pulse  energies  (2  pj  &  4  pj),  best  dissection  was 
accomplished  at  small  spot  separations  (10  pm),  while  for  the  larger  pulse  energy  (8  pj)  wider  spot 
separation  (up  to  20  pm)  still  permitted  excellent  separation. 

These  trends  were  also  identified  with  60  picosecond  laser  pulses  at  energies  from  10-25  pj  and  spot 
separations  from  10-30  pm.  As  seen  in  Figure  3B,  blade-free  dissection  (Grades  3  &4)  was  possible  with 
fewer  parameter  combinations  than  with  the  450  femtosecond  pulses.  Grade  4  dissection  was  accomplished 
only  at  25  pj  and  10-20  pm  spot  separations. 
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ComMi  Flap  Dissection  as  a  Function  of 
Laser  Energy,  Spot  Separation,  and  Pulse  Duration 
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Figure  3:  Dissection  of  corneal  flaps  as  a  function  of  pulse  energy  and  spot  separation  for  450  femtosecond 
(a)  and  60  picosecond  laser  (b)  pulses.  Both  4  and  8  pj  femtosecond  pulses  produced  high  quality 
dissection,  while  25  pj  pulse  energies  were  required  for  picosecond  pulses. 

4.2  SEM  analysis  of  porcine  lamellar  dissections 

Electron  micrographs  comparing  the  surfaces  characterized  by  different  laser  parameter  combinations  revealed 
qualitative  differences,  which  corresponded  to  the  dissection  scale  grades.  Surfaces  created  during  grade  4 
dissections  generally  showed  a  homogenous,  dimpled  pattern  (Figure  4).  Such  surfaces 


50  microns 

Figure  4:  SEMs  of  450  femtosecond  high  grade  (Grade  4)  porcine  cornea  lamellar  dissections  made  with  4 
pj  pulse  energy,  10  pm  spot  separation  (A),  8  pj  pulse  energy,  10  pm  spot  separation  (B),  8  pj  pulse 
energy,  15  pm  spot  separation  (C),  and  8  pj  pulse  energy,  20  pm  spot  separation  (D).  A  more  irregular 
surface  was  seen  with  the  wider  spot  separation  (D),  suggesting  tissue  tearing  during  dissection. 
Magnification  200  X. 

Similar  observations  were  found  for  60  picosecond  pulses  at  the  25  pj  energy.  For  the  10  pm  spot 
separation,  a  surface  pattern  comparable  to  the  best  femtosecond  dissections  was  obtedned.  However,  for  the 
20  pm  spot  separation,  which  also  produced  a  grade  4  dissection,  a  more  irregular  SEM  surface  was  seen 
(Figure  5),  similar  to  the  same  spot  separation  for  450  femtoseconds  at  8  pj  energy  (Figure  4b) 

Surfaces  created  by  femtosecond  and  picosecond  pulses  during  low  (1-2)  grade  dissections  displayed  irregular 
surface  patterns  suggestive  of  tissue  tearing  and  use  of  the  knife  blade  to  complete  the  cuts  (Figure  6).  All 
femtosecond  dissections  were  at  least  grade  2,  while  several  picosecond  parameter  combinations  required 
extensive  mechanical  assistance  (grade  1).  The  blade  cuts  are  more  easily  seen  on  these  grade  1  surfaces, 
although  they  are  still  apparent  on  grade  2  dissections. 


could  be  created  with  femtosecond  pulses  at  both  4  pj  and  8  pj,  with  wider  spot  separations  (up  to  15  pm) 
possible  at  the  higher  energy.  Although  a  dissection  grade  of  4  was  obtained  at  8  pj  and  20  pm  spot 
separation,  the  surface  quality  was  less  homogenous  than  at  the  smaller  spot  separations. 


Figure  5:  SEMs  of  60  picosecond  high  grade  (Grade  4)  porcine  cornea  lamellar  dissections  made  with  25  |J,j 
pulse  energy,  10  |J,m  spot  separation  (A),  25  pj  pulse  energy,  20  pm  spot  separation  (B).  A  more  irregular 
surface  was  seen  with  the  wider  spot  separation  (B),  suggesting  tissue  tearing  during  dissection. 
Magnification  200  X. 
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Figure  6:  SEMs  of  60  picosecond  porcine  cornea  lamellar  dissections  made  with  9.5  pj  pulse  energy,  10 
pm  spot  separation  (A,  grade  1),  16  pj  pulse  energy,  10  p,m  spot  separation  (B,  grade  2),  and  450 
femtosecond  porcine  cornea  lamellar  dissections  made  with  4  pj  pulse  energy,  20  pm  spot  separation  (C, 
grade  2),  and  8  pj  pulse  energy,  30  pm  spot  separation  (D,  grade  2).  Magnification  200  X. 
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4.3  Lenticule  cutting  in  porcine  and  primate  cadaver  eyes 

To  evaluate  a  prototype  refractive  procedure  ex  vivo,  we  utilized  optimized  femtosecond  and  picosecond 
parameters  to  perform  laser  keratomileusis, “  in  which  a  block  of  corneal  tissue  is  removed  after 
photodisrupting  its  perimeter  with  scanned  laser  pulses  (Figure  7). 


Figure  7:  A  schematic  of  laser  keratomileusis,  (a)  A  femtosecond  laser  cuts  a  lenticule  within  the  stroma 
of  the  cornea.  The  laser  is  then  used  to  cut  a  hinged  flap  along  the  perimeter  of  the  lenticule  and  up  the 
surface  of  the  cornea,  (b)  The  flap  is  lifted  and  the  lenticule  is  removed.  Afterward,  the  flap  is  replaced  and 
allowed  to  heal. 

We  first  performed  this  procedure  in  fi'esh  porcine  cadaver  eyes  using  femtosecond  pulses  at  10  ^m  spot 
separation  and  7.5  pJ  energy.  Lenticules  as  thin  as  50  pm  could  be  created  and  removed  with  forceps  alone 
(Figure  8).  Because  similar  results  using  picosecond  pulses  have  not  translated  to  procedures  in  humans, 
we  next  evaluated  the  procedure  in  a  model  more  closely  resembling  the  human  cornea.  Eyes  from 
euthanized  rhesus  monkeys  were  harvested  within  2  hours  of  death  and  subjected  to  laser  keratomileusis. 
Utilizing  similar  parameters  (15  pm  spot  separation  and  7.5  pJ  energy),  successful  procedures  could  be 
performed  using  450  femtosecond  pulses  (Figure  9),  with  excellent  dissection  and  surface  qualities.  In 
contrast,  we  could  not  perform  laser  in  situ  keratomileusis  in  enucleated  primate  eyes  with  picosecond 
pulses  without  extensive  mechanical  dissection,  even  when  parameters  which  had  produced  excellent  results 
in  porcine  eyes  (Figure  10). 
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Figure  8;  SEM  of  femtosecond  laser  keratomileusis  in  enucleated  porcine  eye  at  lOX  (A)  and  200  X  (B) 
magnification.  Note  smooth  surface  made  with  7.5  p.j  pulse  energy,  10  |im  spot  separation.  A  70  |im 
lenticule  is  shown  to  right  of  corneal  hed,  with  flap  to  left. 


Figure  9:  SEM  of  femtosecond  laser  keratomileusis  in  enucleated  primate  eye  at  lOX  (A)  and  200  X  (B) 
magnification.  Note  smooth  surface  made  with  7.5  p.j  pulse  energy,  15  |tm  spot  separation.  A  100  |J.m 
lenticule  is  shown  to  left  of  corneal  hed,  with  flap  to  right. 


Figure  10:  SEM  of  60  picosecond  laser  keratomileusis  in  enucleated  primate  eye  at  lOX  (A)  and  200  X  (B) 
magnification.  Note  irregular  surface  made  with  25  pj  pulse  energy,  15  pm  spot  separation.  A  lenticule 
could  not  be  dissected  free  due  to  interconnections  between  the  anterior  lenticule  surface  and  the  flap. 

5.  DISCUSSION 

Lamellar  corneal  dissection  (using  scanned  patterns  of  ultrashort  laser  intrastromal  photodisruption)  has 
been  advocated  as  a  potential  adjunct  and  alternative  to  current  refractive  surgery  techniques.’®’  While 
picosecond  lasers  were  first  advanced  for  this  purpose,  clinical  results  in  humans  have  been  disappointing. 
Our  findings  indicate  that  many  of  the  shortcomings  of  picosecond  technology  (including  incomplete  tissue 
cutting,  poor  dissection  and  surface  quality)  are  significantly  improved  when  the  pulse  duration  is  further 
reduced  into  the  femtosecond  regime. 

The  clinical  improvements  obtained  with  femtosecond  pulses  can  be  understood  by  considering  the  pulse 
duration  dependency  of  the  photodisruptive  process  itself.  Complete  intrastromal  lamellar  dissection 
requires  contiguous  photodisruption,  in  turn  necessitating  a  match  between  the  plasma  volume,  cavitation 
bubble  size,  and  the  pattern  spot  separation.  If  the  spot  separation  in  a  particular  scanning  pattern  is  too 
large  for  the  chosen  laser  parameters,  tissue  islands  are  left  behind  that  must  be  broken  mechanically.  If 
spot  separations  are  too  small,  the  cavitation  bubble  interferes  with  the  next  photodisruption,  thereby  also 
preventing  contiguous  cutting. 

Photodisruption  begins  with  formation  of  a  plasma,  the  size  of  which  is  largely  dependent  on  the  focal  spot 
diameter  of  the  optical  delivery  system.  Generally,  the  smallest  spot  diameter  is  desired,  since  this  requires 
the  lowest  pulse  energy  to  exceed  the  intensity  threshold  for  plasma  generation.  In  most  delivery  systems, 
there  is  a  practical  optical  limitation  on  spot  diameter  of  8-10  |im,  which  is  constant  over  a  wide  range  of 
pulse  durations  and  energies.  Procedure  time  and  mechanical  limitations  also  impose  constraints  on  how 
small  a  focus  size  can  be  used  and  how  fast  it  can  be  scanned. 

At  these  8-10  pm  focal  spot  sizes,  nanosecond  pulses  from  the  NdrYAG  laser  require  millijoule  level 
energies  to  generate  a  plasma.®  The  secondary  acoustic  effects  spread  over  several  hundred  microns,  making 
contiguous  corneal  photodisruption  impossible. For  picosecond  pulses  operating  at  threshold  energies, 
focal  spots  of  approximately  8-10  pm  predict  final  cavitation  bubble  sizes  of  approximately  25  pm,  a  1:3 
ratio  which  results  in  near-contiguous  photodisruption.  In  contrast,  for  femtosecond  plasma  volumes  of 
approximately  8-10  pm  at  threshold  energies,  cavitation  sizes  are  predicted  at  approximately  half  picosecond 
values,  at  approximately  12  pm.’"*  This  ratio  of  1:1.5  ,  translates  into  closer  placement  of  femtosecond 


thresholds  associated  with  femtosecond  pulses,  which  translate  into  highly  controlled  and  reproducible 
photodisruption.  ® 

While  the  current  results  in  porcine  eyes  are  encouraging,  the  differences  in  picosecond  laser  behavior  in  this 
model  (as  reported  previously)  and  human  clinical  studies  must  be  kept  in  mind.  Our  findings  also 
reveal  that  picosecond  parameters  that  work  in  enucleated  porcine  corneas  do  not  translate  well  into  (freshly 
enucleated)  primate  eyes.  While  femtosecond  parameters  did  cross  over  well  to  the  ex  vivo  primate  model, 
in  vivo  systems  that  closely  reproduce  the  clinical  situation  must  ultimately  be  used. 

The  laser  parameters  outlined  here  provide  attainable  requirements  for  any  potential  clinical  devices. 
Advances  in  ultrafast  laser  design  and  the  development  of  powerful  laser  diodes  has  made  low  cost,  reliable, 
diode-pumped,  femtosecond  laser  systems  possible  for  the  first  time.  Combined  with  sophisticated 
scanning  delivery  systems  now  available,  a  number  of  potential  laser  corneal  procedures  are  possible, 
including  corneal  flap  cutting,  laser  keratomileusis  and  possibly  intrastromal  keratectomy. 
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ABSTRACT 

We  determined  the  wavelength  dependence  of  the  minimum  spot  size  of  a  laser  beam  focused  through  human 
sclera  to  evaluate  the  potential  for  transcleral  glaucoma  surgical  techniques  using  ultrashort-pulsed  lasers.  The 
spectrum  of  the  forward  scattered  light  was  measured  by  collimating  the  incident  and  transmitted  beam  in  a 
spectrophotometer.  This  spectrum  shows  that  sclera  is  highly  scattering  until  1100  nm,  after  which,  the 
transmission  spectrum  is  similar  to  water.  To  measure  the  minimal  spot  size,  a  laser  beam  was  focused  on  the 
back  surface  of  sclera  of  differing  thickness.  The  minimum  spot  at  800  nm,  1060  nm,  1301  nm,  and  1557  nm 
was  imaged.  At  800  nm,  the  spot  size  was  invariant  upon  focal  lens  position,  being  a  thousand  fold  larger  than  ftie 
incident  beam  spot  size.  As  the  wavelength  increased,  the  area  of  the  spot  decreased,  so  that  at  1557  nm,  the 
minimal  spot  size  was  on  the  order  of  the  incident  beam  spot  size. 

Keywords:  sclera,  scattering,  spot  size,  wavelength  dependence,  surgery,  glaucoma 


2.  INTRODUCTION 

Glaucoma  affects  ^proximately  two  million  Americans  while  another  ten  million  are  at  risk  due  to 
statistically  significant  elevated  intraocular  pressure  (lOP).  Medications  can  lower  lOP,  although  they  are 
only  partially  effective,  with  an  approximately  40%  ftulure  rate  over  3  years.’  Laser  and  incisional  surgical 
procedures,  which  increase  fluid  drainage  from  the  eye,  are  generally  reserved  for  medication  failures. 
Argon  laser  trabeculoplasty,  the  most  common  laser  glaucoma  procedme,  may  stretch  the  normal  channels 
in  die  trabecular  meshwork,  thereby  increasing  outflow  of  aqueous  humor.  The  effectiveness  of  this 
procedure  is  short  lived  however,  averaging  about  two  years,  after  which  incisional  surgery  is  often  tried.*’’ 

Filtering  surgery,  the  cutting  of  a  drainage  channel  (fistula)  directly  into  the  sclera,  is  the  most  common 
incisional  surgery  for  the  treatment  of  glaucoma.  The  various  techniques  all  share  a  host  of  complications, 
including  very  low  intraocular  pressure  (hypotony),  cataract  formation,  and  infection.^’’’®  In  addition, 
surgical  failure  often  occurs  over  the  course  of  mondis  to  years  due  to  scarring  and  closure  of  the  fistula. 
Despite  these  problems,  glaucoma  specialists  are  moving  towards  earlier  surgical  intervention  due  to  its 
better  lOP  control  in  most  patients,  increasing  the  need  for  a  safer  and  more  effective  surgical  alternative. 

The  majority  of  post-operative  scarring  responsible  for  fistula  closure  occurs  due  to  wound  healing  in  the 
tissues  overlying  the  sclera,  the  conjunctiva  and  episclera,  which  are  damaged  during  incisional 
procedures.*  A  number  of  laser  surgical  procedures  to  improve  surgical  outcomes  have  been  proposed  and 
tested.  *  These  generally  fall  into  two  approaches: 
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1 .  Ab  intemo  procedures,  which  ablate  tissue  from  the  internal  surface  of  the  sclera,  either  by 
directing  the  light  through  the  cornea  with  a  contact  lens  or  directly  via  a  probe  introduced 
into  the  anterior  chamber  of  the  eye. 

2.  Ab  extemo  procedures,  which  ablate  tissue  from  the  external  surface  of  the  sclera,  either  by 
directing  the  light  through  the  conjunctiva  or  directly  via  a  probe  introduced  under  the 
conjunctiva. 

Despite  these  attempts,  no  laser-based  method  has  delivered  enough  clinical  or  economic  advantages  to 
displace  traditional  surgical  methods.  Ultrashort-pulsed  lasers,  operating  at  appropriate  wavelengths,  may 
permit  unique  transcleral  procedures  that  could  offer  advantages  of  each  of  these  methods. 

Femtosecond  laser  pulses  used  for  cornea  surgery  have  been  demonstrated  to  cause  little  collateral  damage, 
both  due  to  the  low  average  power  ~  no  burning  —  and  to  their  low  pulse  energy  -  small  shock  waves  and 
bubble  oscillations.  ®  In  addition,  subsurface  tissue  photodisruption  is  possible  and  has  been  demonstrated 
in  transparent  comeal  tissue  using  near-infrared  and  infrared  wavelengths.  By  using  a  wavelegnth  that 
is  transmitted  by  the  sclera,  and  is  able  to  be  focused  at  its  back  surface,  it  may  be  possible  to  create  an  ab 
intemo  fistula  without  the  technical  problems  associated  with  contact  lens  or  probe  delivery  methods.  Such 
a  method  would  avoid  disraption  of  the  overlying  conjuctiva  and  episclera,  something  current  ab  extemo 
techniques  cannot  offer. 

As  a  first  step,  this  paper  examines  the  minimum  spot  size  attainable  as  a  function  of  wavelength  when  the 
laser  is  focussed  through  the  sclera  itself  to  its  back  surface.  Our  first  experiment  yielded  forward 
transmission  for  the  wavelengths  between  500  nm  and  2500  nm  using  a  spectrophotometer  to  obtain  an 
initial  indicator  of  scattering  and  minimum  spot  size.  For  the  second  experiment,  we  imaged  the  spot  on 
the  back  surface  of  the  sclera  using  a  vidicon  camera  with  a  microscope  objective,  comparing  wavelengths 
from  800  nm  to  1557  nm. 


3.  METHODS 


3.1  Tissue  Samples 

Human  scleral  sections  (not  suitable  for  transplantation  )  from  regions  adjacent  to  the  cornea  were  cut  to  a 
thickness  ranging  from  0.30  mm  to  0.77  mm  using  a  vibratome.  Four  different  thickness  sections  from 
each  of  three  globes  were  cut  for  a  total  of  twelve  samples.  The  sample  thickness  was  measured  before 
each  experiment  using  a  calipers  and  all  measiuements  were  made  within  two  weeks  of  globe  donation. 

3.2  Experiment  1:  Spectrophotometer  Measurements 

The  forward  scattered  light  as  a  function  of  wavelengtii  was  measured  using  a  spectrophotometer  (Perkin- 
Elmer  Lambda  9  UVAHS/NIR  Spectrophotometer)  by  placing  the  apparatus  shown  in  Fig.  1  inside  the 
sample  chamber.'^”  The  apparatus  collimates  the  incident  incoherent  broadband  light  using  two  irises  with 
a  2.16  mm  diameter  opening  separated  by  approximately  3.0  cm.  The  tissue  was  mounted  to  a  slide  (1mm 
thick)  using  a  drop  of  saline  solution.  The  forward  scattered  light  was  collected  using  two  irises  separated 
by  3  cm  with  2.16  mm  opening.  The  tissue  was  centered  between  the  two  sets  of  irises  separated  by 
approximately  2.3  cm. 

For  each  of  the  twelve  samples,  three  measurements  of  the  absorption  spectrum,  500-2500  nm,  were  taken, 
then  converted  to  total  transmission.  This  measurement  was  taken  to  identify  trends  in  forward  scattering 
to  narrow  the  search  for  optimal  wavelengths  to  penetrate  the  sclera. 


Slide 


Figure  1:  The  tissue  holder  for  the  spectrophotometer  only  allows  collimated 
light  to  strike  the  sclera  and  only  collects  the  forward  transmitted  light. 
The  incident  beam  is  a  converging  incoherent  broadband  source. 


3.3  Experiment  2:  Spot  Size  Measurements 

Figure  2  shows  the  experimental  setup  used  to  measure  the  minimum  spot  size.  The  laser  light  is  incident 
from  the  left  onto  a  variable  attenuator  to  prevent  camera  saturation.  The  beam  then  passes  through  a  soft 
aperture  (about  2  m  of  Coming  Flexcor  1060  or  Newport  F-SV  620  for  800  nm  optical  fiber)  to  insure  that 
the  beam  from  each  source  is  ftmdamental  Gaussian  mode.  The  emergent  beam  was  collimated  (Newport 
1015LD,  f  =  7.4  mm),  then  focused  using  a  6.24  nun  focal  length  lens  with  a  NA  of  0.40  (Thor  Labs 
3501 10-C).  The  front  surface  of  a  0.15  mm  thick  microscope  cover  slip  was  placed  at  the  focal  spot.  This 
surface  corresponds  to  the  back  surface  of  the  sclera.  The  vidicon  camera  (Electrophysics  Model  7290A) 
was  positioned  such  that  the  object  plane  corresponded  to  the  back  surface  of  the  tissue  while  using  a  20X 
0.40  NA  microscope  objective.  The  minimum  spot  size  was  recorded  as  a  640  X  480  24  Bit  RGB  bitmap. 
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Figure  2:  The  minimum  spot  size  was  determined  by  placing  the  back  surface  of  the  sclera  in  the  focal  plane  of  a 
microscope  camera.  The  focusing  lens  was  moved  until  the  minimum  spot  size  was  obtained. 


The  tissue  was  placed  on  the  cover  with  a  drop  of  saline.  The  lens  position  that  produced  the  minimum 
spot  size  at  the  back  surface  of  the  sclera  was  determined  by  shaking  the  tissue  horizontally  in  the  focal 
plane  and  moving  the  lens  until  the  smallest  spot  occurred.  The  Electrophysics  Model  7290A  camera  has  a 
slow  response  time  so  that  shaking  the  sample  allowed  us  to  average  over  a  portion  of  the  tissue. 

The  shaker  was  DC  motor  with  an  off  center  circular  cam  that  provided  about  1  nun  of  motion.  This  image 
was  recorded,  then  the  shaker  was  stopped  and  several  stationary  images  of  the  speckle  were  taken.  The 
above  procedure  was  preformed  using  four  different  lasers:  a  800  nm  Ti:Sapphire  oscillator  operating  at 
100  hfflz  and  50  fs  pulses,  a  1060  nm  Nd:Glass  oscillator  at  100  MHz  and  300  fs  pulses,  a  1301  nm 
continuous  wave  laser  diode,  and  a  Er:Glass  fiber  laser  operating  with  100  fs  pulses  at  100  MHz. 
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3.4  Image  Analysis 


The  camera  images  collected  were  corrected  for  the  nonlinear  camera  response  and  normalized  such  that 
the  peak  intensity  was  the  same  for  all  images.  The  fiill-area  half-maximum  (FAHM)  sizes,  as  defined  as 
the  area  of  the  spot  that  is  greater  than  or  equal  to  the  maximum  pixel  intensity,  were  calculated  for  each 
image  as  a  comparison  between  the  different  wavelengths.  The  FAHM  is  the  relevant  measure  since  the 
area  of  the  spot  is  directly  proportional  to  the  energy  for  photodisruption. 


4.  RESULTS 


4.1  Forward  Scattering  of  Sclera 

Twelve  absorption  spectra  of  half  thickness  sclera  (average  0.34  mm)  were  averaged  and  are  shown  in  Fig. 
3.  Veiy  little  forward  transmission  occurs  below  1 100  nm.  The  graph  shows  transmission  dips  that  follow 
the  water  dips  at  1450  nm  and  1900  nm.  Several  wavelengths  of  interest  are  listed  in  Table  1. 


Transmission  Spectrum  of  HaH  Thickness  Sclera 


Figure  3:  The  forward  transmission  of  sclera  is  the  jagged  curve.  The  smooth  curve  is  the  percent  transmission 
through  0.78  mm  of  saline  solution.  See  Table  1  for  the  numbers  noted  in  the  graph. 


Table  1:  Forward  transmission  of  half  thickness  of  sclera. 


Number 

Wavelength 

(nm) 

Laser 

Forward 

Transmission 

(%) 

Saline 

Transmission 

(%) 

1 

800 

TirSapphire 

2.15 

100 

2 

1060 

Nd:Glass 

4.02 

100 

3 

1301 

Laser  Diode 

6.22 

100 

4 

1557 

EnGlass 

6.77 

59.5 

5 

1639 

NA 

8.96 

76.3 

6 

1695 

NA 

11.02 

77.9 

7 

1777 

NA 

12.11 

63.5 

8 

1845 

NA 

12.14 

55.3 

9 

2229 

NA 

9.01 

31.8 

The  trends  in  spectrophotometer  measurements  were  consistent  between  samples.  Sclera  displayed  the  same 
transmission  minima  as  saline  at  1450  nm  and  1900  nm  (as  can  be  expected  since  sclera  is  80%  water). 

However,  at  wavelengths  below  1100  nm,  where  water  is  transparent  at  this  thickness,  the  forward  transmission 
of  sclera  increased  linearly,  in  contrast  to  expected  results  and  implying  that  signiflcant  scattering  occurs  for 
these  shorter  wavelengths.  At  longer  wavelengths,  the  forward  transmission  increased  except  for  the 
absorption  peaks  of  water  at  1450nm  and  1900  mn.  Local  maxima  occur  at  1695  nm,  1777  nm,  1845  nm,  and 
2229  nm.  The  first  and  last  of  these  peaks  correspond  to  the  water  transmission  peaks  as  shown  in  Fig.  3. 

4.2  Minimum  Spot  Size  Measurements 

The  minimum  size  (defined  as  the  spot  size  without  the  tissue)  are  shown  in  Figure  4.  The  difference  in  sizes  is 
most  likely  due  to  slightly  different  beam  divergences  and  the  wavelength.  The  size  of  the  speckle  image  at  800 
nm  was  invariant  with  translation  of  the  lens.  Figure  5  shows  the  minimum  spot  size  obtainable  at  different 
wavelengths  through  full  thickness  of  sclera.  As  expected  in  highly  scattering  tissue,  the  pattern  is  dominated  by 
speckle.  At  800  nm,  the  average  spot  size  is  roughly  1000  times  the  size  of  the  minimum  spot.  The  average 
sizes  appear  to  decrease  for  1301  nm  and  1557  nm,  respectively. .  Figures  6  display  the  moving  averaged 
images  of  the  spot  size  through  the  full  thickness  of  the  sclera. 
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Figure  4:  The  minimum  spot  size  at  the  tested  wavelengths  are  shown.  The  minimum  spot  size  is  the  size  of  the 
beam  when  focused  on  to  the  fi-ont  surface  of  the  cover  slip. 
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Typical  speckle  patterns  through  full  thickness  of  sclera. 
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Figure  7:  The  FAHM  for  the  averaged  measurements  (shaken)  and  the  speckle  (stationary)  patterns  are  shown.  The 
FAHM  were  normalized  to  the  FAHM  of  the  minimum  spot  size  without  the  tissue. 


Figure  7  shows  the  FAHM  of  the  shaken  tissue.  The  FAHM  of  the  speckle  can  be  less  than  the  minimum 
spot  size.  The  data  at  800  mn  suggests  severely  peaked  and  separated  speckle  since  the  FAHM  of  the 
speckle  is  up  to  two  orders  of  magnitude  less  than  the  average  (Figure  5).  As  the  wavelength  increases,  the 
average  FAHM  becomes  more  predictable,  i.e.  increases  with  increasing  tissue  thickness,  and  the  average 
FAMH  and  speckle  FAHM  begin  to  converge. 

4.3  Prototype  Trans-tissue  Procedure 

To  demonstrate  how  a  transcleral  channel  might  be  created  using  subsurface  photodisruption,  we  next 
performed  a  transcomeal  fistula  using  a  1060  mn  450  fs  laser.  The  method  is  schematized  in  Figure  8a.  In 
this  study,  we  left  the  channel  partial  thickness  to  show  that  no  damage  to  the  external  surface  was  created. 
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A  similar  experiment  in  sclera  using  the  seime  laser  was  able  to  obtain  consistent  photodisruption  at  depths 
of 250  pm  or  less.  This  corresponds  with  the  average  FAHM  of  Figure  7. 


Figure  8:  (a)  Procedure  used  to  cut  the  cornea  fistula,  (b)  A  hole  cut  in  cornea  from  the  back  surface  to  the  front 
surface  using  1060  nm  femtosecond  pulses.  Note  that  the  top  (external  surface)  layer  is  intact. 


5.  DISCUSSION 

Analogous  to  recent  work  using  ultrashort-pulsed  lasers  for  subsurface  comeal  photodismption,  we  have 
begun  to  characterize  laser  parameters  for  transcleral  surgical  procedures.  Beginning  with  laser 
wavelength,  we  first  defined  the  scattering  properties  of  the  scleral  tissue.  In  contrast  to  the  cornea,  which 
transmits  up  to  90%  of  visible  wavelengths,  the  white  sclera  is  an  excellent  scatterer  of  visible  light. 
Surprisingly,  the  cornea  and  sclera  a  similar  in  structure.  Both  are  composed  of  collagen  fibers  (optical 
index  n  =  1.47)  embedded  in  a  muccopolysaccharide  ground  substance  (n  =  1.345).  The  main  differences 
are  the  average  center  to  center  fibril  spacing  and  the  regularity  of  the  stracture.  Fourier  analysis  reveals 
that  the  center-to-center  spacing  in  cornea  and  sclera  are  59  nm  and  285  nm,  respectively,  with  the 
variations  being  greater  in  sclera.’^'”  Bragg  diffraction  theory  predicts  that  scattering  begins  once  the  half 
wavelength  of  light  in  the  tissue  is  on  the  order  of  the  index  variations  in  the  medium.'®  For  cornea  and 
sclera,  scattering  is  expected  to  be  significant  for  wavelengths  shorter  than  166  nm  and  800  nm, 
respectively,  using  this  theory.  Since  the  peaks  of  the  spatial  frequency  spectrum  of  sclera  are  quite 
broad,'®  significant  scattering  at  longer  wavelengths  may  also  occur. 

Since  the  average  spot  size  decreases  by  about  two  orders  of  magnitude  from  800  run  to  1557  run,  longer 
wavelengths  may  be  perferable  for  transcleral  applications.  Considering  only  scattering,  subsurface 
photodisruption  at  800  nm  would  require  at  least  energies  1000  times  greater  than  surface  photodisruption. 
Pulses  of  such  high  energy  may  reach  threshold  above  the  focal  spot,  resulting  in  photodismption  near  the 
surface  of  the  tissue.  At  1557  nm,  the  pulse  energy  needs  to  be  only  about  four  times  greater  than  at  the 
surface.  These  results  require  expansion  in  ex  vivo  and  in  vivo  models  to  determine  laser  energy, 
repetition  rate  and  pulse  width  parameter  specifications.  The  preliminary  data  in  this  paper  suggests  that 
transcleral  glaucoma  surgery  may  be  possible  at  longer  wavelengths  and  ultrashort  pulse  widths. 
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ABSTRACT 

Retinal  lesions  produced  by  ultrashort  laser  pulses  in  the  pico-  and  femtosecond  range  were  examined  by  electron  microscopy. 
Retinal  pigment  epithelial  (RPE)  cells  that  contained  fractured  and  striated  melanosomes  typically  exhibited  severe  damage  to 
the  other  components  of  the  cell.  However,  having  observed  RPE  cell  damage  without  coincident  fractured  melanosomes,  it 
is  thought  that  melanosome  fracture  itself  is  not  responsible  for  the  damage  that  occurs  within  the  RPE  cell.  Nevertheless, 
the  percentage  of  melanosomes  fractured  per  lesion  seems  to  parallel  the  severity  of  damage  within  that  lesion  site.  No  trend 
existed  between  percentage  of  melanosomes  fractured  and  the  peak  power  of  laser  delivery.  However,  with  decreasing  laser 
pulsewidth,  there  was  a  decrease  in  the  percentage  of  melanosomes  showing  fracture. 

Keywords:  retinal  pigment  epithelium,  laser,  picosecond,  femtosecond,  laser  induced  breakdown,  electron  microscopy, 
Brach’s  membrane 


1.  INTRODUCTION 

As  the  prevalence  of  ultrashort  laser  systems  increases  in  research  and  industry,  safety  standards  must  be  established 
based  upon  a  clear  understanding  of  their  potential  for  harmful  intaaction  with  human  tissue.  Furthermore,  by  clarifying  the 
mechanisms  of  damage  through  with  these  technologically  advanced  laser  systems  act ,  one  may  discover  potentially 
therapreutic  appUcations. 


2.  PREVIOUS  STUDIES 

In  1987,  Bimgruber  et  al.'  concluded  that,  in  a  manner  similar  to  previous  longer  wavelength  systems,  the  energy  of 
ultrashort  laser  pulses  was  absorbed  in  large  part  by  the  melanin  pigment  contained  within  the  retinal  pigment  epithelium. 
After  producing  visible  lesions  in  grey  chinchilla  rabbits,  they  found  that  no  lesions  could  be  detected  ophthalmoscopically  or 
with  fluorescein  angiogrr^hy  when  applying  ultrashort  pulses  to  the  eyes  of  nonpigmented  albino  rabbits.  Goldman  et  al.^ 
conducted  studies  of  30  ps  pulses  and  theorized  that  damage  was  due  in  large  part  to  acoustic  (elastic)  and  shock  (inelastic) 
waves  generated  by  a  very  sharp  temperature  gradient  between  the  absorbing  melanin  granule  and  its  surroundings.  They  felt 
that  these  shock  waves  caused  mechaitical  damage  to  membranes,  and  could  thus  lead  to  cell  death.  Roider  et  al.^  theorized 
that  a  single  pulse  in  the  femtosecond  range  has  such  short  time  for  energy  deposition  (i.e.  the  pulse  duration  is  well  below 
the  thermal  relaxation  time)  that  there  is  a  small  window  of  opportunity  for  thermal  spread.  Consequently,  the  site  of  the 
lesion  remains  confined.  Cain  et  al."  reported  that  self-focusing  and  laser  induced  breakdown  (LIB)  were  piossible  factors  in  the 
damage  merhanism  by  concluding  that  ultrashort  laser  pulses  would  reach  a  "flash"  endpoint  at  which  LIB  would  occur  and  a 
plasma  would  be  formed.  With  expansion  of  the  plasma,  a  supersonic  shock  wave  and  cavitation  bubble  would  be  generated 
and  could  lead  to  damage  of  the  surrounding  retina.  They  also  stated  that  during  such  an  event,  this  cavitation  bubble  could 
course  and  re-expand  a  number  of  times.  Using  strobe  illumination  and  a  high  speed  camera,  Kelly  and  Lin^  were  able  to 
capture  images  of  this  transient  bubble  formation  inside  RPE  cells.  Below  LIB  threshold  however,  it  is  thought  that  the 
mechanism  of  damage  results  from  thermal  deposition  of  energy  into  melanosomes  with  subsequent  steam  bubble  formation 
within  the  RPE  cell.  Gerstman  et  al.®  proposed  that  the  temperature  rise  at  the  melanosome  causes  vaporization  of  the 
immediate  surrounding  medium,  thus  creating  a  bubble  that  expands  outward.  In  a  light  microscopy  study,  Toth  et  al.’ 
further  explored  these  issues  by  examining  ultrashort  laser  lesions  with  the  light  nucroscope.  They  found  that  with  the  use  of 
femtosecond  pulses,  a  narrow  column  of  damage  was  formed  within  the  retina  at  lower  energies  than  with  picosecond  pulses. 
This  supported  the  concept  proposed  by  Cain  et  al."*  that  LIB  occured  within  a  narrow  site  of  plasma  formation  when 
ultrashort  pulses  are  utilized. 
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3.  PURPOSE 


This  report  is  focused  upon  establishing  the  relationships  between  laser  energy/pulsewidth  and  the  electron 
microscopic  appearance  of  the  lesion  at  the  level  of  retinal  pigment  epithelial  (RPE)  cells,  Bruch’s  membrane,  and  choroidal 
vessels;  paying  particular  attention  to  the  effects  on  melanosomes  and  the  RPE/Bruch’s  membrane  interface. 

4.  METHODS 

The  treatment  and  procedures  used  in  this  study  conformed  to  the  Use  of  Animals  in  Ophthalmic  and  Vision 
Research  and  Federal  Guidelines.  Animals  involved  in  this  study  were  procured,  maintained,  and  used  in  accordance  with  the 
recommendations  of  the  NIH  “Guide  for  the  Care  and  Use  of  Labcoatory  Animals”  (National  Institute  of  Health  Publication 
No.  85-23,  revised  1985)  the  Office  of  Protection  from  Research  Risk  “Public  Health  Service  Pohcy  on  Humane  Care  and 
Use  of  Laboratory  Animals”  (revised  1986),  and  the  Animal  Welfare  Act. 

Two  mature  Macaca  mulatta  weighing  3-4  kg  were  maintained  under  standard  laboratory  conditions  with  monitoring 
and  care  during  laser  lesion  placement,  enucleation,  and  euthanasia  as  previously  described.  (Cain  et  al.®)  Single  laser  pulses 
with  a  range  of  energies  ftom  0.58-9.50  uJ  were  dehvered  in  a  3X3  grid  to  the  macular  area  of  each  of  three  eyes  at  one  of  the 
three  laser  pulsewidths:  90  fs,  3ps,  or  60  ps  pulses  as  previously  described  by  Cain  et  al.®  All  laser  energies  reported  in  this 
study  were  the  energies  dehvered  to  the  cornea  and  measured  as  a  percentage  of  beam  dehvered  to  a  beamsplitter.  Of  the  27 
retinal  laser  lesions  contained  within  the  three  macular  grids,  twelve  sites  were  processed  for  TEM  evaluation.  Of  these 
twelve,  four  sites  were  produced  by  60ps  laser  pulses,  five  sites  from  3ps  pulses,  and  three  sites  from  90  fs  pulses. 

The  eyes  were  incised  anterior  to  the  equator  and  immersed  in  a  3%  gluteraldehyde  and  O.IM  sodium  cocodylate 
buffer  immediately  foUowing  enucleation.  The  posterior  eye  cup  was  cut  away  from  the  anterior  segment  after  10  minutes  and 
replaced  in  the  fixative.  The  macular  area  was  later  dissected  and  embedded  in  Spurr’s  resin.  Sections  of  approximately  60  nm 
thickness  were  obtained  from  the  region  of  laser  delivery.  Each  section  was  stained  with  2%  uranyl  acetate  and  Sato  Lead 
Stain. 

Each  laser  lesion  was  examined  using  an  electron  microscope  (IEM-1200EX,  JEOL,  Peabody,  MA)  at  800  -  6000X 
magnification.  Photographs  of  these  high  magnification  images  were  taken  with  the  microscope’s  camera  and  arranged  in  a 
montage  for  each  lesion  site. 

Upon  examinahon,  particular  notice  was  paid  to  the  general  appearance  of  the  RPE  cells,  the  occurrence  of  fiactured 
or  striated  melanosomes,  the  integrity  of  Bruch’s  membrane,  and  whether  or  not  hemorrhage  or  coagulation  had  occurred 
within  the  choriocapUlaiis.  A  fractured  melanosome  was  defined  as  any  melanosome  with  clear  evidence  of  rupture,  (see 
figure,  black  arrow)  The  number  of  fractured  melansomes  per  lesion  site  was  counted  and  expressed  as  a  ratio  over  the  total 
number  of  melanosomes  per  lesion  site.  A  melansome  was  designated  as  striated  if  it  had  a  striped  appearance  and  lacked  the 
homogenous  electron  density  that  normally  appears  on  TEM.  (see  figure,  white  arrow)  Melanosomes  without  obvious 
fiacture  but  alterations  in  their  normally  uniform  contour  were  also  noted  and  deemed  irregular. 

5.  RESULTS 

Severe  damage  to  retinal  pigment  epithelial  cells  was  observed  in  9  of  12  lesions.  At  1  hour  after  laser  delivery, 
these  included  the  disappearance  of  apical  microvilli  (7/12),  condensation  of  nuclear  material  (9/12),  “smudged”  appearance  of 
cytoplasm  (likely  resulting  from  denaturation  of  cellular  protein  material)  (9/12),  and  the  formation  of  vacuoles  in  the  basal 
region  of  the  RPE  cells(9/12).  (see  figure)  At  24  hours  after  laser  delivery,  these  same  findings  were  observed  along  with  an 
apparent  “stacking”  of  RPE  cells.  The  swelling  and  vacuolization  of  mitochondria  at  the  base  of  all  RPE  cells  may  not  have 
been  the  direct  result  of  laser  exposure  but  rather  was  likely  due  to  expected  cell  degeneration  during  time  of  experimentation 
prior  to  fixation. 

RPE  cells  show  a  characteristic  “lifting  off’  of  Bruch’s  membrane  in  all  of  the  90fs  lesions.  This  separation  is 
characterized  by  the  RPE  cell  becoming  detached  with  the  basal  lamina  of  the  RPE  and  zonula  occludens  remaining  intact.  It 
appears  that  it  is  the  basal  infolding  that  is  disrupted  because  remnants  of  the  basal  infolding  are  visible  atop  the  intact  basal 
lamina.  The  3ps  and  60ps  lesions  did  not  show  this  same  separation  but  did  demonstrate  mild  elevation  of  the  RPE  cell. 
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perhaps  from  due  to  basal  vacuole  formation.  With  3ps  lesions,  only  the  highest  energy  lesion  was  clearly  elevated  (1/5). 
Conversely,  with  60ps  lesions,  all  but  the  highest  energy  lesion  were  elevated  above  Bruch’s  membrane  (3/4). 

Fractured  melanosomes  were  found  in  9  of  the  12  lesions.  It  should  be  noted  that  these  were  not  the  same 
aforementioned  9  lesions  in  which  RPE  ceU  damage  was  observed.  Fractured  melanosomes  appeared  in  the  apical  region  more 
frequently  while  none  were  found  basally.  Striated  melanosomes,  which  were  visible  in  9  of  12  lesions,  were  typically  found 
centrally,  with  a  few  seen  apical  and  basal.  Irregular  melanosomes  were  seen  in  7  of  12  lesions. 

No  apparent  relationship  was  found  between  melansome  fracture  and  peak  power.  The  highest  percentage  of  fractured 
Tnpian«!nmp.s  (36.4%)  was  Seen  in  a  lesion  created  with  the  lowest  peak  power  (3.77x10“  W),  while  no  fractures  were  seen  in 
select  lesions  produced  with  higher  peak  powers  (1.93x10^  -  1.05x10^  W).  The  lesion  produced  with  the  highest  peak  power 
(2.71x10^  W)  contained  a  midrange  percentage  of  fractures  (6.3%). 

A  trend  did  exist  between  melanosome  fracture  ratio  and  decreasing  laser  pulsewidth.  As  the  pulsewidth  decreased, 
the  percentage  of  melanosomes  showing  fracture  also  decreased.  The  highest  percentage  of  fractures  were  found  in  60ps 
lesions  (12.0%),  then  in  3ps  lesions  (6.2%),  and  the  least  in  90fs  lesions  (4.6%).  A  similar  trend  was  seen  with  the  mean 
percentages  of  striated  melanosomes. 

There  was  no  evidence  of  disruption  of  Bruch’s  membrane  in  any  of  the  12  lesions.  Furthermore,  no  hemorrhage 
from  the  choroidal  circulation  was  seen.  The  clumping  of  erythrocytes  within  choroidal  capillaries  seen  in  segments  of  3ps 
lesions  was  thought  to  be  artifact  introduced  during  specimen  handling. 

6.  DISCUSSION 

Fractured  melanosomes  were  seen  in  9  of  12  lesions,  but  as  noted,  not  the  same  9  lesions  in  which  severe  RPE  cell 
injury  was  found.  Glickman  et  al.’  theorized  that  melanosome  fracture  resulted  in  a  release  of  melanin  radicals  which  could 
cause  oxidative  damage  to  cell  proteins  and  act  as  a  prerequisite  event  to  cell  death.  The  appearance  of  cell  injury  without 
corresponding  melanosome  fracture,  as  well  as  the  observation  of  a  fractured  melanosome  within  an  RPE  cell  that  showed 
almost  no  signs  of  damage,  dispute  this  hypothesis  that  fracture  must  occur.  Nevertheless,  our  results  support  the 
relationship  between  melanosome  fracture  and  the  extent  of  damage  within  the  RPE  ceU  itself.  In  other  words,  the  percentage 
of  fractured  melanosomes  paraUels  the  severity  of  the  lesion  and  laser  pulses  that  are  likely  to  produce  a  high  percentage  of 
fractured  melanosomes  also  induce  severe  damage  to  the  rest  of  the  RPE  cell. 

Striated  melanomes  were  found  within  9  of  12  lesions.  Breathnach  et  al."  previously  described  striated  melanosomes 
as  those  still  undergoing  progressive  melanization;  whereas  mature  melanin  granules  were  "practically  uniformly  electron 
opaque".  In  1975,  Goldman  et  al.^  described  melanosomes  with  similar  features  resulting  from  ultrashort  laser  damage.  They 
theorized  that  “the  energy  not  involved  in  the  temperature  rise  at  the  surface  of  the  melanin  would  be  available  to  break  the 
cross  linkages  between  the  fibers  upon  which  the  melanin  was  deposited  in  the  developing  melanosome.  The  breakage  of 
these  bonds  would  separate  and  expose  the  fibers,  and  produce  a  striated  appearance...”.  In  our  study,  a  few  striated 
melanosomes  were  observed  in  the  control  rrticrogr£q)h  of  normal  RPE,  taken  a  distance  away  from  laser  delivery.  One 
possible  explanation  for  this  observation  is  a  species  variation  between  humans  and  other  primates.  However,  no  electron 
microgr^hs  from  previous  publications  regarding  monkey  RPE  were  found  to  show  striated  melanosomes  such  as  these.  It 
is  also  possible  that  the  tissue  removed  as  “normal  control”  was  not  adequately  distant  from  the  sites  of  laser  lesions  and 
therefore  cannot  be  thought  of  as  a  pure  control. 

Irregular  melanosomes  with  apparent  shape  alterations  were  seen  in  7  of  12  lesions.  Lacking  any  systematic  pattern 
of  ^pearance,  melanosomes  with  irregular  contour  alterations  do  not  seem  to  represent  a  specific  type  of  ultrastractural 
melanosome.  Rather,  they  are  merely  the  consequence  of  the  level  of  resolution  and  the  two  dimensional  restrictions  of 
histological  sectioning  and  imaging  and  are  therefore  of  httle  significance.  It  is  also  possible  that  these  irregularly  shaped 
melanosomes  represent  artifact  that  is  introduced  in  the  processing  of  the  tissue  specimens. 

Finding  no  breaks  in  Bruch’s  membrane  was  expected  since  the  previous  light  microscope  study  (Toth  et  al.^)  had 
not  seen  any  bleeding  from  the  choroidal  circulation.  Several  previous  studies  of  laser  effects  found  extensive  retinal  and  RPE 
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damage  and  yet  observed  no  damage  to  Bruch’s  membrane.  (Bimgruber  et  al.’;  Courant  et  al.'®;  Roider  et  al.";  Wallow  et  al.'^; 
Marshall  et  al.‘^) 

The  appearance  of  a  characteristic  detachment  of  RPE  cells  above  Bruch's  membrane  in  90fs  lesions  supports  the 
concept  of  a  cavitation  bubble  forming  and  receding  immediately  after  laser  irradiation  at  energies  above  LIB  threshold.  Such 
phenomena  could  generate  physical  forces  that  would  lead  to  this  "lifting  off  of  RPE  cells.  Goldbaum  et  al.’"  reported  that 
adhesion  between  basal  lamina  and  basal  infolding  was  stronger  than  the  integrity  of  the  RPE  cell  itself.  Thus,  in  order  to 
produce  the  observed  detachment,  some  force  must  have  acted  within  or  above  the  RPE  cells.  This  is  again  consistent  with 
the  concept  of  a  rapidly  forming  and  degenerating  cavitation  bubble  forming  within  the  RPE  cell  upon  laser  irradiation. 
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ABSTRACT 

Several  laser  systems  are  currently  under  investigation  for  the  purpose  of  removing  hard  dental  tissues. 
However,  either  undesired  thermal  side  effects  or  the  lack  of  efficiency  have  already  been  demonstrated 
in  most  cases.  In  this  paper,  advantages  and  limitations  of  using  ultrashort  laser  pulses  with  either 
picosecond  or  femtosecond  durations  are  discussed.  The  major  advantages  associated  with  these  pulse 
durations  is  the  ability  to  produce  very  precise  cavities  without  significant  thermal  side  effects.  Even 
disruptive  effects  due  to  shock  wave  generation  seem  to  be  negligible  at  moderate  pulse  energies  close  to 
the  ablation  threshold.  The  quality  of  these  cavities  is  found  to  be  superior  to  the  quality  achievable 
with  other  laser  systems.  Moreover,  a  spectroscopical  analysis  of  the  laser-induced  plasma  sparks 
enables  an  on-line  health  diagnosis  of  the  irradiated  volume.  Limitations  arise  from  the  development 
of  a  suitable  delivery  system  and  from  the  cost  of  generating  ultrashort  laser  pulses. 


1.  INTRODUCTION 

Since  the  early  stage  of  laser  development  several  attempts  have  been  made  to  apply  this  tool  for 
certain  procedures  in  dentistry,  as  well.  However,  no  satisfying  solution  has  yet  been  provided  for  one 
of  its  most  demanding  applications:  the  usage  of  lasers  in  the  therapy  of  dental  decay.  Replacing  or 
supporting  conventional  drilling  machines  by  lasers  in  order  to  achieve  a  more  accurate  and  painless 
treatment  of  dental  decay  has  been  a  challenge  for  the  past  decades.  Since  pain  is  usually  caused 
by  vibration  and  heat  due  to  friction,  it  can  be  avoided  by  using  a  contactless  operating  laser  that  is 
working  in  the  “cold”  ablation  range.  Usually  continuous  wave  (cw)  lasers  and  pulsed  lasers  with  pulse 
durations  down  to  the  microsecond  range  generate  a  lot  of  heat  in  the  region  of  the  pulp  chamber 
during  the  ablation  process.  This  is  due  to  the  fact  that  in  this  time  frame  of  pulse  duration  heat 
diffusion  plays  a  very  important  role  in  the  interaction  mechanism.^  The  first  lasers  used  in  dentistry 
are  representatives  of  this  group:  these  are  the  ruby  lasers  introduced  by  Stern  and  Sognnaes^,  and 
Goldman  et  al.^,  as  well  as  the  CO2  lasers  investigated  by  Stern  et  al.^  and  Frame.^  They  show  very 
strong  thermal  side  effects,  since  the  temperature  increase  in  the  inner  tooth  reaches  values  higher 
than  10°  C,  thereby  injuring  the  pulp. 
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Meanwhile  several  experiments  have  been  performed  using  alternative  laser  systems.  Among  these  are 
the  Er:YAG  laser  and  excimer  lasers,  especially  the  ArF  excimer  laser.  The  Er:YAG  laser  was  first 
used  in  dentistry  by  Hibst  and  Keller.®  The  wavelength  of  the  Er:YAG  laser  at  2.94  fxm  matches  the 
resonance  frequency  of  the  vibrational  oscillations  of  water  molecules  contained  in  the  teeth,  thereby 
strongly  enhancing  the  absorption  of  Er:YAG  radiation.  The  absorption  leads  to  small  microexplosions 
breaking  the  hydroxyapatite  structure.  The  coincidence  of  thermal  and  mechanical  ablation  effects  has 
led  to  the  term  “thermomechanical”  interaction.'^  However,  as  reported  by  Niemz  et  al.®  and  Frentzen 
et  ah,®  cracks  up  to  300  /rm  deep  may  be  induced  by  this  type  of  interaction  that  can  potentially  be 
the  origin  for  new  caries  development.  Liesenhoff  et  al.^®  found  that  the  ArF  laser  at  a  wavelength 
of  193  nm  shows  only  very  little  thermal  side  eflPects.  Heat  induced  ruptures  are  reduced  due  to  the 
shorter  pulse  duration  of  typically  15  ns.  However,  the  ablation  rate,  i.e.  the  ablated  volume  per  pulse, 
is  very  low  according  to  Rechmann  et  al.^^  This  ineffectiveness  and  the  general  risks  of  UV  radiation 
are  the  major  disadvantages  concerning  its  use  in  dentistry,  although  it  may  be  of  importance  for 
other  medical  applications. 

A  new  approach  was  made  by  Niemz  et  al.®  when  first  applying  picosecond  laser  pulses  to  dental 
applications.  In  this  so-called  plasma-induced  ablation, a  localized  microplasma  is  induced  at  the 
laser  focus.  With  pulse  durations  of  a  few  picoseconds,  pulse  energies  up  to  1  mJ  and  a  focal  spot  size 
of  about  30  /rm  the  power  density  at  the  focus  reaches  values  >  10^^  W/cm^.  Hence,  the  amplitude  of 
the  electric  field  exceeds  10'  V/cm  and  optical  breakdown  occurs.  The  ionized  material  is  ablated,  and 
thermal  as  well  as  mechanical  damage  is  negligible  when  using  pulse  energies  close  to  the  threshold  of 
plasma  ignition.  Recently,  a  further  step  forward  was  made  when  using  femtosecond  lasers  to  ablate 
hard  dental  tissues. In  general,  the  results  of  these  studies  show  great  promise  that  the  initial 
difficulties  of  “laser  dentistry”  can  finally  be  overcome  when  choosing  the  correct  laser  parameters, 
since  femtosecond  pulses  as  well  as  picosecond  pulses  produce  ablations  superior  in  quality  to  those 
achievable  by  longer  pulse  durations.  The  question  remains  whether  these  highly  sophisticated  ultra- 
short  pulse  laser  systems  -  usually  consisting  of  an  oscillator  laser  and  a  regenerative  amplifier  -  will 
be  able  to  form  a  marketable  alternative  to  conventional  drilling  machines. 


2.  MATERIALS  AND  METHODS 

The  picosecond  laser  system. 

The  Nd:YLF  laser  system  is  designed  as  a  two  stage  combination  (Fig.  1),  consisting  of  a  laser  oscillator 
and  a  regenerative  amplifier,  to  provide  laser  pulses  with  durations  down  to  30  ps  and  energies  up 
to  1  mJ  at  a  wavelength  of  1053  nm.  The  oscillator  NdrYLF  rod  is  pumped  by  a  temperature  tuned 
1  Watt  diode  laser  (DL)  using  beam  shaping  and  collimating  optics  (CO).  The  Nd:YLF  crystal  itself  is 
coated  with  a  high  reflecting  mirror  on  the  surface  pointing  to  the  diode  laser.  An  acousto-optic  mode 
locker  (AOM)  is  placed  near  the  flat  10  %  output  coupler  (OC).  For  active  amplitude  modulation  an 
amplified  80  MHz  signal  is  applied  to  this  device,  generating  a  train  of  short  laser  pulses  with  typical 
durations  of  25  ps  each.  A  real-time  autocorrelation  system  allows  continuous  supervision  of  the  pulse 
width.  For  the  purpose  of  selecting  the  1053  nm  transition,  a  Brewster  polarizer  (BP)  is  added  to  the 
cavity.  At  the  half-wave  plate  (HWP)  the  160  MHz  pulse  train,  consisting  of  0.2  nJ  pulses,  experiences 
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a  90®  rotation  of  the  polarization  vector.  Using  a  4  %  reflecting  mirror  (M3)  and  a  polarizing  beam 
splitter,  the  oscillator  pulses  are  then  injected  into  the  regenerative  amplifier  unit  described  by  Bado 
et  al.^^  The  76  mm  amplifier  Nd:YLF  rod  is  pumped  by  a  single  flashlamp,  controlled  by  a  standard 
laser  power  supply  (model  204A,  Quantronix  Inc.).  The  cavity  employs  two  highly  reflecting  mirrors 
with  a  radius  of  curvature  of  1  m  each.  Applying  a  2  kV  voltage  signal  with  up  to  1  kHz  repetition  rate 
to  a  LiNbOs  Pockels  cell  (PC)  provides  half-wave  retardation  per  round-trip.  In  combination  with  the 
double-pass  half-wave  retardation  of  the  intracavity  quarter- wave  plate  (QWP),  a  selected  oscillator 
pulse  is  seeded  and  trapped  in  the  amplifier  unit.  The  driving  of  the  Pockels  cell  is  synchronized  to 
the  mode  locking  process  by  feeding  the  80  MHz  signal  into  a  special  divider  and  timer  logic.  After 
about  100  roundtrips  in  the  cavity  the  seeded  pulse  reaches  its  saturation  limit.  The  pulse  energy 
can  be  boosted  up  to  1  mJ,  corresponding  to  an  amplification  of  10®  of  the  oscillator  output  energy. 
Installation  of  an  aperture  (A)  restricts  the  laser  operation  to  the  fundamental  TEMqo  mode.  At 
maximum  gain  the  Pockels  cell  driver  switches  back  to  0  V,  causing  no  retardation.  The  polarization 
vector  is  now  rotated  by  90®  as  the  pulse  double-passes  the  quarter-wave  plate  and  the  Pockels  cell  in 
the  left  part  of  the  cavity.  Consequently,  the  amplified  pulse  is  then  reflected  at  the  polarizing  beam 
splitter  and  dumped  out  of  the  regenerative  amplifier.  Mirror  M3  is  now  transmitting  96  %  of  the 
amplified  pulse  energy.  Using  mirror  M4,  the  pulse  train  is  finally  injected  into  the  application  unit. 
Autocorrelation  of  these  pulses  shows  that  their  pulse  duration  has  slightly  increased  to  about  30  ps 
due  to  dispersion  inside  the  amplifier  cavity. 


Fig.  1:  The  picosecond  Nd:YLF  laser. 
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For  the  sake  of  automating  these  experiments  an  application  unit  was  developed.  This  device  consists 
basically  of  delivering  optics  and  a  computer  controlled  three-axes  translation  stage.  After  expanding 
the  Nd:YLF  laser  beam  four  times  by  the  lenses  L2  and  L3,  it  is  tightly  focused  onto  the  tissue  sample 
by  the  lens  L4.  This  focusing  lens  is  made  of  Si02  and  has  a  focal  distance  of  100  mm.  The  focus 
spot  is  measured  with  the  knife  edge  method  and  has  a  diameter  of  about  30  //m.  Stepping  motors 
connected  to  the  translation  stage  allow  precise  spot-to-spot  movements  of  the  tissue  within  1  ^m.  A 
software  package  gives  the  user  a  choice  of  different  ablation  patterns.  Primarily,  square  geometries 
have  been  selected. 


The  femtosecond  laser  system. 

In  order  to  investigate  the  effects  of  femtosecond  laser  pulses  on  hard  dental  tissues,  a  Ti:Sapphire 
laser  source  (model  Spitfire,  Spectra-Physics)  was  used.  This  laser  is  able  to  emit  laser  pulses  with 
durations  of  130  fs  at  a  central  wavelength  of  780  nm.  For  these  studies,  the  pulse  energy  is  adjusted 
to  50  /iJ  at  a  repetition  rate  of  1  kHz.  The  laser  beam  is  focused  to  a  spot  of  30  ^m  in  diameter.  The 
studies  on  human  enamel  with  this  laser  system  were  performed  at  the  Laserzentrum  Hannover  by 
the  order  of  Dr.  med.  A.  Kasenbacher  (Traunstein,  Germany). 


Diagnostic  system. 

For  the  purpose  of  analyzing  dental  decay,  a  spectroscopical  arrangement  was  set  up  as  shown  in 
Fig.  2.  The  laser- induced  plasma  spark  is  optically  imaged  onto  the  entrance  pupil  of  a  spectrometer. 
The  readout  of  the  photomultiplier  tube  (PMT)  is  triggered  by  the  laser  pulses.  The  spectra  are 
recorded  between  400  nm  and  600  nm.  A  BBO  crystal  is  used  to  normalize  the  relative  intensities  at 
the  wavelength  of  the  second  harmonic  of  the  laser. 


BBO 


PC  Oscilloscope 


Spectrometer 


Fig.  2:  Spectroscopy  of  laser-induced  plasma  sparks. 
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Tooth  preparation. 

The  teeth  were  obtained  from  the  Dental  School  of  the  University  of  Heidelberg  with  a  protocol 
approved  by  an  institutional  review  board.  The  molars  were  always  kept  in  humid  environment  to  avoid 
cracking  due  to  dryness.  Artificial  caries  was  prepared  according  to  Silverstone^^  by  using  synthetic 
saliva  containing  1  %  hydroxymethylcellulosis  buffered  at  pH  4.8  with  acetic  acid.  After  the  ablation 
experiments  the  teeth  were  dried  in  the  vacuum  chamber  of  an  exsiccator  and  coated  with  a  30  nm 
thick  gold  layer.  The  morphologic  changes  induced  by  the  laser  radiation  were  studied  using  a  scanning 
electron  microscope  (model  1810,  Amray  Inc.). 

3.  RESULTS 

Figures  3a-b  show  scanning  electron  micrographs  (SEM)  of  cavities  generated  by  the  Nd:YLF  picose¬ 
cond  laser.  Fig.  3a  represents  a  cavity  ablated  in  sound  enamel,  whereas  Fig.  3b  shows  a  cavity  within 
artificial  caries.  Both  cavities  have  lateral  dimensions  of  1  x  1  mm^  and  a  depth  of  about  400  fj,m.  They 
were  created  by  distributing  1  mJ  laser  pulses  onto  40  lines  over  the  tooth  surface  with  400  lasered 
spots  per  line,  and  repeating  this  procedure  ten  times  for  the  cavity  in  Fig.  3a  and  once  for  the  cavitj 
in  Fig.  3b.  Thus,  a  total  number  of  160  000  laser  shots  was  necessary  for  the  cavity  in  sound  enamel 
and  only  one  tenth  of  this  number  was  needed  to  generate  the  cavity  in  carious  enamel.  Hence,  at  a 
pulse  energy  of  1  mJ  and  a  repetition  rate  of  1  kHz,  a  cavity  as  shown  in  Fig.  3a  is  ablated  within 
a  total  time  of  160  seconds.  Therefore,  the  efiiciency  for  removing  sound  enamel  with  the  described 
laser  system  is  about  0.15  mm^  per  minute  and  millijoule.  By  contrast,  the  cavity  shown  in  Fig.  3b 
requires  a  duration  of  only  16  seconds  and  is  characterized  by  an  efiiciency  of  1.5  mm^  per  minute  and 
millijoule,  respectively.  Both  cavities  have  a  very  precise  geometry,  and  especially  the  edges  are  very 
clean  and  sharp.  The  bottom  surfaces  of  cavities  within  sound  enamel  are  slightly  rougher  than  the 
surfaces  in  carious  enamel.  The  overall  roughness  is  of  the  order  of  10  ^m  and  ideally  suited  for  most 
filling  materials. 


Fig.  3:  Sound  enamel  (a)  and.  carious  enamel  (b)  ablated  with  the  Nd:YLF  picosecond  laser. 


intensity 


By  contrast,  Fig.  4  shows  a  cavity  in  a  human  tooth  produced  by  the  Ti:Sapphire  femtosecond  laser. 
The  cavity  has  lateral  dimensions  of  approximately  1x1  mm^  and  a  depth  of  about  315  /xm.  It  was 
created  by  distributing  12  scans  over  the  tooth  surface  with  a  total  of  approximately  two  million  pulses 
at  an  energy  of  50  /xJ  each.  Therefore,  the  efficiency  for  removing  sound  enamel  with  this  femtosecond 
laser  system  is  about  0.18  mm^  per  minute  and  millijoule. 


Fig.  Sound  enamel  ablated  with  the  Ti:Sapphire  femtosecond  laser. 

The  results  of  plasma  spectroscopy  are  summarized  in  Figs.  5a-b.  Typical  spectra  obtained  from 
plasma  sparks  on  either  sound  enamel  (Fig.  5a)  or  carious  enamel  (Fig.  5b)  are  shown.  Due  to  the 
demineralization  associated  with  dental  decay,  a  strong  decrease  in  intensity  of  the  calcium  transitions 
becomes  apparent. 


Fig.  5:  Plasma  spectrum  of  sound  enamel  (a)  and  carious  enamel  (b). 
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4.  DISCUSSION 


The  ablation  experiments  performed  with  both  the  NdrYLF  picosecond  laser  and  the  TirSapphire  fem¬ 
tosecond  laser  prove  that  cavity  preparation  is  possible  with  ultrashort  laser  pulses  without  inducing 
severe  thermal  or  mechanical  damage.  This  is  a  little  surprising,  because  deep  cracks  inside  the  tooth 
due  to  shock  wave  generation  were  predicted  by  Frentzen  and  Koort,^  when  using  ultrashort  laser 
pulses.  Such  shock  wave  effects  were  observed  indeed,  when  focusing  the  Nd:YLF  laser  onto  90  fiia. 
thick  microscope  glass  plates  as  reported  by  Pioch  et  However,  no  tooth  alterations  due  to  shock 
waves  -  like  cracking  or  structural  changes  -  have  been  observed  so  far  when  applying  moderate  pulse 
energies  close  to  the  ablation  threshold.  These  observations  can  probably  be  attributed  to  the  different 
morphologic  structures  of  glass  and  teeth.  Acoustic  waves  are  strongly  damped  by  the  organic  matrix 
of  teeth,  while  glass  easily  splinters  when  being  exposed  to  high  pressure  waves. 

The  reproducibility  of  clean  and  regularly  shaped  cavity  preparation  in  teeth  with  both  picosecond 
and  femtosecond  lasers  is  very  high,  since  the  cavities  obtained  in  all  teeth  investigated  were  of  a 
similar  high  quality.  Hence,  strong  experimental  evidence  is  given  that  these  lasers  are  suitable  for 
clinical  applications.  The  observed  roughness  of  the  cavity  bottoms  should  enable  a  proper  adhesion 
of  filling  materials. 

When  discussing  the  limitations  of  a  new  technique,  one  should  first  look  at  potential  delivery  systems. 
For  a  clinical  application  of  the  lasers  discussed  above,  optical  fiber  delivery  and/or  an  even  faster 
scanning  mirror  device  need  to  be  developed.  When  aiming  at  treatment  durations  comparable  to  the 
conventional  drill,  energy  outputs  of  a  few  hundred  microjoules  at  repetition  rates  of  up  to  10  kHz 
should  be  sufficient  for  clinical  applications.  Both  Nd:YLF  and  Ti:Sapphire  lasers  with  repetition  rates 
of  5  kHz  and  higher  are  already  commercially  available.  Another  limitation  is  the  cost  of  such  highly 
sophisticated  laser  systems.  Laser  therapy  alone  might  not  justify  more  expensive  treatment  costs. 
However,  the  present  paper  indicates  that  the  mechanism  of  plasma-induced  ablation  enables  selective 
removal  of  dental  decay,  as  well.  On  the  one  hand,  the  ablation  rate  of  dental  decay  is  higher  than  that 
of  sound  enamel.  On  the  other  hand,  the  generated  plasma  spark  can  be  spectroscopically  analyzed  in 
real-time.  With  this  additional  option  a  computer  controlled  picosecond  or  femtosecond  laser  system 
could  possibly  decide  in  the  near  future,  whether  all  dental  decay  has  been  removed  or  not,  thereby 
minimizing  the  amount  of  hard  dental  tissue  to  be  ablated.  Since  the  cavity  size  is  then  reduced  -  even 
after  a  refilling  -,  the  number  of  indications  for  a  dental  crown  or  prothesis  might  decline.  In  these 
cases,  i.e.  when  the  advantages  of  minimally  invasive  surgery  become  significant,  the  problem  of  cost 
might  be  overcome. 
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ABSTRACT 

Plasma  luminescence  spectroscopy  was  used  for  precise  ablation  of  bone  tissue  without  damaging 
nearby  soft  tissue  using  an  ultrashort  pulse  laser  (USPL).  Strong  contrast  of  the  luminescence 
spectra  between  bone  marrow  and  spinal  cord  provided  the  real  time  feedback  control  so  bone 
tissue  is  selectively  ablated  while  preserving  the  spinal  cord. 

Keywords  :  ultrashort  pulse  laser,  plasma,  luminescence  spectroscopy,  feedback  control 

1.  INTRODUCTION 

As  lasers  are  accepted  as  possible  replacements  for  conventional  surgical  tools  in  many  medical 
fields,  a  focus  of  research  has  been  on  more  effectively  removing  the  biological  tissues  while 
reducing  mechanical  and  thermal  damage  to  the  surrounding  tissues.  Recently,  many  studies  have 
shown  that  ultrashort  pulse  lasers  (<  1  ps)  can  meet  this  need’’®.  In  this  short  pulse  domain,  the 
ablation  process  is  distinguished  from  the  conventional  longer  pulse  ablation  which  mostly 
depends  on  the  photothermal  heating  in  the  media.  As  the  pulse  width  becomes  shorter, 
multiphoton  ionization  dominates  during  the  ablation  process  which  leads  to  the  plasma  optical 
breakdown.  The  plasma  is  ejected  at  a  rate  ~10’  cm/s  which  is  faster  than  the  shock  wave  speed 
~10®  cm/s  and  the  even  slower  heat  diffusion.  Therefore,  the  mechanical  and  thermal  damage  is 
minimized  with  these  short  pulses.  Additionally,  since  the  plasma  is  reflective  and  absorptive,  for 
longer  pulses  such  as  1  ns  pulses,  most  of  the  incoming  energy  is  reflected  back  or  channeled  into 
plasma  heating  instead  of  being  absorbed  by  the  tissue  and  used  for  tissue  removal.  On  the  other 
hand,  the  energy  deposition  time  for  ultrashort  pulses  is  much  shorter  than  the  plasma  expansion 
time  and  therefore  the  ablation  is  highly  efficient  (~l)xm/pulse)^. 

Luminescence  is  generated  from  a  high  temperature  plasma  during  ablation.  Strong 
luminescence  signals  for  calcium  based  plaque  in  atherosclerotic  tissue  has  been  reported’  ®.  Our 
studies  are  focused  on  a  feasibility  test  of  using  the  luminescence  spectroscopy  as  a  tool  for  optical 
feedback  during  micro-spinal  surgery  where  die  calcium  based  bone  tissue  needs  to  be  removed 
while  preserving  the  soft  nerve  tissue  which  stems  from  the  spinal  cord.  The  results  suggest  that 
the  hard  tissue  can  be  selectively  ablated  with  high  accuracy. 

2.  RESULTS  AND  DISCUSSIONS 

Our  ultrashort  pulse  laser  ablation  system  consists  of  four  separate  lasers.  An  82  MHz  Ti-Sapphire 
actively  mode-locked  laser  (Spectra  Physics,  Model  #  3960)  is  pumped  by  a  5  W,  frequency 
doubled  Nd:YAG  laser  (Spectra  Physics  Model :  Millenia)  mnning  at  532  nm.  The  mode-locked 
laser  pulse  has  duration  of  100  fs  at  790  nm.  Its  pulse  is  amphfi^  by  a  Ti-Sapphire  regenerative 
amplifier  (Positive  Light,  Model  :  Spitfire)  through  a  chirped  pulse  amplification  (CPA)  process. 
This  amplifier  is  pumped  by  a  10  W,  527  nm  Nd:YLF  laser  (Positive  Light,  Model :  Merlin).  The 
final  pulse  duration  is  about  150  fs  running  at  1  kHz  and  its  amplified  energy  is  more  than  1 
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mJ/pulse  at  790  nm.  The  ablation  rate  is  approximately  1  mm/s  using  a  1  kHz  beam  train.  The 
focused  beam  size  was  approximately  100  |im  (FWHM)  with  TEMqo  mode. 

A  fresh  porcine  spine  was  frozen  and  cut  transversely  to  expose  the  bone  marrow  and  spinal 
cord.  The  cross  section  was  cleaned  with  a  soft  brush  to  remove  the  debris  of  the  soft  tissue  left  in 
the  bone  marrow  during  cutting.  The  laser  pulses  were  focused  onto  the  cross  sections  of  spinal 
cord  and  bone  marrow  to  generate  the  plasma  luminescence. 

The  luminescence  signal  was  collected  by  a  1  mm  diameter  optical  fiber  and  was  delivered  to 
the  spectrometer  and  charge  coupled  device  (CCD)  camera.  The  source  light  was  normally  incident 
onto  the  tissue  and  the  fiber  was  placed  20°  fi'om  the  normal  direction  and  5  mm  away  from  the 
tissue  surface.  A  mechanical  shutter  for  the  source  laser  and  a  CCD  shutter  were  controlled  by  a 
pulse  generator.  Fig.  1  shows  the  luminescence  spectra  from  both  bone  and  spinal  cord.  As  seen 
from  the  figure,  not  only  the  overall  luminescence  intensity  is  strong  but  strong  calcium  lines  are 
observed  from  the  bone  luminescence  spectrum.  On  the  other  hand,  the  luminescence  from  the 
spinal  cord  is  much  smaller  and  there  are  no  calcium  line  features. 

The  overall  luminescence  intensity  decreases  as  the  ablation  front  moves  deeper  into  the 
ablation  hole  because  less  ablated  material  is  ejected  and  the  distance  between  the  ablation  front  and 
detector  becomes  larger  as  discussed  in  the  previous  study'®.  Fig.  2  shows  the  temporal  changes 
of  the  absolute  luminescence  intensities  inte^ated  over  each  100  shots  during  1  sec,  1  kHz  pulse 
train  As  expected,  the  total  intensity  of  luminescence  decreases  rapidly  as  the  ablation  front  moves 
deeper  into  the  sample  while  the  luminescence  intensity  does  not  vary  significantly  for  spinal  cord. 
After  the  first  600  shots,  the  overall  luminescence  from  spinal  cord  becomes  stronger  than  that 
from  the  bone  tissue.  From  this  study,  it  was  concluded  that  by  comparing  the  luminescence  ratio 
between  616  nm  and  575  nm,  we  can  accurately  discriminate  the  two  tissue  types  better  thm 
comparing  the  absolute  intensities.  The  luminescence  ratio  between  the  two  wavelengths  is 
consistent  for  both  bone  and  spinal  cord  despite  the  fact  that  the  absolute  intensity  of  the  bone 
changes  dramatically  during  1000  shots  as  shown  in  Fig.  3.  The  ratios  between  these  two 
wavelengths  (616  nni/575  nm)  were  4.8  (±  1.06)  for  bone  and  1.4  (±  0.27)  for  spinal  cord. 

In  a  more  practical  setup  as  shown  in  Fig.  4,  this  beam  was  delivered  and  focused  onto  the 
tissue  using  an  articulated  arm  which  is  composed  of  seven  separate  high  damage  threshold  mirrors 
and  one  focusing  lens.  The  luminescence  signal  was  collected  by  a  200  pm  optical  fiber  which  is 
attached  on  the  handpiece  of  the  articulated  arm  and  connected  to  a  1x2  fiber  coupler.  This  signal 
was  detected  by  two  photomultipher  tubes  (PMT)  which  are  equipped  with  616  nm  and  575  nm 
bandpass  filters  respectively.  To  remove  the  intense  laser  source  light,  an  additional  short  wave 
pass  filter  with  cut-off  wavelength  at  650  nm  were  added  to  each  PMT's  providing  optical  density 
of  7  for  the  source  light.  The  typical  luminescence  signals  for  these  two  detectors  are  shown  in 
Fig.  5(a)  for  bone  and  in  Fig.  5(b)  for  spinal  cord.  The  initial  strong  peak  corresponds  to  the  laser 
source  light  which  is  too  strong  to  be  filtered  even  by  the  OD  7  filters.  Each  detected  signal  was 
gated  and  integrated  immediately  after  the  strong  laser  signal  so  that  only  the  plasma  luminescence 
is  collected.  The  integrated  signed  was  compared  in  a  computer  and  a  TIL  signal  was  generated  to 
determine  if  the  tissue  is  bone  or  spinal  cord  and  to  control  the  laser  shutter.  When  the  laser  hit  the 
spinal  cord,  the  computer  sends  an  "off  signal  to  the  mechanical  shutter  to  close  it.  After  the 
ablation  stops,  the  shutter  reopens  in  0.5  sec  so  that  the  surgeon  can  keep  ablating  the  bone. 

Since  luminescence  spectroscopy  requires  a  small  fraction  of  the  tissue  to  be  ablated,  it  is  of 
concern  how  to  minimize  the  damage  to  tide  soft  tissue.  As  mentioned  earlier,  the  ablation  rate  for 
this  short  pulse  width  is  approximately  1  jim/pulse.  The  future  goal  of  the  study  will  focus  on 
limiting  the  damage  to  less  than  5  |xm  which  is  believed  to  cause  negligible  damage  to  nerve  tissue. 
To  accomplish  this,  a  fast  electronics  package  is  under  development.  Currently,  the  possible 
maximum  damage  is  between  10  - 15  |im. 


93 


3.  CONCLUSION 


We  demonstrated  that  a  safe  and  precise  microsurgery  system  using  ultrashort  pulse  laser  is 
feasible.  The  two  wavelength  comparison  technique  provides  accurate  selective  tissue  ablation. 
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ABSTRACT 


Three  en2yines  differing  in  their  structural  composition  were  irradiated  by  UV  lasers  to  study  the  effect  of  temperature, 
protein  concentration  and  addition  of  small  molecules  on  their  sensitivity  to  radiation  exposure.  The  laser  -  induced  effects 
were  due  to  the  structural  complexity  of  the  protein  molecules  and  depended  on  the  dose  applied,  the  wavelength  and  the 
density  of  irradiation.  The  multi-enzyme  2-oxoglutarate  dehydrogenase  complex  (20GDC)  was  subjected  to  pronounced 
irradiation  -  induced  changes  whereas  the  response  of  the  two  other  enzymes  was  less  significant.  Reduction  of  the  protein 
levels  in  irradiated  samples  was  important  under  the  XeCl  laser  coercion  and  the  effects  depended  on  the  doses  applied.  The 
laser  irradiation  effects  are  suggested  to  be  realized  by  means  of  conformational  changes  in  the  protein  molecules  and 
intermolecular  association  -  dissociation  processes. 

Keywords:  UV  laser  irradiation,  2-oxoglutarate  dehydrogenase  complex,  fluorescence. 


1.  INTRODUCTION 

Laser  irradiation  of  enzymes  causes  significant  changes  in  their  activity.  It  has  been  shown  that  a  chromophore  -  assisted 
laser  inactivates  93%  of  fl-galactosidase  activity  and  80%  of  alkaline  phosphatase  as  well  as  87%  of  acetylcholinesterase 
activity.  Although  thermal  denaturation  and  photochemical  mechanisms  for  proteins  inactivation  were  postulated',  the 
precise  nature  of  the  laser  -  mediated  damage  to  the  protein  fimction  has  not  been  established. 

Three  enzymes  differing  in  the  complexity  of  their  structural  compositions  were  studied  in  connection  with  laser-  induced 
changes  in  their  functional  activities.  Thiamine  triphosphatase  (ThTPase:  EC  3.6.1.28),  characterized  by  the  simplest 
structure  presented  by  one  polypeptide  chain  (Mr  28  000  ),  alcohol  oxidase  (AO:  EC  1.1.3.13),  built  of  four  identical 
subunits  (Mr  300  000)  having  coenzyme  FAD  (flavine  adenine  dinucleotide)  in  its  composition,  and  2  -oxoglutarate 
dehydrogenase  complex  as  an  ensemble  of  three  separate  enzymes:  2  -  oxoglutarate  dehydrogenase  (Elo:  EC  1. 2.4.2)  , 
dihydrolipoamide  succinyltransferase  (E2o:EC  2.3.1.61)  and  dihydrolipoamide  dehydrogenase  (E3  :  EC  1.8. 1.4)  , 
possessing  TPP,  lipoic  acid  residues  and  FAD  incoiporated  as  cofactors  into  Elo,  E2o  and  E3,  accordingly,  have  been 
examined  for  their  activity  under  the  effect  of  the  UV  laser  irradiation. 

The  fluorescence  of  2-oxoglutarate  dehydrogenase  complex  from  bovine  heart  after  the  irradiation  by  the  XeCl  and  the 
nitrogen  lasers  was  studied.  The  fluorescence  intensity  of  bound  flavine  adenine  dinucleotide  was  increased  after  high  doses 
of  the  excimer  laser  treatment. 


2.  MATERIALS  AND  METHODS 

The  sources  of  coherent  UV  irradiation  were  a  pulse  nitrogen  laser  and  an  excimer  XeCl  laser.  The  electron  discharging 
XeCl  laser  designed  at  the  Grodno  State  University  served  as  a  source  of  UV  irradiation  (X  =308  nm).  The  laser  performs 
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the  following  parameters:  the  generating  energy  of  the  pulse  was  around  2J  with  a  pulse  frequency  of  1  Hz  and  halfwidth 
equal  to  50  ns.  The  specially  designed  XeCl  laser  enables  us  to  vary  the  pulse  duration  from  25  to  100  ns  with  generating 
energy  of  0.6  mJ.  The  protein  samples  were  subjected  to  laser  pulses  with  energy  of  25,  50  and  150  mJ  and  60  ns  halfwidth. 
The  pulse  was  generated  by  separation  of  the  laser  beam  producing  the  energy  required  in  the  experiment.  The  beam  was 
focused  with  a  long  -  focal  quartz  lens  (F=lm)  and  directed  vertically  on  the  tubes  containing  protein  solutions  to  obtain  a 
spot  on  the  sample  with  diameter  of  5  mm. 


Standard  nitrogen  laser  (LGP  -  505)  was  used  as  source  of  the  UV  irradiation  (X,=337  nm)  to  generate  a  pulsed  laser  beam 
with  a  pulse  width  of  10ns  at  a  frequency  ranging  from  12  to  1000  Hz.  The  protein  samples  were  subjected  to  1  min,  3  min 
and  5  min  of  laser  pulses  of  0.25  mJ  per  pulse  at  a  frequency  of  16  Hz  and  a  power  4  MW  that  corresponded  to  the  doses  of 
irradiation  equal  to  240  mJ,  720  mJ  and  1200  niJ.  Use  of  different  doses  of  irradiation  by  the  two  lasers  was  caused  by  the 
facts  that  eximeric  laser  generates  twofold  higher  pulse  power  then  nitrogen  one  and  the  activity  of  the  XeCl  laser  emanation 
was  higher.  The  experiments  done  have  shown  that  the  lowest  dose  of  nitrogen  laser  irradiation  when  some  changes  of  the 
functional  activity  of  20GDC  become  noticeable  corresponded  to  the  1  min  of  the  protein  complex  treatment. 

20GDC  was  extracted  from  bovine  heart  and  purified  by  the  method  of  Stanley^  with  our  modifications^.  Enzyme  activity 
was  assayed  by  monitoring  NADN  formation  at  340nm  and  30®C,  alcohol  oxidase  was  isolated  from  Candida  boidinii  cells 
by  a  specially  developed  method'*  and  thiamine  triphosphatase  was  purified  by  a  method  developed  for  the  enzyme  from 
bovine  bram^.  The  computations  were  performed  by  using  linear  least  square  regression^.  The  kinetic  parameters  were 
calculated  using  Lineweaver-  Burk  plots’. 


Steady-state  fluorescence  spectra  of  the  protein  solutions  were  recorded  by  the  SDL-2  and  DFS-52  spectrometers  (LOMO, 
Russia).  Fluorescence  decay  measurements  were  performed  with  a  nanosecond  pulse  fluorimeter  *.  The  samples  were 
excited  using  a  nitrogen  flash  lamp  generating  exciting  pulse  with  a  half  width  of  1  ns.  A  PC  486DX  was  used  for  the 
operation  and  data  processing.  The  decay  data  analysis  (deconvolution)  was  performed  according  to  the  method  taking  into 
accoxmt  the  finite  width  of  the  profiles  of  the  excitation  pulse  and  the  response  function  of  the  apparatus.  ■ 


Experimentally  obtained  fluorescence  decay  curve  I(t)  is  considered  as  a  convolution  of  the  decay  law  of  the  sample  F(t) 
wi&  the  instrument  response  function  L(t) 

t 

I{t)=  \L(t')F{t-t')df  (1) 

0 


The  fluorescence  decay  law  is  usually  simulated  as  a  sum  of  exponentials 


(2) 


where  cxj,  and  Xj  are  the  decay  amplitude  and  lifetime  respectively. 

•  •  13 

Using  non-linear  least  squares  method  the  fitting  parameters  oq,  xj  were  determined  by  minimizmg  , 

where  the  index  i  sums  over  the  appropriate  chaimel  region.  !(/,)  and  Y(t,)  denote,  respectively,  the  observed  and  calculated 
values  corresponding  to  the  i  chaimel,  and  W(ti)  is  the  corresponding  statistical  weight,  v  represents  the  number  of  degrees 
of  freedom  for  the  fluorescence  decay  surface. 

W(tj)  is  defined  as  a  magnitude,  reciprocal  to  the  variance  of  the  observed  fluorescence  intensity  in  i  channel. 

^—  =  var /(/,;.  (4) 

W(t,) 

We  used  a  modified  expression  for  statistical  wei^ts  to  take  into  account  the  mfluence  of  noise  on  the  instrument  response 
function  L(t). 

—  =  var  I(tJ  +  Jvnr  L(t')  ■  F^(t,  -  f)dt' .  (5) 

^\h/  0 

This  modification  is  useful  when  the  average  lifetime  is  less  then  or  equal  to  the  excitation  pulse  duration. 
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To  monitor  the  dependence  of  the  apparatus  response  function  on  the  wavelength  a  “shift  method”''  was  used. 

The  protein  samples  of  OGDC,  AO  and  ThTPase  were  prepared  by  dissolving  of  the  protein  either  in  0.02  M  phosphate  or 
Tris  -  buffer  (pH  7.4). 


3.  RESULTS  AND  DISCUSSION 

The  radiation  studies  carried  on  purified  en2ymes  were  performed  using  an  UV  laser  equipment  generating  308  and  337  nm 
laser  light.  All  the  irradiation  procedures  were  performed  at  4®  C  (excluding  the  experiments  on  increasing  temperature 
effects).  The  possible  effects  of  the  protein  concentration  in  irradiated  sample,  temperature  and  the  addition  of  small 
molecules-  dithiothreitol  and  sodium  azide  (NaNj)  on  the  functional  activity  of  three  enzymes  of  interest- were  examined. 

Thiamine  triphosphatase  possessing  the  simplest  composition  was  unaffected  by  any  of  the  lasers  used.  The  response  of 
alcohol  oxidase,  a  tetrameric  flavoprotein,  to  excimer  laser  light  exposure  was  noticeable.  It  was  established  that  the  effects 
obtained  depended  on  the  protein  concentration  It  was  shown  that  when  the  protein  level  in  the  sample  was  equal  to  0.25 
mg/ml  the  laser  irradiation  at  a  dose  of  50  and  200  mJ  diminished  the  AO  activity  in  both  the  cases.  The  kinetic 
investigations  provide  evidence  for  the  Ymax  being  80%  of  the  initial  value  and  the  BCm  did  not  change.  In  increasing  AO 
concentration  to  0.5  mg/ml  the  effect  of  these  dose  irradiation  was  not  detected.  However,  soon  after  the  treatment  of  more 
concentrated  preparations  ( after  storing  of  the  enzyme  during  24  hours )  the  AO  activity  decreased  quickly.  The  preliminary 
introduction  of  the  alcohol  oxidase  stabiliser  -  NaN3  -  to  the  alcohol  oxidase  solution  did  not  only  protect  AO  fi'om 
inactivation  but  prolonged  the  enzyme  lifetime.  The  further  increasing  of  the  AO  concentration  to  1  mg/ml  led  to  some 
activation  of  the  enzyme  (up  to  130%). 

Being  the  ensemble  of  tree  constituent  enzymes,  the  multi  -  enzyme  2-oxoglutarate  dehydrogenase  eomplex  was  mostly 
affected  by  UV  laser  irradiation.  It  has  been  shown  that  nitrogen  laser  caused  a  significant  decrease  of  the  OGDC  activity 
(Table  1).  The  effect  did  not  depend  on  the  protein  concentration  in  the  preparations  of  the  complex.  The  irradiation  of  the 
enzyme  solution  at  240C  resulted  in  a  loss  of  12%  of  its  activity  even  in  the  case  of  short  (Imin)  treatment  by  the  beam 
possessing  medium  power  of  the  pulse  (720  mJ). 

The  effect  of  the  temperature  increase  was  unnoticed  with  using  of  XeCl  laser,  vdiereas  the  protein  level  was  an  important 
factor.  The  high  dose  of  the  laser  irradiation  (150  mJ)  and  reduced  protein  concentration  (up  to  0.12  mg/ml)  caused  a 
pronounced  loss  of  the  OGDC  activity  (see  Table  1).  Even  lower  doses  of  the  laser  irradiation  resulted  in  marked  inhibition 
of  the  diluted  preparation  of  the  complex. 


Table  1.  Activity  of  the  2-oxoglutarate  dehydrogenase  complex 


affected  by  different  doses  of  UV  laser  irradiation  (t=4°C) 


Enzyme  concentrations  (mg/ml) 

ACTIVITY,  % 

XeCl  laser  ( mJ) 

Nitroi 

?en  laser  (mJ) 

25 

50 

150 

240 

720 

1200 

20GDC,  0.7 

114 

108 

104 

100 

108 

60 

20GDC,  0.7  +  0.5  mM  DTT 

- 

128 

118 

108 

112 

100 

20GDC,  0.25 

156 

100 

- 

100 

100 

70 

20GDC,  0.12 

40 

40 

20 

- 

- 

- 

20GDC,  0.12  0.5  mM  DTT 

100 

128 

118 

- 

- 

- 

The  inhibition  of  the  complex  activity  after  the  treatment  of  the  enzyme  by  high  doses  of  both  the  lasers  used  in  the 
experiments  was  almost  abolished  when  a  reducing  agent  dithiotreitol  at  a  concentration  of  0.5  mM  was  presented  m  the 
OGDC  solution,  and  even  some  activation  was  noticed  after  the  complex  was  subjected  to  excimer  and  nitrogen  lasers 
irradiation.  The  effect  was  characteristic  both  for  more  concentrated  preparations  and  diluted  ones  .Interestingly,  the  low 
doses  of  XeCl  laser  irradiation  (25  mJ)  resulted  in  activation  of  the  OGDC  in  preparations  with  lower  (0.25  mg/ml)  protein 
concentration. 
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The  2  -oxoglutarate  dehydrogenase  complex  isolated  from  bovine  heart  possesses  two  types  of  active  sites  with  different 
substrate  affinity,  the  property  we  have  reported  earlier'"' .  We  therefore  ran  an  experiment  to  test  the  effect  of  the  type  of 
laser  irradiation  and  the  dose  applied  on  the  kinetics  of  20GDC  and  possible  changes  in  the  S0.5  and  V^ax  values  were 
determined.  The  low  doses  of  the  XeCl  laser  (25  mJ)  applied  for  diluted  preparations  of  the  complex  produced  a  3-fold 
decrease  of  the  S0.5  values  estimated  for  the  slowly  functioning  active  sites  localized  on  the  OGD  dimers  (from  250  to  33 
mM),  simultaneously  increasing  this  index  calculated  for  the  sites  with  high  affimty  to  2-oxoglutarate  (from  9  to  14  mM)  and 
producing  a  100%  increase  of  the  V^ax  of  the  catalysis  by  these  type  of  active  sites.  The  effects  obtained  on  the  concentrated 
20GDC  were  realized  by  the  mean  of  S0.5  values  decreasing  determined  for  quickly  -  reacting  sites  (20%)  without  noticeable 
changes  in  the  Vmax- 

The  inhibition  of  OGDC  activity  by  nitrogen  laser  irradiation  was  attended  by  7.5  -  fold  increase  of  the  S0.5  values 
determined  for  the  active  sites  binding  2  -  oxoglutarate  with  high  affinity  and  by  2  -  fold  decrease  of  V„ax-  Thus,  the  laser 
irradiation  affects  the  affinity  of  the  OGD  to  the  substrate  as  well  as  the  velocity  of  the  reaction  of  2  -  oxoglutarate 
decarboxylation. 

The  third  component  of  20GDC,  dihydrolipoamide  dehydrogenase,  catalysing  the  terminating  reaction  of  the  complex,  was 
more  sensitive  to  the  effects  of  irradiation.  The  treatment  of  E  3  by  high  doses  of  the  nitrogen  laser  (1200mJ)  inhibited  50% 
of  the  enzyme  activity.  The  excimer  laser  effects  on  the  activity  of  the  third  component  of  OGDC  depended  on  the  dose  of 
irradiation  applied,  and  the  most  significant  inhibition  (45%)  was  registered  in  response  to  the  protein  irradiation  by  308  nm 
laser  light  with  a  dose  of  1 50  mJ. 

It  is  suggested  that  the  radiation  -  caused  changes  of  multi  -  enzyme  complex  functional  activity  can  occur  at  least  partially 
as  a  result  of  changing  of  conformational  states  of  OGDC  protein  molecules  .  Intermolecular  association  -  dissociation 
processes  may  also  take  place.  The  mechanisms  can  be  realized  via  breaking  of  the  intersubunit  interactions  after  a  high  dose 
irradiation.  Radiation  damage  to  protein  can  also  occur  by  free  radical  reactions  as  a  result  of  free  radical  attacks  destroying 
subunits  of  the  enzymes  and  producing  unexpected  target  sizes  of  some  protein'^.  So  DTT  protecting  OGDC  from  a  loss  of 
activity  can  act  as  a  free  radical  scavenger. 

The  next  step  of  the  work  was  focused  on  the  study  of  the  laser  irradiation  influence  on  the  fluorescence  of  tryptophane 
residues  and  bound  flavine  dinucleotide  in  order  to  check  whether  20GDC  protein  molecules  underwent  some 
conformational  rearrangements  induced  by  the  irradiation. 

20GDC  was  shown  to  possess  significant  intrinsic  fluorescence  (the  main  parameters  of  it  are  listed  in  the  Table  2).  Excited 
at  295  nm,  in  the  region  of  tryptophane  residues  absorption,  the  protein  complex  exhibits  fluorescence  with  a  maximum  at 
330  nm.  This  short-wave  position  of  the  spectrum  may  be  a  result  of  tryptophane  residues  localization  in  low-  polar  regions 
of  the  protein  globule,  wftich  does  not  enable  water  to  penetrate.  The  other  possible  reason  is  a  sigmficant  slowing  of 
relaxation  processes  in  chromophores  after  excitation.  The  polarization  degree  P,  determined  by  segmental  mobility  of 
aminn  add  residues,  is  characterized  by  low  values  (around  0.12),  confirming  our  first  suggestion.  Obtained  small  long¬ 
wave  shift  of  spectrum  (~2  nm,  see  Table  2)  when  using  “red  edge  excitation”  indicates  microenvironment  heterogeneity  of 
the  excited  tryptophanyls. 


Table  2.  Parameters  of  2-OGDC  intrinsic  fluorescence.  ^*=295  nm.  Protein  concentration  was  0.3  mg/ml. 


Sample 

Intensity 

^Tnax5 

P 

cp^> 

2-OGDC 

1.00 

330 

0.12 

0.050 

2-  OGDC  " 

- 

332 

0.13 

- 

2-  OGDC^’ 

0.90 

330 

0.12 

0.045 

^ex  =  305  nm 

after  the  XeCl  laser  irradiation  (150  mJ) 

quantum  yield  was  determined  relatively  to  tryptophane  with  (p=0.20. 


The  fluorescence  decay  was  nonexponential  (Table  3,  Fig.l)  and  might  be  described  by  three  exponential  terms  with 
fluorescence  lifetime  t,=  0.89  ns,  T2=  2.6  ns,  Z3=  5.4  ns  under  conditions  of  satisfactory  values  for  statistical  criteria. 


Figure  1.  Fluorescence  decay  curve  of  20GDC. 

A-  experimental  and  calculated  decay  curve:  dots 
correspond  to  experimental  data,  solid  curve  is  the 
result  of  approximation  with  three-exponential 
function;  B-  distribution  of  weighted  residuals; 

C-  autocorrelation  function  of  weighted  residuals. 


Figure  2.  Fluorescence  spectra  of  native  20GDC  (1) 
and  after  the  XeCl  laser  irradiation  using  the  doses  of 
25  (2)  and  150mJ(3). 


Table  3.  Fluorescence  decay  lifetimes  of  OGDC  approximated  by  single,  bi-,  and  three- 
exponential  functions.  297  nm,  340  ran.  Standard  deviation  of  t  determination  did  not 
exceed  5%.  (t)=  aiTi  +a2'C2+a3T3,  (a)+a2+a3=l) 


a. 

-Cnns 

0.2 

T2,  ns 

02 

T3,ns 

<x>,  ns 

1.00 

3.30 

- 

- 

- 

- 

- 

13.8 

0.37 

1.20 

0.63 

4.60 

- 

- 

- 

1.04 

0.48 

0.88 

0.31 

2.64 

0.21 

5.38 

2.40 

0.97 

Nonexponential  character  of  tryptophane  fluorescence  decay  is  a  typical  property  of  proteins'^.  Most  of  the  studied  proteins 
demonstrate  the  same  character  of  the  intrinsic  fluorescence  decay  presented  as  a  sum  of  three  exponentials.  There  may  be 
three  equilibrated  conformations  of  tryptophanyls  in  respect  to  a  peptide  bond  where  they  possess  different  values  for  the 
fluorescence  decay  constants.  The  protein  structure  is  not  static  and  undergoes  significant  fluctuations'®.  As  a  result,  the 
function  of  decay  is  occasionally  presented  by  three  continuous  chromophore  distribution  according  to  the  constants  of 
emissive  transition'^.  Such  distributions  are  fi-equently  wide  enough  and  are  overlapped,  and  only  one  or  two  chromophore 
distributions  can  be  determined  using  mathematical  calculations'*.  Solving  the  problem  of  the  distribution  type  is  a 
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complicated  task  and  need  further  experiments.  The  data  obtained  from  the  fluorescence  decay  curves  confirm  the 
assumption  about  high  heterogeneity  of  microenvironment  of  emitting  tryptophanyls. 

20GDC  comprises  the  incorporated  FAD  as  a  coenzyme  -wWch  possesses  the  property  of  intrinsic  fluorescence  with  a 
maviTniim  of  the  spectTum  at  530  nm  (Fig.2).  The  fluorescence  intensity  of  the  native  protein  is  low  and  is  equal  only  to  1% 
of  the  free  coenzyme  fluorescence  recorded  under  the  same  conditions. 

The  XeCl  laser  irradiation  of  the  protein  resulted  in  some  decrease  of  tryptophane  fluorescence  of  the  2-oxoglutarate 
dehydrogenase  complex  (Table  2)  producing  no  changes  in  the  emission  spectrum,  its  localization  and  polarization  degree. 
The  FAD  fluorescence  was  more  sensitive  to  the  irradiation.  The  FAD  emission  intensity  increased  twofold  when  treated  by 
a  high  dose  (150  mJ)  of  the  excimer  laser  (Fig.2). 

The  data  obtained  indicate  that  the  active  site  localized  on  the  E3  component  of  20GDC  had  a  pronounced  sensitivity  to  the 
XeCl  laser  irradiation.  These  results  are  in  agreement  with  the  kinetic  studies  vihich  showed  a  loss  of  this  constituent  enzyme 
activity  depending  on  the  irradiation  dose  applied.  It  is  known  that  the  fluorescence  of  FAD  incorporated  into  the  protein  is 
effectively  quenched  by  the  neighboring  amino  acid  residues  (tyrosine,  tryptophane  or  cysteine)  possessing  groups  which 
can  be  electron  donors'®.  The  laser  irradiation  causes  conformational  changes  in  the  enzyme  structure,  thus  affecting  the 
structure  of  the  active  site.  These  processes  resulted  in  decreased  quenching  and  enhanced  fluorescence  of  the  protein 
molecule.  Moreover,  one  or  several  tryptophane  residues  possibly  localized  in  the  vicinity  of  the  active  site  contribute  to  the 
changes  in  the  protein  intrinsic  fluorescence.  Since  20GDC  is  a  polytryptophane  protein,  we  have  not  foimd  noticeable 
changes  in  a  number  of  fluorescence  parameters,  such  as  spectrum  localization,  polarization  degree  and  others. 
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Abstract 

Because  of  low  operating  speed  and  excessive  collateral  damage, 
lasers  have  not  succeeded  in  replacing  conventional  tools  in  many 
surgical  and  dental  applications.  Recent  developments  now  allow 
the  new  generation  of  cimplified  ultrashort  pulse  lasers  to  operate 
at  high  repetition  rates  and  high  single  pulse  energies.  A 
Titanium: sapphire  Chirped  Pulse  Regenerative  Amplifier  system 
operating  at  1  KHz  and  50  fs  pulse  duration,  was  used  to 
demonstrate  ultrashort  pulse  ablation  of  hard  and  soft  tissue. 
Maximum  ablation  rates  for  enamel  and  dentin  were  approximately 
0.650  pm  /pulse  and  1.2  pm  /  pulse  respectively.  Temperature 
measurements  at  both  front  and  rear  surface  of  a  1  mm  dentin  and 
enamel  slices  showed  minimal  increases.  Scanning  electron 
micrographs  clearly  show  that  little  thermal  damage  is  generate  by 
the  laser  system.  If  an  effective  delivery  system  is  developed, 
ultrashort  pulse  system  may  offer  a  viable  alternative  as  a  safe, 
low  noise  dental  tool. 

Introduction. 

Using  lasers  to  ablate  hard  tissue  offers  the  potential  for  rapid, 
precise,  accurate  operation  with  minimal  thermal  and  mechanical 
damage  to  surrounding  tissue.  Lasers  have  therefore  been  used  for 
a  variety  of  sensitive  hard  and  soft  tissue  applications  such  as 
mid-ear  bone  ablation,  implantation  of  joint  prostheses  in 
orthopedic  surgery,  removal  of  malignant  tissue  and  other 
procedures  involving  bone,  nail,  cartilage  and  other  hard  and  soft 
tissue. 
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In  the  field  of  dentistry,  laser  applications  include  treatment  of 
oral  malignancies  and  periodontal  diseases,  dental  caries 
prevention,  control,  and  removal  of  hard  tissue  for  cavity  and 
restoration,  preparation,  and  sterilization  (1  -  9).  Clearly,  the 
greatest  anxiety  generators  in  the  dentist's  office  are  the  dental 
drill  and  the  accompanying  anesthesia  shots  that  are  required  for 
most  hard  tissue  procedures.  Indeed,  the  fear  of  pain  during 
dental  procedures  is  so  great  that  it  deters  patients  from  seeking 
help  while  the  dcimage  to  healthy  tooth  structure  continues  to 
spread  (10,11).  According  to  a  survey  by  Smith  et  al  (10),  40%  of 
the  population  reported  that  they  have  fear  of  dental  treatment, 
with  17%  reporting  high  anxiety. 

In  a  recent  review  ( 9 )  Wigdor  et  al  reported  that  many  patients 
are  so  fearful  of  dental  procedures  that  they  frequently  wait 
until  the  pain  from  infection  is  so  great  that  they  are  willing  to 
endure  the  perceived  painful  treatment  to  gain  relief.  The  long 
tern  goal  of  dental  laser  research  has  thus  long  been  the 
development  of  a  new  modality  to  replace  the  dental  hand  piece. 

In  spite  of  their  promise,  lasers  has  not  been  able  to  replace  the 
drill  or  scalpel  in  many  surgical  and  dental  procedures.  The 
limitation  of  existing  laser  systems  will  be  outlined  below, 
essentially,  however,  there  are  two  interaction  characteristics 
which  limit  the  success  of  surgical  and  dental  procedures;  low 
removal  rates  and  high  collateral  damage. 

Low  removal  rates  is  generally  a  characteristic  of  short  pulse 
lasers  (-ns).  It  typically  involves  material  removal  on  the  order 
of  a  single  to  a  few  pm  in  depth  and  low  residual  heating.  Since 
pulse  duration  is  relatively  short,  thermal  seed  electron  induce 
plasma  generation  through  subsequent  avalanche  process.  The 
plasma  mediated  ablation,  results  in  plasma  shielding  as  pulse 
energy  is  increased  which  followed  by  a  sharp  drop  in  per-pulse 
ablation  efficiency  with  energy.  As  a  consequence,  all  high 
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intensity  short  pulse  lasers  are  limited  in  their  ablation  rates 
to  the  pm/pulse  range. 

Longer  pulse  ablation  such  as  the  100  to  300  ps  pulsed  free- 
running  normal  mode  solid  state  lasers  results  in  ablation  on  the 
order  of  several  tens  of  a  micrometer  per  pulse.  Ablation  with 
these  systems,  however,  is  also  characterized  by  a  considerably 
larger  residual  pulse  energy,  larger  amount  of  material  ejection 
and  a  larger  stress  transient. 

It  is  thus  possible  to  generalize  and  assert  that  under  the 
current  status  of  laser  technology  (where  ultrashort  pulse  lasers 
are  not  considered),  we  have  a  choice  between  a  lower  ablation 
rate,  higher  precision  nanosecond  systems  with  lower  thermal  and 
mechanical  signature,  and  that  of  the  considerably  higher  ablation 
rate  longer  pulse  system  also  characterized  by  higher  residual 
heat  and  large  transient  stresses. 

When  a  practical  system  is  considered  several  concerns  must  be 
addressed  simultaneously:  Tissue  selectivity,  ablation 
efficiency,  thermal  and  mechnical  damage,  confinement  of  ablative 
effects,  ease  of  delivery,  ablated  surface  morphology, 
microhardness  and  changes  in  mineralization.  The  perfect  laser 
tool,  of  course,  need  to  address  and  give  adequate  solutions  to 
all  of  these  problems,  since  any  one  of  these  potential  problem 
could  jeopardize  the  tooth  or  result  in  an  inferior  treatment.  In 
addition,  an  adequate  practical  tool  must  be  sufficiently 
inexpensive  to  be  affordable  to  private  practitioners. 

Previous  studies  of  ultrashort  pulse  regime  with  hard  dental 
tissue  included  pulse  durations  from  the  100  fs  to  the  ns  range 
(12,  13,  14,  15).  This  paper,  describes  recent  experiments  using 
USPL  systems  operating  at  very  high  peak  power.  It  demonstrates 
the  technical  advantage  of  USPL  over  other,  existing  laser  systems 
as  well  as  conventional,  non-laser  technologies.  The  problems 
presented  by  the  current  high  cost  of  USPL  system,  and  the 


relatively  large  size  of  current  USPL  systems  must  naturally  be 
addressed. 

Material  and  Methods 

Ultrashort  pulse  Laser  system. 

For  plasma  mediated  ablation  studies,  laser  pulses  generated  by  a 
Ti: Sapphire  Chirped  Pulse  Amplifier  (CPA)  system  with  an 
frequency-doubled,  Nd:YAG  pump  laser.  The  800  nm  Ti: Sapphire 
oscillator  was  pumped  by  a  Millennia  Laser  (Spectra  Physics,  CA) 
used  also  pumped  by  an  Nd:YAG  source.  Output  compressed  pulses, 

50  fs  long  and  up  to  500  pJ  in  energy  at  1  KHz  were  used.  Various 
attenuator  allowed  determination  of  threshold  fluence  and  ablation 
characteristics  over  the  entire  energy  range.  A  spot  size  of  100 
pm  allowed  studies  of  fluences  ranging  from  0.4  J/cm^  to  4  J/cm^. 


Grating  Stretcher 


which  are  stretched  using  a  1200  1/mm  grating  pair.  The  pulses  are 
amplified  to  1  mJ  in  the  regenerative  amplifier  at  kilohertz 
repetition  rates.  The  pulses  are  recompressed  by  a  1200  1/mm 
grating  pair  which  results  in  50  fs,  500  pJ  pulses. 


Ablation  studies  were  performed  using  thin  dentin  slices.  Dentin 
slices  0.5  to  1.0  mm  thick  were  cut  from  the  middle  section  of 
freshly  extracted  third  molars,  in  a  plane  perpendicular  to  the 
occlusal-cervical  direction.  Teeth  were  treated  with  0.5M  EDTA 
for  2  minutes  to  remove  the  smear  layer,  then  stored  in  10%  thymol 
solution  until  treatment.  Tissue  was  sliced  parallel  to  the  crown 
and  washed  with  demineralized  water.  Slices  included  an  external 
ring  of  enamel  tissue  on  the  outside  and  dentin. 

Ablation  rates  were  determined  by  monitoring  the  ablated  slice 
with  a  video  camera.  A  frame  by  frame  analysis  allowed  detection 
of  ablation  start  point  as  well  as  that  of  the  perforation 
endpoint . 

A  thermal  camera  (Inframetrics  Model  600  Imaging  Radiometer, 
Bedford,  MA)  with  a  scan  speed  of  32  frames  per  second  was  used 
for  temperature  measurements.  The  center  temperature  adjustment 
ranges  from  20  to  1,500  “c.  The  camera  field  of  view  could  be 
adjusted  from  1:1  to  8:1,  and  was  maintained  at  2:1  or  4:1.  The 
camera  detector  was  a  liquid  nitrogen  cooled  HgCdTe  with  spectral 
sensitivity  from  8  to  12pm.  The  thermal  camera  was  placed  so  that 
the  surface  temperature  (ST)  could  be  viewed  directly.  The  camera 
was  operated  in  the  point  mode  so  that  both  the  entire  temperature 
field  could  be  viewed  with  false  color  presentation  and  a  single 
point  (specified  by  the  operator)  numerical  value  could  be  read 
directly  with  the  camera  specified  accuracy  of  ±  0.5  degrees. 
Temperatures  were  measure  on  the  back  surface  of  the  1  mm  dentin 
slices  while  the  laser  was  ablating  the  front  surface.  Camera 
emissivity  was  calibrated  by  conparing  temperature  of  dentin  slice 
half  painted  with  black  ink.  The  emissivity  of  specification  for 
the  camera  software  was  adjusted  until  the  temperature  reading  of 
the  non  painted  section  matched  that  of  section  painted  with  the 
black  ink.  The  camera  was  further  checked  against  a  standard 
thermometer  readings  of  water  reservoir  which  temperature  was 
gradually  varied. 


Results 
Ablation  rates 

The  ablation  rates  in  enamel  and  dentin  at  all  fluence  levels 
studied  are  summarized  in  figure  2  through  6 .  The  threshold  for 
plasma  formation  for  both  enamel  and  dentin  was  0.4  J/^m  As  is 
evident  from  the  figure  2  and  figure  3,  both  enamel  and  dentin 
show  a  clear  jDattern  of  diminishipa  amount  of  material  remove^  oer 
pulse  as  the  pulse  fluence  increases.  From  the  ablation  threshold 
at  approximately  0.5  J/cm^,  ablation  rate  (AR)  increased  rapidly  to 
about  1  pm/pulse  at  3  J/cvc?  where  ablation  for  both  tissue  types 
were  very  similar.  While  the  ablation  threshold  appears  higher 
for  interaction  with  enamel,  this  is  really  a  consequence  of  our 
method  of  measuring  ablation  depth  per  pulse  (slice  perforation) 
and  the  stalling  that  occurs  earlier  in  enamel  unless  higher 
fluences  are  used.  Beyond  that  level,  only  a  very  small  increase 
occurs  in  AR  with  fluence 

For  comparison,  figure  4  illustrate  ablation  with  higher  pulse 
energies  and  fluences  but  longer  pulse  duration,  namely  that  of 
350  fs  (and  1.054  pm)  and  one  nanosecond.  As  can  be  seen  from  the 
figure,  the  ablation  pattern  for  the  350  fs  pulses  is  similar  to 
that  of  the  50  fs  one.  Once  the  pulse  duration  crosses  the 
ultrashort  threshold  (-  10  ps)  the  ablation  fluence  threshold  for 
the  1.054  pm  1  ns  pulses  is  considerably  higher  (20  J/cm^),  and  the 
ablation  characteristics  change. 

As  figure  4  shows,  the  nanosecond  pulse  ablation  at  34  J/cm^ 
illustrated  the  AR  differences.  Here,  3  J/cm^  is  well  below  the 
ablation  threshold  of  either  enamel  or  dentin.  Ablation  threshold 
for  the  nanosecond  system  was  obse3rved  at  approximately  20  J/cm^ 
but  is  not  shown  here  since  it  was  not  the  object  of  this  study. 

For  the  same  fluence  level  in  the  nanosecond  regime  AR  are  very 
different  for  enamel  and  dentin  (dentin  ablation  is  almost  a 
factor  of  four  greater) .  The  comparison  thus  shows  the 
considerably  enhance  efficiency  of  the  CPA  pulses  (approximately  1 
pm/mJ  at  1  J/cm^  for  both  enamel  and  dentin)  in  comparison  to  the 


34  J/cm^  Ins  pulse  with  0.03  pm/mJ  for  ablation  of  Enamel,  and  0.12 
pm/mJ  of  dentin.  It  also  shows  the  lack  of  material  sensitivity 
of  the  fs  ablation,  especially  when  compared  to  the  ns  regime. 


FLUENCE  (J/cm^) 


Figure  2. 

Dentin  ablation  rates  as  a 
function  of  50  fs  pulse  fluence. 


Figure  3 . 

Enamel  ablation  rates  as  a 
function  of  50  fs  pulse  fluence 


Figure  5  and  Figure  6  show  the  dependence  of  the  per-pulse 
ablation  rate  on  the  sample  thickness.  As  was  discussed 
previously  (16)  the  ablation  spot  size  changes  as  a  function  of 
penetration  depth  into  the  tissue.  To  investigate  if  that  change 
also  translated  into  a  per-pulse  ablation  rate  changes  various 
slice  thickness  were  tested. 


As  both  figure  5  and  6  show,  essentially  no  statistically 
significant  variation  were  noted  in  either  enamel  or  dentin.  On 
the  other  hand,  at  the  maximum  pulse  energy  of  400  pJ  per  pulse 
the  largest  enamel  slice  that  was  penetrated  was  380  pm  thick, 
while  the  largest  dentin  slice  that  was  penetrated  was  900  pm 
thick. 


in 


Figure  4  Dentin  (triangles)  and  enamel  (squares)  ablation  rates 
as  function  of  fluence  for  350  fs  pulses  (empty  symbols)  and  1  ns 
pulses  ( full  symbols ) . 
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Figure  5 

Dentin  ablation  rates  as  a 
function  of  ablated  slice 
thickness 


Figure  6. 

Enamel  ablation  rates  as  a 
function  of  ablated  slice 
thickness 


Temperature  measurements . 

Temperature  observed  at  the  back  surface  of  a  1mm  thick  dentin 
showed  an  increase  of  up  to  45  °C  when  the  beam  was  applied 
continuously  for  a  period  of  several  minutes.  This  increase, 
however,  was  observed  only  when  the  laser  beam  ablation  penetrate 
the  1mm  dentin  slice  and  the  interaction  was  taking  place  very 
close  to  the  back  (i.e.  monitored)  surface.  When  the  ablation  was 
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confined  to  interactions  at  the  front  surface  alone,  temperature 
increases  at  a  surface  1  mm  away  were  confined  to  only  5  to  10  ^C. 

Scanning  electron  microscopy 

The  ablation  fluence  of  3  J/cm2  exhibits  surface  characteristics 
obseirved  in  slightly  longer  pulses  of  350  fs  (15).  The  lack  of 
crumbling  and  clean  crater  edge  illustrate  this  point  and  the  high 
localization  of  the  interaction  processes  with  ultrashort  pulse. 
The  extremely  smooth  damage-free  characteristics  of  the  crater 
walls  is  clearly  evident.  Ablation  crater  walls  demonstrated 
small  granules  a  few  micrometer  in  size.  Exposed  dentinal  tubules 
and  the  lack  of  significant  amount  of  melted  dentin  clearly 
indicate  a  lack  of  significant  heating  and  the  mechanical  nature 
of  the  ablation  process  with  fs  pulses. 

Enamel  ablation  is  similar  in  appearance  to  that  of  dentin.  The 
main  difference  is  that  in  the  range  of  0.5  to  5  J/cm^  dentin 
ablation  results  in  larger  deeper  craters  (see  discussion  of 
ablation  rate  below) .  Morphologically,  ablated  surface  are  very 
smooth,  free  of  thermal  damage  or  cracking. 

Conclusion 

Negligible  collateral  damage  and  highly  efficient  ablation  rate 
per  pulse  in  the  ultrashort  pulse  regime  was  demonstrated.  A 
comparison  of  50  fs  and  350  fs  to  pulses  1  ns  long,  clearly 
demonstrated  the  ultrashort  pulse  regime  reduced  sensitivity  to 
material  composition.  The  strong  dependence  of  ultrashort  pulse 
multi-photon  interaction  processes  on  threshold  fluence,  along 
with  scanning  electron  micrographs  illustrate  a  sharply  defined 
ablation  threshold  and  the  highly  precise  spatial  nature  of  the 
interaction.  Analysis  of  our  temperature  increases  during  the 
interaction  confirms  expectations  of  relatively  small  temperature 
increases  and  confinement  of  thesnnal  effects.  Thus,  the  potential 
for  an  efficient,  precise,  and  minimally  damaged  ultrashort  pulse 
operation  was  demonstrated,  offering  hope  for  practical 
application  of  such  devices  in  the  near  future. 
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with  Ultrashort-Pulse  Lasers 


Optical  parameters  measurements  by  collimated  light  transmission 
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Maria  Quarto,  Giuseppe  Roberti,  Fulvia  Villani 

Faculty  of  Medicine  and  Surgery,  University  of  Naples  “Federico  11” 

ABSTRACT 

Some  medical  ^plications  involving  the  interaction  between  light  and  biological  tissues  require  both  the  knowledge  of 
optical  characteristics  of  tissues  and  a  realistic  treatment  of  light  transport  into  them.  In  this  work  we  describe  a 
transillumination  technique  testing  water  solutions  of  Intralipid  by  the  transmitted  radiation  intensity.  The  experimental 
apparatus  includes  a  diode  laser  (Xemiss  ==  820  nm),  a  detection  fiber,  a  PMT  and  a  Digital  Signal  Analyzer. 

The  performances  of  this  very  simple  and  cheap  system  are  comparable  to  those  ones  obtained  with  more  sophisticated 
apparatus;  the  results  show  that  this  technique  could  represent  a  preliminary  step  toward  the  realization  of  an  user-fiiendly 
and  cheap  laser  system  for  measuring  optical  parameters  of  tissues. 

Keywords:  transillumination.  Intralipid 


1.  INTRODUCTION 

The  research  on  the  breast  cancer  by  transillumination  represents  one  of  the  most  promising  diagnostic  application  of  light. 
The  comprehension  of  phenomena  that  govern  the  interaction  between  light  and  biological  tissues  has  an  important  part  in 
the  use  of  transillumination. 

To  describe  the  problem  of  light  propagation  in  biological  tissues,  two  different  theories  have  been  formulated:  the 
analytical  theory  and  the  transport  theory.  The  analytical  fiieory,  that  bases  on  the  Maxwell’s  equations,  represents  the  most 
rigorous  approach  allowing,  as  a  rule,  to  take  into  account  the  effects  connected  to  the  multiple  scattering,  the  dif&action  and 
the  interference.  However,  the  unwareness  of  dielectric  characteristics  of  tissues  and  the  mathematical  difficulties  make  it 
useless  in  the  practice.  In  the  transport  theory  the  energy  transported  by  light  is  highly  considered  in  spite  of  light  wave 
nature.  In  this  case,  the  optical  behaviour  of  biological  tissues  is  fixed  by  the  scattering  coefficient  fis»  fb®  absorption 
coefficient  pa  and  tb®  phase  function  p(0)  [1].  In  practice,  a  general  solution  of  the  transport  equation  isn’t  available  yet.  It  is 
demonstrated  that,  in  particular  conditions,  it  is  possible  to  get  rough  but  simple  enough  analytical  solution.  In  particular,  the 
diffusion  approximation  is  used  to  define  the  distribution  of  light  intensity  in  a  medium  where  the  diffusion  is  dominant;  this 
happens  for  the  biological  tissues  in  the  so-called  “therapeutic  window”,  700-1300  nm  [2,3].  In  fact,  in  this  spectral  range, 
where  the  largest  penetration  of  light  in  the  biological  tissues  (about  10  cm)  is  obtained,  the  scattering  coefficient  is  much 
higher  than  the  absorption  one  [3];  on  this  account  we  can  assert  the  light  attenuation  in  this  spectral  region  is  mainly  due  to 
the  scattering  that  origins  fi’om  the  changes  of  the  refiactive  index  of  coimective  tissue  and  cellular  costituents,  as 
membrans,  cell  organelles  and  nuclei  [4]. 

In  this  work,  we  present  measurements  of  spatial  distribution  of  transmitted  light  through  samples  irradiated  at  820  nm.  The 
technique  presented  in  this  work,  although  very  simple,  could  be  used  as  an  alternative  to  the  time-resolved 
transillumination.  In  fact  the  light  detection  is  performed  with  an  optical  fiber  whose  numerical  aperture  allows  to  select  the 
less  diffiised  photons.  The  effect  of  such  a  detection  could  be  considered  as  that  one  obtained  by  reducing  the  time-gate  in 
time-resolved  techniques  and  allows  to  detect  only  the  photons  that  covered  the  shortest  paths. 

2.  MATERIAL  AND  METHOD 

The  experimental  apparatus  used  for  the  trasmittance  measurements  is  represented  in  Fig.  1.  A  continous  laser  beam  emitted 
by  a  diode  laser  (Hamamatsu  PLP02,  X  =  820  nm,  CW  =  O.lmW),  after  having  been  collimated  thanks  to  an  aspheric  lens 
(NA  =  0.50,  effective  focal  =  2.00),  hits  the  sample  placed  in  a  fixed  position  on  a  sample-holder. 

The  transmitted  radiation  was  collected  by  a  multimodal  optical  fiber  (diameter  core  =  600  jJm,  NA  =  0.4)  whose  head  was 
moved  in  a  vertical  plane  ortogonal  to  the  laser  beam  direction  by  two  computer  controlled  micrometric  traslators  (Physik 
Instrument,  travel  range  5  cm,  risolution  1  pm).  The  motion  was  performed  along  the  vertical  direction  at  a  value  of  the 
horizontal  coordinate  corresponding  to  the  center  of  the  cuvette,  for  a  travel  equal  to  3  cm  at  600  pm/s  speed.  The  detection 
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has  been  carried  out  by  a  PMT  (Hamamatsu  R406)  and  a  Digital  Signal  Analyser  (Tektronix  DSA602,  sampling  rate  -  2 
Gsample/s,  bandwidth  =  1  GHz). 


diode 

laser 

cd±b  aspherte  lens 


Fig.  1  Experimental  set-up. 


The  samples  have  been  obtained  mixing  suitable  quantities  of  distilled  water  and  Intralipid  10%  (Farmacia,  Italy)  and  placed 
in  plastic  cuvette  (1  cm  thick).  The  power  measured  on  each  sample  has  been  normalized  to  the  power  transmitted  through 
the  same  cuvette  filled  with  distilled  water.  Since  the  refiactive  index  of  Intralipid  solutions  is  very  close  to  that  one  of  the 


water,  the  normalization  allows  to  disregard  the  reflections  at  the  glass-air  and  glass-sample  interfaces  due  to  the 
discontinuity  of  the  refiactive  index.  To  obtain  the  scattering  and  absorption  coefficients  of  the  samples,  the  spatial  profile  of 
trasmitted  intensity  has  been  compared  with  the  solutions  of  the  dif^ion  equation  in  case  the  sample  is  irradiated  by  a 
collimated  beam  with  a  small  diameter  and  gaussian  spatial  profile  [5]. 

The  Intralipid  attenuation  coefficient  pp  h^s  been  determined  thanks  to  the  Beer’s  law  by  absorption  measurements  on 
diluted  solutions  placed  in  0.5  cm  thick  cuvette  to  reduce  multiple  scattering  effects.  On  the  other  hand,  since  Intralipid  is  a 
highly  scattering  medium  (pa  is  several  ordersof  magnitude  smaller  than  ps),  the  pt  value  could  be  approximed  to  the  pg 
value.  The  value  obtained  for  the  scattering  coefficient  is  pg  =  (227  +  11)  cm  .  As  far  as  the  absorption  coefficient  is 
concerned,  we  choose  the  value  equal  to  0.003  cm  reported  in  literature  [6].  Finally,  for  the  anisotropy  factor,  we  choose 
the  value  of  0.75,  obtained  by  the  extrapolation  of  literature  data  at  800  nm  [6,7,8]. 


3.  RESULTS  AND  DISCUSSION 

The  measurements  of  the  power  of  the  light  transmitted  by  the  analyzed  water  solutions  of  Intralipid  as  a  function  of  the 
coordinate  r  which  represents  the  vertical  position  of  the  fiber  with  resjrect  to  the  intersection  between  the  direction  of  laser 
beam  and  the  scanning  plane  are  reported  in  Fig.  2.  The  data  are  normalized  as  described  above.  The  plots  refer  to  four 
Intralipid-water  concentration  ratio  (1:10,  1:15,  1:20  and  1:40).  The  solid  line  superimposed  on  the  experimental  data 
represents  the  solution  of  the  diffusion  equation  obtained  by  using  the  pg  and  pa  values  obtained  by  the  linear  scaling  of  the 
vdues  reported  in  the  previous  paragraph  as  a  function  of  the  concentration.  The  plots  show  a  good  agreement  between  the 
experimental  and  theoretical  data.  The  normalized  light  power  behaviour  shows  that  the  beam  is  highly  attenuated  mainly 
owing  to  the  diffusion  that  spreads  the  photons  away  fixrm  the  input  direction,  the  beam  broadening  decreasing  with 
decreasing  Intralipid  concentration. 
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withdrawn  or  are  unavailable. 

[3255-02]  Transparent  media 

A.  L.  Gaeta,  J.  Ranka,  Cornell  Univ. 

[3255-04]  Pressure  and  temperature  evolution  induced  by  ultrashort-laser-pulse  ablation 
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